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INTRODUCTION 

I 

This document describes the Langley Research Center 16- Foot Transonic 
Tunnel standard set of equations. The engineering units necessary for these 
equations are computed on site from the raw data millivolts or counts. These 
quantities with additional constants are used as input to the program for 
computing the forces and moments and the various coefficients. 

This document supersedes NASA Technical Memorandum 86319, 
Computations for the 16-Foot Transonic Tunnel, Revision 1, January, 1987. 

This document is intended to be a companion document to NASA Technical 
Memorandum 102750, A User’s Guide to the Langley 16-Foot Transonic Tunnel, 
Revision 1, September 1990. 

The equations are grouped into modules, so that only the required modules 
need be used. The modules are as follows: 



A. Wind Tunnel Parameters 

B. Jet Exhaust Measurements 

C. Skin Friction Drag 

D. Balance Loads and Model Attitudes 

E. Internal Drag (or Exit-Flow Distributions) 

F. Pressure Coefficients and Integrated Forces 

G . Thrust Removal Options 

H. Turboprop Options 

I . Inlet Distortion 

Individual customizing of these equations for a specific job application is 
permitted through the use of code constants. These equations do not cover all 
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possible jobs; however, they are coded so that modifications of selected equations 
may be easily carried out. 

The format of this document is arranged so that the module designations 
correspond to the Appendix designations in which the respective calculations 
equations are given. 


WIND TUNNEL PARAMETERS 

The wind tunnel parameters are computed from the required static and 
total pressure measurements. The Reynolds number, dynamic pressure and 
tunnel total temperatures are computed. When the tunnel Mach number is 
computed, a polynomial fit from the 1990 wind tunnel calibration is used to correct 
the ratio of static pressure to total pressure used in the Mach number calculation. 
These wind tunnel parameters are stored for use by other modules. Refer to 
Appendix A for calculations. 

JET EXHAUST MEASUREMENTS 

Jet exhaust information is calculated for the primary, secondary and 
tertiary flow conditions. 

The primary flow conditions for each engine, up to a maximum of four, are 
calculated. The various parameters that are computed are mass flow and ideal 
thrust for each engine. The average nozzle pressure ratio and average total 
temperature over all the engines is obtained. The total mass flow is derived from 
chamber, flowmeter, and/or venturi measurements. Discharge coefficients for the 
total system are computed as well as the ideal thrust. For the primary flows, a 
dual air supply system is used in providing inputs for the mass flow parameters. 

For the secondary and tertiary flows, the mass flows and other parameters 
are computed. Refer to Appendix B for calculations. 


2 


SKIN FRICTION DRAG 

The skin friction drag for the model is computed in addition to any 
empennage skin friction drag. Refer to Appendix C for calculations. Information 
from the wind tunnel parameters is used. Drag from the various components as 
well as total drag is computed. 

BALANCE LOAD AND MODEL ATTITUDES 
The balance computations for the force and moment coefficients for up to 
five balances may be computed from this module. Allowances for the method of 
attaching the balances are made. The measured forces and moments are 
corrected for balance interactions. Then an allowance is made for high onler 
interactions and momentum tares. The forces and moments are rotated to the 
desired axis and the final correct coefficients are computed as well as the angle of 
attack and sideslip angles. Refer to Appendix D for calculations. 

INTERNAL DRAG 

The internal drag and various forces on the engines are computed using the 
equations given in Appendix E. The result of these computations are used in the 
balance computations of module D to correct the force measured by the balances. 

PRESSURE COEFFICIENTS AND INTEGRATED FORCES 
Pressure coefficients are computed by using the equations given in 
Appendix F. Various integrated forces due to the pressures are calculated 
including hinge moment coefficients. 
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THRUST REMOVAL 

Various thrust removal coefficients may be computed according to specified 
flags which specify the model setup. Various configurations are permitted which 
may include two balances. Reference Appendix G for calculations. 

TURBOPROP OPTIONS 

The drag and thrust coefficients due to the propeller and jet engine are 
computed as well as the combined totals. Horsepower and efficiency of the engines 
are derived with other quantities. Reference Appendix H for calculations. 

INLET DISTORTION 

Inlet engine face pressure distortion and mass flow rates are computed by 
using the equations given in Appendix I. Various profiles or gradient across the 
engine face which result from the airstream entering into the inlet are determined 
from the measured pressures and calculated ratios. 
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APPENDIX A 
Tunnel Parameters 


Nomenclature 

Required Constants 

Atmospheric Pressure 

Mach Number 

Tunnel Static Pressure 

Tunnel Total Pressure 

Tunnel Dynamic Pressure.... 

Dew Point 

Tunnel Total Temperature 
Reynolds Number 


..7 

..9 

..9 

10 

10 

11 

11 

12 

12 
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MODULE A TUNNEL PARAMETERS 


SYMBOL 

MACH 

MCODE 


PO 

PO/PTO 

PTANKG 

PTANKH 

PTG 

PTH 

PTKSON 

PTO 

PTSON 

QO 

REFL 

RN 

RN/FT 

RT(J) 


NOMENCLATURE 
Free stream Mach number. 

Mach number calculation code. 

■1, PTANKG and PTH are needed. 

=2, PTANKH and PTH are needed. 

=3, PTANKG and PTG are needed. 

=4, PTANKH and PTG are needed. 

=5, PTKSON and PTSON are needed. 

Tunnel static pressure, Ibs/sq. in. 

Ratio of tunnel static pressure to total pressure. 
Tunnel tank pressure measured by gage, lbs/sq. in. 
Tunnel tank pressure measured by Ruska, lbs/sq. in. 
Tunnel total pressure measured by gage, lbs/sq. in. 
Tunnel total pressure measured by Ruska, lbs/sq. in. 
Tunnel tank pressure measured by Digiquartz, 
lbsVsq. in. 

Tunnel total pressure, lbs/sq. in. 

Tunnel total pressure measured by sonar manometer, 
lbs/sq. in. 

Dynamic pressure, lbs/sq. in. 

Reference length, feet. 

Reynolds number based on reference length. 

Reynolds number per foot. 

Tunnel total temperature measurements, *F, where 
J = probe number. 
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SYMBOL 

T(J) 


TTO 


NOMENCLATURE 

Constants required from project engineer, where 
J = probe number. If probe is bad or does not exist, then 
its value should be set to 0.0. If no correction is made to 
the temperature probe, then its value should be set to 1.0. 
Tunnel total temperature, *F. 


APPENDIX A 

Module A 

Tunnel Parameters 


A. Required Constant* 

1. MCODE (default value = 2) must be provided if values other than PTANKH 
and PTH are used to compute Mach number. 

2. The constants used in determining tunnel total temperature are T2, T3, 

T4 and T5 which must equal 0.0 or 1.0. 

One-tunnel temperature measurement 

T2 = 1.0, T3 = T4 = T5 = 0.0 (Eq.A-1) 

Two-tunnel temperature measurements 

T2 = T3 = 1.0, T4 = T5, = 0.0 (Eq. A-2) 

Note that the numbers 2 through 5 correspond to resistance thermometer 
numbers normally used. 

3. A reference model length, REFL, must be given in units of feet to compute 
model Reynolds number. 

B. Atmosph eric Pressure 

Atmospheric pressure calculation may be handled in the standard program 
for quantities. Its inclusion (if required) and method of obtaining (dialed-in 
optional digital channel or measured by gage in analog channel) is left 
optional to the project engineer. However, measuring atmospheric pressure 
with a gage is recommended rather than entering this pressure reading into 
an analog channel since it is possible for significant variations to occur 
during the course of a tunnel run. 
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C. Mach Number 

1. MCODE indicates which measurements are to be used for Mach number 
calculation (see nomenclature on page A-l). The default value of MCODE 
is 2. Multiple options are provided to allow for the possibility of instrument 
failure during a test. If the digital MCODE input is 1 to 5, then digital 
value overrides the C-card value. If the digital value is zero, then the "C H 
value overrides. The reference pressures may also change. 


If MCODE = 1 

PO/PTO = FXPTANKG/PTH) (Eq. A-3) 

If MCODE = 2 

PO/PTO = F(PTANKH/PTH) (Eq. A-4) 

If MCODE = 3 

PO/PTO = F(PTANKG/PTG) (Eq. A-5) 

If MCODE = 4 

PO/PTO = F(PTANKH/PTG) (Eq. A-6) 

If MCODE = 5 

PO/PTO = F(PTKSON/PTSON) (Eq. A-7) 


whore F(«) is the fitted polynomial from the 1990 calibration 
D. Tunnel Static Pressure 

PO calculation automatically depends on MCODE. No input is required from 
the project engineer. The normal procedure (internal constant MCODE = 2) 
uses PTH for computation. 
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If MCODE £ 2 


PO * (PO/PTO)PTH 


(Eq. A-9) 


If MCODE = 3 or 4 
PO = (PO/PTO)PTG 


(Eq.A-10) 


If MCODE = 5 

PO = (PO/PTO)PTSON 


(Eq.A-11) 


E. Tunnel Tot al Pressure 

Pit) to calculation automatically depends on MCODE. No input is required 
from the project engineer. le normal procedure (MCODE = 2) uses PTH. 


If MCODE $ 2 

PTO = P™ (Eq. A- 12) 

If MCODE = 3 or 4 

PTO = PTG (Eq. A- 13) 

If MCODE = 5 

PTO = PTSON (Eq. A- 14) 


F. Tunnel Dynamic Pffftfmrff 

Tunnel dynamic T . essure is computed as follows: 

If MACH < .1 

QO = PO (Eq. A- 15) 
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If MACH 2; .1 


QO = 0.7 * PO * MACH I. 2 (Eq. A-16) 

G. Dew Point 

Dew point calculation may be handled in the standard program for quantities. 

Its inclusion, channel location, and name are left optional to the project 

engineer. 

H. Tunnel Total Temperature 

1. Provision is made for four individual tunnel total temperature 
measurements. They may be either thermocouples or resistance 
thermometers; however, the appropriate equation must be specified for the 
standard program for quantities. Note that resistance thermometer one 
(1) (strut head) should not be used. If resistance thermometers are used, 
their calibrations are included internal to the program. 

2. The constants required from the project engineer are T2, T3, T4, and T5 
(0.0 or 1.0). 

TTQ _ (RT2 * T2) + (RT3 » T3) + (RT4 * T4) + (RT5 » T5) 

T2 + T3 + T4 + T5 (Eq.A-17) 

I. Reynolds Number 

1. The constant required from the project engineer is REFL. 
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181193 * 10 * PTO * MACh( TTO + 658.27 + 39.72 MACH 2 


RN/FT = 


(TTO+459.67) Z ^L+0.2 MACH 


] 


(EqA-18) 


RN * RN/FT * REFL 


(Eq. A- 19) 


The derivation of the formulas in Appendix A can be found in Ames Research 
Staff, Equations, Tables, and Charts for compressible flow, NACA Report 1135 
(1953). 
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APPENDIX B 

Jet Exhaust Measurements 


Nomenclature 

Required Constants 25 

Test for Exhaust Model 27 

Compute Common Constants 27 

Individual Engine Measurements 27 

Total Exhaust System Properties 33 

Secondary Flow Measurements 39 

Tertiary Flow Measurements 41 


15 




MODULE B JET EXHAUST MEASUREMENTS 


A 




£ 



SYMBOL 

AENG(I) 


AREF 

AT(I) 

AVRI(L) 


C* 

CDSI(L) 

CFI(M) 

CFICHR(M) 

FI(M) 

FICHR(M) 


FIENG(I) 


NOMENCLATURE 

Flow area to be used for determining each engine mass-flow rate 
from plenum chamber measurements, where I = engine 
number. This area is generally based on the area of the plenum 
orifice nozzles (AENG(I) = (orifice areaV2 for twin engines), 
sq. in. 

Model reference area used for coefficients, sq. in. 

Throat area of each engine, where I = engine number, sq. in. 
Area of throat of in-line (not MCV) venturi, where L = venturi 
number, sq. in. 

Critical area, sq. in. 

Discharge coefficient, where L = venturi number. 

Ideal thrust coefficient based on mass-flow rate, where M = air 
system. 

Ideal thrust coefficient based on mass-flow rate obtained from 
plenum chamber measurements, where M = air system. 

Ideal thrust of total primary exhaust system based on measured 
mass-flow rate, lbs., where M = air system. 

Ideal thrust of total primary exhaust system based on mass-flow 
rate obtained from plenum chamber measurements, lbs., where 
M = air system. 

Ideal thrust of individual engines (where I = engine number (up 
to 4)) based on mass-flow rate obtained from individual plenum 
chamber measurements, lbs. 


r 
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SYMBOL 

FM(M) 

FMS 

g 

GAMJ 

IAIR(I) 

ICH(I) 


INTFMl(M) 

INTFMS 

KAE(I) 

KBL 

KCH(I) 

KI1 

KI2 

KI3 

KJ1 

KT2 

KJ3 

KJ4 


NOMENCLATURE 

Primary exhaust flow air flowmeter frequency, hertz, where 
M * air system. 

Secondary flow air flowmeter frequency, hertz. 

Acceleration due to gravity, 32.174 feet per second 
Ratio of specific heats for primary exhaust flow. 

Air system for each engine, where I = engine number. Must be 1 
or 2, default = 0. 

Intercept to be used for determining each engine mass-flow rate 
from plenum chamber measurements, where I = engine 
number. 

Flowmeter number for primary flow air flowmeter, where 
M = air system. 

Flowmeter n umber for secondary flow air flowmeter. 

Constant used in chamber mass-flow calculation, used if second 
order curve fit is required, where I = engine number. 

If set to 1, tertiary flow computation is done. If set to 0, tertiary 
flow computation is omitted. 

Slope to be used for determining each engine mass-flow rate from 
plenum chamber measurements, where I = engine number. 
Internally computed constant. 

Internally computed constant. 

Internally computed constant. 

Internally computed constant (function of GAMJ). 

Internally computed constant (function of GAMJ). 

Internally computed constant (function of GAMJ). 

Internally computed constant (function of GAMJ). 
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NOMENCLATURE 


SYMBOL 

KJ5 


VLIJl 


KPAV(I) 

KPBUJ) 

KPCH(I) 

KPS 

KPS(J) 

KPT(I,J) 

KPTBI/J) 

KPTS(J) 

KR(I,J) 


Internally computed constant (function of GAMJ). 

If no correction is to be made to the pressure probe, then its value is 
set to LO. If the probe is faulty or does not exist, then its value is set 
feoQ. 

Constants used to determine average primary jet total pressure 
ratio from all engines, where I = engine number. These 
constants must equal 0.0 or 1.0. See note. 

Constants used to determine average static pressure in tertiary 
duct, where J = probe number. Must equal 0.0 or 1.0. See note. 
Break pressure for calculation of WPENG(I) for second order 
equations, lbs/sq. in. 

Secondary flowmeter constant (internally computed). 

Constants used to determine average static pressure in 
secondary air duct, where J = probe number. Must equal 0.0 or 
1.0. See note. 

Constants used in computing jet total pressure, where I = engine 
number and J = probe number. These constants must equal 0.0 
or 1.0. See note. 

Constants used to determine average total pressure in tertiary 
duct, where J = probe number. Must equal 0.0 or 1.0. See note. 
Constants used to determine average total pressure in secondary 
air duct, where J = probe number. Must equal 0.0 or 1.0. 

See note. 

Rake constant for each probe in each engine, where I * engine 
number and J * probe number. Must be equal to 0.0 or 1.0. See 
note. 
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SYMBOL 

KSEC 

KSW(M) 


KTAV(I) 

KTT(I,J) 

KV 

KVA(J.M) 

KVARI(L) 


MBLDOT 

MCV(M) 


NOMENCLATURE 

If set to 1, secondary flow computation is done. If set to 0, 
secondary flow computation is omitted. 

Switch for chamber, venturi or flowmeter, where M = air system. 
= -1, Multiple Critical Venturi mass-flow calculation. 

= 0, Flowmeter mass-flow calculation. 

= 1, Chamber mass-flow calculation. 

= 2, In-line (not MCV) venturi mass-flow calculation. 

Constants used to determine average primary jet total 
temperature from all engines, where I = engine number. These 
constants must equal 0.0 or 1.0. See note. 

Constants used in determining primary jet total temperature, 
where I = engine number and J = probe number. These constants 
must equal 0.0 or 1.0. See note. 

Venturi constant, used to account for different venturi 
calibrations. It includes venturi throat area and discharge 
coefficient. 

Constants used to determine average static pressure of multiple 
critical venturi, where J = probe number and M = air system. 
Constants used in the computation of in-line (not MCV) venturi 
weight flow rate, where L = 1 to 4 represents values of Pt/P at 
A/A* of venturi to convert measured static pressure at throat to a 
total pressure and L = 5 to 8 represents averaging factors (must 
be 0.0 or 1.0). 

Ratio of tertiary weight flow to tertiary ideal weight flow 
Venturi meter number, where M = air system. 
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symbol 

MDOT(M) 

MDOTCH 

MSDOT 

NPTE(I) 

NTTE(I) 

NUMENG 

PBL(J) 

PBLAVE 

PCH(I) 

PCHOKE 

PFM(M) 

PFMS 

PS(J) 

PSEC 

PTBUJ) 

PTBLAV 


NOMENCLATURE 

Primary mass-flow rate as measured by flowmeter, slugs/sec., 
where M = air system. 

Primary mass-flow rate as computed from plenum chamber 
measurements, slugs/sec. 

Secondary flow mass-flow rate, slugs/sec. 

Number of total pressure probes in each engine, where 
I = engine number. (Internally computed). 

Number of total temperature probes in each engine, where 
I = engine number. (Internally computed). 

Number of engines in model (maximum of 4). NUMENG = 0 for 
aerodynamics model (no other constants required). 

Static pressure measurements in the tertiary duct (up to 4), 
where J = probe number, lbs/sq. in. 

Average static pressure in the tertiary duct, lbs/sq. in. 

Individual engine-plenum-chamber total pressure, I = e ng ine 
number, lbs/sq. in. 

Primary jet-total-pressure ratio for choked flow. 

Pressure measured at primary flow flowmeter, lbs/sq. in., where 
M = air system. 

Pressure measured at secondary flow flowmeter, lbs/sq. in. 

Static pressure measurements in the secondary flow duct (up to 
4), where J = probe number, lbs/sq. in. 

Average static pressure in the secondary flow duct, lbs/sq. in. 
Total pressure measurements in the tertiary duct (up to 4), 
where J = probe number, lbs/sq. in. 

Average total pressure in the teritary duct, lbs/sq. in. 
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SYMBOL 

PTB/PTJ 

PTB/PTO 

PTENG(I) 

PTENGGVPO 

PTENGO(I) 

PTJ(I.J) 

PTJ/PO(M) 

PTS(J) 

PTS/PTJ 

PTS/PTO 

PTSEC 

PTV 

PV 

PV1 

PV2 

PV/PTV 


NOMENCLATURE 

Ratio of tertiary total pressure to primary jet total pressure. 

Ratio of tertiary toti il pressure to free-stream total pressure. 
Average primary jet total pressure in each engine, where 
I = engine number, lbs/sq. in. 

Ratio of average primary jet total pressure in each engine to 
tunnel static pressure, where I = engine number. 

Ratio of average primary jet total pressure in each engine to 
tunnel static pressure, where I = engine number. 

Individual primary jet total pressure measurements, where 
I = engine number and J = probe number, lbs/sq. in. 

Average primary jet total pressure ratio (all engines), where 
M = air system. 

Total pressure measurements in the secondary flow duct, where 
J = probe number, lbs/sq. in. 

Ratio of secondary flow total pressure to primary jet total 
pressure. 

Ratio of secondary flow total pressure to free-stream total 
pressure. 

Average total pressure in the secondary flow duct, lbs/sq. in. 
Tertiary venturi total pressure, lbs/sq. in. 

Tertiary venturi static pressure, lbs/sq. in. 

Averaged multiple critical venturi static pressure upstream of 
venturi throat, lbs/sq. in. 

Averaged multiple critical venturi static pressure downstream of 
venturi throat, lbs/sq. in. 

Ratio of tertiary venturi static pressure to tertiary total pressure. 
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SYMBOL 

PVEN(I,M) 


PVRI(L,M) 

RDUCT 

RJ 

RNMCV(M) 
R S 
RV 

TCH(I) 

TFM(M) 

TFMS 

THETBL 

THETSE 

TTBL 

TTENG(I) 

TTJ(I,J) 

TTJAVCKM) 

TTSEC 

TTV 

TV(M) 


NOMENCLATURE 

Multiple critical static pressure, where 1=1 and 3 are upstream 
and 1 = 2 and 4 are downstream of venturi throat, lbs/sq. in., 
where M = air system. 

In-line (not MCV) venturi static pressure, lbs/sq. in., where 
L = venturi number and M = air system. 

Venturi throat Reynolds number 

Gas constant for primary flow, ft/degree Rankine. 

Venturi Reynolds number, where M = air system 
Gas constant for secondary flow, ft/degree Rankine. 

Gas constant for tertiary flow, ft/degree Rankine. 

Individual engine-plenum chamber total temperature, 

I = engine number, *F. 

Temperature at primary flowmeter, ‘F, where M = air system. 
Temperature at secondary flowmeter, *F. 

Tertiary flow corrected mass-flow ratio. 

Secondary flow corrected mass-flow ratio. 

Total temperature of tertiary flow, *F. 

Average primary jet total temperature in each engine where 
I = engine number, *F. 

Individual primary jet total temperature measurements where 
I = engine number and J = probe number, *F. 

Average primary jet total temperature (all engines), *F, where 
M = air system. 

Secondary flow total temperature, *F. 

Temperature at the tertiary venturi, *F. 

Multiple critical venturi temperature, *F, where M = air system. 
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SYMBOL 

TVRI(L.M) 

VIS 

VRATIO(M) 

WI(M) 

WIBL 

WIENG(I) 

WMCV(M) 

WMCV/WI(M) 

WP(M) 

WPBL 

WPCHR(M) 

WPCHR/WKM) 

WPENG(I) 

WPSEC 


NOMENCLATURE 

Temperature at the in-line (not MCV) venturi, where L = venturi 
number, *F, and M * air system number. 

Free-stream viscosity, lb. secVsq. ft. 

Ratio of multiple critical venturi static pressures (should be less 
than 0.93), where M * air system number. 

Ideal weight flow of primary flow, lbs/sec., where M s air system 
number. 

Ideal weight flow of tertiary flow, lbs/sec. 

Ideal weight flow of each individual engine primary flow, where 
I = engine number, lbs/sec. 

Multiple critical venturi weight flow rate, lbs/sec., where M = air 
system. 

Ratio of multiple critical venturi weight flow rate to ideal weight 
flow rate, where M = air system. 

Measured weight flow of air primary flow flowmeter or venturi, 
lbs/sec., where M = air system. 

Tertiary weight flow rate obtained from venturi, lbs/sec. 

Total primary flow weight flow rate obtained from plenum 
chamber measurements, lbs/sec, where M * air system. 
Discharge coefficient of total primary flow system as obtained 
from plenum chamber measurements for entire system, where 
M = air system number. 

Primary flow weight flow rate of each engine obtained from 
plenum-chamber measurements, where I = engine number, 
lbs/sec. 

Secondary flow weight flow rate, lbs/sec. 
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SYMBOL 

WP/WI(M) 

WFEAVIE(I) 

WPVRI(M) 

WPVRI/WI(M) 

WVRI(L.M) 

Z 

ZS 


nomenclature 

Primary flow discharge coefficient using flowmeter or venturi 
weight flow rate for entire system, where M * air system. 
Discharge coefficient of each individual engine as obtained from 
plenum-chamber measurements, where I = engine number. 
Sum of in-line (not MCV) venturi weight flow rates, lbs/sec., 
where M = air system. 

Ratio of summation of in-line (not MCV) venturi weight flow rate 
to ideal weight flow rate, where M = air system. 

In-line (not MCV) venturi weight flow rate, lbs/sec., where 
L = venturi number, and M = air system. 

Primary flowmeter constant. (Internally computed). 

Secondary flowmeter constant. (Internally computed). 
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Jet Exhaust Measurements 


A. Required Constanta 

1. All constants are initialized to a value of zero. The project engineer needs 
to supply only those constants which are required for the quantities to be 
computed. In addition, by logical use of combinations of these constants, 
several options are available to the project engineer. One of these options is 
discussed later. 

2. NUMENG - number of engines in model. NUMENG » 0 for aerodynamics 
model (no other constants are required). 

3. KR(I,J) -Rake constant for each probe in each engine, where I = engine 
number and J = probe number. 

If no correction is to be made to the total pressure probe, then its value is set 

equal to 1.0. If the probe is faulty or does not exist, then its value is set equal to 

0 . 0 . 


• Example: Two engines; five probes in the first, and three probes in the 

second. 

Engine 1 is corrected to integrated rake values, engine 2 probes 
are uncorrected. 

KRG,1)= 1.051 
KR(1,2) = .986 
KR(1,3) = .972 
KR(1,4) = .987 
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1. The constant required from the project engineer is NUMENG (0 to 4). 
IF NUMENG = 0, skip module B. 


1. The constants required from the project engineer are GAMJ and RJ. 


KJ2 = 


GAMJ *64,348 
(GAMJ-l)RJ 


(Eq.B-1) 


„ TO _ I 2 *(GAMJ)*(RJ) 

V (GAMJ- !)• 32.1 74 


(Eq. B-2) 


KJ4 = 


/ GAMJ- 1 
1 GAMJ 


(Eq.B-3) 


KJ5 = 


GAMJ 


(Eq.B-4) 


PCHOKE 


■K 


GAMJ 

GAMJ - iY|gamj-i 


(Eq. B-5) 


1. This permits computation for four separate engines with the following 
instrumentation in each engine: 


JKHNIIONAUT HI 


PRECEDING PAGE BUNK NOT FILMED 


a. jet total pressures 

b. jet total temperatures 

c. chamber pressure 

d. chamber temperature 

2. Jet total pressure 

a. Jet total pressure will always be called PTJ(I,J), where I = engine 
number and J = probe number. 

b. An example of representing the third measurement (probe 3) of jet tota? 
pressure in engine 2 is named PTJ(2.3). 

c . The constants required from the project engineer are KR(I,J) and 
KPT(I,J). 


PTENG(I) = 


N1 T (1) PTJ(I ,J)*KR(I,J) 

J= WrE?I5 

I KPT(I,J) 
j^i 


(Eq.B-6) 


PTENGO(I) = 


PTENG(I) 

PO 


(Eq.B-7) 


3. Jet total temperature 

a. Jet total temperature measurements are always called TTJ(I,J), where 
I = engine number and J = probe number. 

b. An example of the first measurement (probe 1) of jet total temperature in 
engine 3 is named TTJ(3.1). 
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c. The constants required from the project engineer are KTT(I,J). 


TTENG(I) = 


NTTEO) 

I TTJ(I, J) * KTT(I, J) 
j=i 



I KTT(I,J) 

J=1 


(Eq. B-8) 


4. Chamber weight flow for each engine. 

a. The constants required from the project engineer are KCH(I), ICH(I), 
KAE(I), AT(I), AENGU) and KPCH(I). 

KJ 1 = 0.5316 + (PTENG(I) + 16.9) ((1.581 - 0.00834(TTENG(I) - 60))/105) 1 (Eq. B-9) 

If PCH(I) < KPCH(I) 
then 


AENG(I) * PCH(I) * KJ 1 * 


WPENG(I) = 


ICH(I) + KCH(I) * PCH(I) + KAE(I) * PCH(I) 
VTCH(I) + 459.67 


If PCH(I) > KPCH(I) then 


AENG(I) * PCH(I) * KJ 1 ‘ 


WPENG(I) = 


ICH(I + 4) + KCH(I + 4) * PCH(I) + KAE(I + 4) * PCH(I) 


VTCH(I) + 459.67 


(Eq.B-10) 


1 Reimer, Robert M.: Computation of the Critical Flow Function, Pressure Ratio, and Temperature Ratio for Real 
Air AS ME Paper 62-WA 177 Journal of Basic Engineering, Trans. ASME 
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5. Ideal weight flow for each engine. 

a. The nozzle choke total pressure ratio is calculated internally is 
called PCHOKE. 

b. The constant required from the project engineer i 3 AT(I). 

If PTENGOd) is greater than FCHOKE, use equation B-ll. 


WIENG(I) = 


[KJ1] « [ PTENG(I)] * [ AT(1)1 
VTTCNG(I)+ 459.67 


(Eq. B-ll) 


If PTENGOd) is less than or equal to PCHOKE, use equation B-12. 


KI1= 


KJ2 

(TTENG(I) + 459.67) 


1 


PTENGOd) 



(Eq. B-12) 


If KI1 is less than 0, KI1 = .0001 
then 


WIENG(I) = Vial ♦ AT(I) * PTENG(I) 


PTENGOd) 



(Eq. B-13) 


Note to the project engineer: If the engine is shrouded, then a local static 

pressure in the nozzle shroud should be used rather than PO. The engineer must 
supply a new equation for KI1 and WIENGd). 
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6. Discharge coefficient for each engine based on chamber weight flow. 


WPE/WIE(I) = 


WPENG(I) 

WIENG(I) 


(Eq.B-14) 


If WIENGG) = 0, WPE/WIEG) = 0 


7. Ideal thrust for each engine based on chamber weight flow. 


KI2 = [TTENG(I) + 459.67] * 



( 1 
^PTENGO(I) 



(Eq. B-15) 


If KI2 is less than 0, KI2 = .0001 

FIENG(I) = [KJ3] * [WPENG(I)] * [VKI2] 

(Eq. B-16) 

8. In-line venturi: weight flow for each Venturi, in each air system. The 
equations given below are for critical flow venturi and are intended to be 
very general. 

Air system number = IAIR(I) = (M) 

A(I) * {[VKRI(4,I) * (TVRI(L,M) + 459.67) + VKRI(3,I)] ♦ (TVRI(L,M) + 459.67) + 
VKRI(2,I)} * (TVRI(L,M) + 459.67) + VKRI(l.I) 

(Eq. B-17) 

A(I) where I * 1 to 4 are constants which go into the compressibility term, C*. 
As seen, a 3rd order equation capability exists. Values of VKRI(1,I) to 
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VKRI(4,I) can be input using 'C' cards to allow use of most any critical 
venturi. 


C* = [(A(4) * PVRI(L,M) + A(3)) * PVRI(L.M) + A(2)] * PVRI(L,M) + A(l) 

(Eq.B-18) 


TS = (TVRI(L,M) + 459.67)/1.2 


VIS = 6.086248 * lO* 8 * (TS) 1 5/(TS + 198.6) 


(Eq.B-19) 


(Eq.B-20) 

Individual venturi mass flow is then computed using 

WVRI(L M) = 1V1^(^ M)*K7AR1(L) * AVR)jL)*^ «C**CDSI 

Vg * RJ * (TVRI(L.M) + 459.67) 


(Eq. B-21) 


M OTE - CDSI represents the discharge coefficient of individual venturi. It is 
obtained using an iterative scheme based on venturi throat Reynolds 
number. A table of CD versus RDUCT is required for each venturi. 
RDUCT is computed using 


RDUCT = WVRI(L)/(AVRI(L) * VIS) 

(Eq. B-22) 

Because of the complexity of this computation, an example is included. 
The following information is contained within the data reduction program 
when using the twin critical venturis which measure total mass flow in 
the groundstand (B1234). 
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VKRI(4,1)= 0.0 
VKRI(3,1)= -1.43545E-8 
VKRI(2,1)= 1.36243E-5 
VKR1(1,1)= 0.68166 
VKRI(4 > 2)= 0.0 
VKRI(3,2)= 4.49456E-10 
VKRI(2,2)= -6.06496E-7 
VKRI(1,2)= 2.14835E-4 

KVARI(l) = 1.0040 
KVARI(2) = 1.0039 


VKRI(4,3) = 0.0 
VKRI(3,3) = 1.64438E-13 
VKRI(2,3) = -1.90568E-10 
VKRI(1,3) = 5.4424E-8 
VKRI(4,4) = 0.0 
VKRI(3,4) = 0.0 
VKRI(2,4) = 0.0 
VKRI(1,4) = 0.0 


AVRI(l) = .272009 KVARI(5) = 1.0 
AVRI(2) = .264481 KVARI(6) = 1.0 


Only the KVARI and AVRI constants are required to be input by an 
engineer. Both venturis use the same CD versus RDUCT relationship, 
which is not a table lookup but simply a second order equation. Of course 
table lookup could be used in lieu of the equation. 

The CDSI equation for twin critical venturis in groundstand: 

CDSI = 0.993607 + 3.5062E-4*(EDUCT) - 1.1269E-5 * (RDUCT 2 ) 

where RDUCT = WVRI(L,My(AVRI(L) ♦ VIS * 1.0E06) 

E- Total Exhaust Syste m Properties 
1. Average total pressure ratio. 

a. The constants required from the project engineer are KPAV(I) and 
IAIR(I). 


(M) = IAIR(I) 


numenq „ pTENGOd)] 

PTJ / PCXM) = >*> NUMENff^— - 

1=1 


(Gq. B-23) 


2. Average total temperature. 

a. The constants required from the project engineer are KTAV(I) and 
IAIR(I). 


(M) = LAIRG) 


TTJAVG(M) = 


N T NG [KTAV(I) * TTENG(I)] 

1=1 NUMENC 

I KTAV(I) 

i=i 


(Eq.B-24) 


3. Total weight or mass flow. 

a. Each air system weight flow is in units of lb/sec. 

b. Each air system mass flow is in units of slugs/sec. 

c. The constants required from the project engineer are: 

(1) INTFMl(M) and MCV(M), where M is air system number. 

(2) KSW selects mass flow computation 

= 1; chamber flow 
= 0; flowmeter 
ss-1; multiple critical venturi 
s 2; in-line venturis 

If KSW = l (chamber mass flow calculation) 
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(M) = IAIRG) 


WPCffiUM) = NU ! EN $PENG(I) 

1=1 

(Eq.B-25) 


MDOTch - WPCHR CM) 
32.174 


(Eq.B-26) 


If KSW = 0 (air model with flowmeter) 

Z and KP are determined from standardized flowmeter tables and from 
INTFM 1 (M) constant. 


WP(M) = 


(FM1(M)) * (PFM(M) * (144. ) 
(RJ) * (Z) * (KP) * (TFM(M) + 459.67) 


(Eq. B-27) 


MDOT(M) = 


WP(M) 

32.174 


(Eq. B-28) 


If KSW = -1 (venturi mass flow calculation), the multiple critical venturi code, 
MCV, is decoded to derive those venturi present 


PV1 = 


KVA1* PVENg M) + KVA3 ♦ PVEN(3,M) 
KVA1+KVA3 


PV2 = 


KVA2 ♦ PVEN(2,M) + KVA4 ♦ PVEN(4,M) 
KVA2+KVA4 


(Eq. B-29) 
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VRATIO(M) = — — 
PV1 


(Eq.B-30) 


A(D = (CVK(4,I) ♦ TV(M) + VK(3,I)) * TV(M) + VK(2,I)) * TV(M) + VK(1,I) 

(Eq. B-31) 

C* = ((A(4) * PV1 + A(3)) * PV1 + A(2)) * PV1 + A(l) 

TS = (TV(M) + 459.67Y1.2 
VIS = 6.086248 * 10* * (TS)l-5/(TS + 198.6) 


(Eq. B-32) 
(Eq.B-33) 
(Eq. B-34) 


WMCV(M) = I PV 1 * AREAV(I) * (O) ♦ f 32.174 \ y2 , CDm 

I l(TV(M) + 459.67)RjJ CD(I) 


(Eq. B-35) 


ARMCV = I AREAV(I) 

I 

(Eq. B-36) 

The above summations are over the venturi present. CD(I) is computed by 
linear interpolation from a table of CD vs RNMCV 
where 

RNMCV(M) = WMCV(M)/(ARMCV*VIS) 

(Eq. B-37) 

An iterative scheme is used until successive computations of WMCV differ 
by a desired accuracy. 

4. If KSW = 2 (in-line venturis) 


WVFKM) = £ WVRI(L.M) * KVARI(L + 4) 


(Eq.B-38) 


5. Ideal weight flow (total). 

a. Ideal weight flow of each air system is computed 


WI(M) = 


NUMENG 

I 

1=1 


WIENG(I) 


(Eq.B-39) 


6. Discharge coefficient for each air system. 

a. The discharge coefficient using weight flow from a flowmeter or a 
venturi and the discharge coefficient using weight flow from chamber 
measurements are computed. 


If KSW = 2 

KSW = 1 
KSW = 0 
KSW = -1 


WP(M) = WPVRI(M) 
WP(M) = WPCHR(M) 
WP(M = WP(M) 
WP(M) = WMCV(M) 


MDOT(M) = 

32.174 


WP/WI(M) = - WP(M) 
WI(M) 


(Eq. B-40) 


(Eq.B-41) 


WPCHR / WI(M) = WPCHR(M) 

WI(M) 


(Eq.B-42) 
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WMCV / WI(M) = 


WPVRI / WI(M) = 


WPVRI(M) 

WI(M) 


(Eq. B-43) 


(Eq. B-44) 


If WI(M) = 0; WP/WI(M) = WPCHR(M)/WI(M) = WMCV(MVWKM) * 
WPVRI(MyWI(M) = 0 


7. Ideal thrust for each air system. 

a. The ideal thrust, FICHR(M), and ideal thrust coefficient CFICHR(M) 
are obtained from chamber weight flow. 

b. The ideal thrust, FI(M), and ideal thrust coefficient CFI(M) are 
obtained from flowmeter or venturi measured weight flow. 

c. Note that MACH, PO and QO are from Module A and M = air system. 

d. The constant required from the project engineer is AREF. 


FICHR(M) = 


NUMENG 

I 


FIENG(I) 


l=i 


(Eq. B-45) 


If MACH < .1, 


CFICHR(M) = 


FICHR(M) 

(PO)*(AREF) 


If MACH ;> .1, 


(Eq.B-46) 
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CFICHR(M) = 


FICHR(M) 

(QO)*(AREF) 


KI3 = (TTJAVG + 459.67) * 


.-f.-L-f 4 ' 

VPTJ/PO(M)J 


(Eq.B-47) 


(Eq.B-48) 


If KI3 < 0; KI3 = .0001 


FI(M) = (KJ3) * (WP(M)) * (VKI3) 


(Eq. B-49) 


If MACH < .1, 


CFI(M) = 


FI(M) 

(PO) * ( AREF) 


(Eq. B-50) 


If MACH * .1, 


CFI(M) = 


FI(M) 

(QO)*(AREF) 


(Eq. B-51) 


F. Secondary Flow Measurementa 
If KSEC = 0, skip equations B-52 through B-58. 
1. Secondary passage total pressure. 
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a. The total pressure measurements PTS(J) in the secondary air passage 
(up to 4) are used to compute the average secondary passage total 
pressure. 

b. The constant required from the project engineer is KPTSKJ). 

I KPTS(J)*PTS(J) 

PTSEC=s^1 — 2 

I KPTSKJ) 

J=i 

(Eq.B-52) 


2. Secondary passage static pressure. 

a. Static pressure measurements PS(J) in the secondary air passage (up to 
4) are used to compute the average static pressure in the secondary air 
passage. 

b. The constant required from the project engineer is KPSKJ) 

I KPSKJ) *PS(J) 

PSEC = 

I KPSKJ) 

J=i 

(Eq.B-53) 


3. Secondary duct total temperature 

a. The total temperature TTSEC in the secondary duct is handled in the 
standard program for quantities. 

4. Secondary mass flow 

a. The constants required from the project engineer are RS, KPS, ZS, 
INTFMS. KPS and ZS are determined internally from INTFMS 
constant. 
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* 


WPSEC = 


(FMS)*(PFMS)* (144,0) 

(RS) * (ZS) * (KPS) * (TFMS + 459.67 ) ,lhs/6ec 


(Eq.B-54) 


MSDOT= 


WPSEC 
32.174 ’ 


slugs / sec 


(Eq.B-55) 


5. Pumping characteristics 


PTS / PT J = 


PTSEC 

(PTJ / PO(M)) * (PO) 


(Eq. B-56) 

PTS/PTO = ^ TSEC 
PTO 

(Eq. B-57) 


6. Corrected mass flow ratio 


THETSE = MSDQT f (TTSEC-f 459.67) *RS 
MDOT(M) V (TTJAVG(M) + 459.67) * RJ 


(Eq. B-58) 


G. Tertiary Flow Measurempnt fl 

If KBL = 0, skip equations B-59 through B-64. 

1. Tertiary duct total pressure. 

a. The total pressure measurements PBL(J) in the tertiary duct (up to 4) 
are used to compute the average tertiary duct total pressure. 

b. The constant required from the project engineer is KPTBUJ). 
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I KPTBL( J) * PTBL( J) 

PTBLAV * — ? 

I KPTBL(J) 

J=i 

(Eq. B-59) 


2. Tertiary duct static pressure. 

a. Static pressure measurements PBL(J) in the tertiary duct (up to 4) are 
used to compute the average static pressure in the tertiary duct. 

b. The constant required from the project engineer is KPBL(J). 

I KPBL( J) * PBL(J) 

PBLAVE = ^ 5 

I KPBL(J) 

J=i 

(Eq.B-60) 


3. Tertiary duct total temperature. 

a. Total temperature in the tertiary duct TTBL is handled in the standard 
program for quantities. 


4. Tertiary mass flow. 

a. Venturi total pressure, PTV, and venturi static pressure, PV, are 
required. 

b. Tertiary weight flow is in units of lbs/sec. 

c. Tertiary mass flow is in units of slugs/sec. 

d. The constants required from the project engineer are RV, KV. 


PV/PTV = 


PV 

PTV 
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048117 JTTV +45967 

WPBL = 0. 13594KPTV - PV) / (PTV / 14.696)] * (PTV / 14.696) / ^ — 


(Eq. B-61) 


0.857143 0.2857 

WIBL = 0.72167 *(PV/ PTV) *\(PV/PTV) - 1* PTV /VTTV+ 459.67 

(Eq. B-62) 


MBLDOT = WPBIVWIBL 


5. Pumping characteristics. 


(Eq.B-63) 


PTB/PTJ = 


PTBLAV 

(PTJ / PO(M)) * (PO) 


(Eq. B-64) 


PTB/PTO = 


PTBLAV 

PTO 


(Eq.B-65) 
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APPENDIX C 
Skin Friction Drag 


Nomenclature 

Required Constants 

Test for Skin Friction Calculation 

Fuselage Skin Friction Drag 

Empennage Skin Friction Drag... 
Total Skin Friction Drag 
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MODULE C SKIN FRICTION DRAG 


f 

[ 


h 

r 

i 


i 


t 


> 


w 

\ 


r 

I 


SYMBOL 

AREF 

AWET(I) 


CDF 

CDFAFT 

CDFF 

CDFHT 

CDFNOZ 

CDFR(I) 

CDFVT 

CDFW 

FL(I) 


NOMENCLATURE 

Model reference area used for coefficients, sq. in. If module B 
is used, this constant is already specified. 

Model wetted areas, sq. in. 

Where AWET(1) = total fuselage wetted area. 

AWET(2) = fuselage wetted area up to metric break 
AWET(3) = fuselage wetted area up to nozzle connect 
station. 

AWET(4) = wing wetted area. 

AWET(5) = vertical tail vetted area. 

AWET(6) = horizontal tail wetted area. 

AWET(7) = optional, for additional body. 

Total skin friction drag coefficient. 

Afterbody plus nozzle skin friction drag coefficient. 

Total fuselageawkin friction drag coefficient. 

Horizontal tails (canards) skin friction drag coefficient. 

Nozzle skin friction drag coefficient. 

Individual skin friction drag coefficients calculations. 

Vertical tails(s) skin friction drag coefficient. 

Wing skin friction drag coefficient. 

Model reference lengths, feet. 

Where FL(1) = fuselage length. 

FL(2) = fuselage length up to metric break. 

FL(3) = fuselage length up to nozzle connect station. 

FL(4) = wing mean aerodynamic chord. 
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SYMBOL 


PORMF(I) 


KFAFT 


KFF 

KFNOZ 


NOMENCLATURE 

FL(5) = vertical tail mean aerodynamic chord. 

FL(6) = horizontal tail mean aerodynamic chord. 

FL(7) s optional. 

Form factors 

Where FORMF(l) * fuselage. 

FORMF(2) = wing. 

FORMF(3) = vertical tail. 

FORMF(4) = horizontal tail. 

FORMF(5) = optional. 

Constant used to include proper terms in total skin friction drag 
term, CDF. Must equal 0.0 or 1.0. If the relevant term is to be 
incorporated to the total skin friction drag, set to 1.0; otherwise, 
set to 0.0. 

See KFAFT. 

See KFAFT. 


47 


APPENDIX C 
Module C 
Skin Friction Drag 


Skin friction drag is computed by the method of Frankl and Voishei 2 for 
compressible, turbulent flow on a flat plate. 

A. Required Constanta 

All constants are initialized to a value of 0.0 except FORMF(I) which is 
initialized to a value of 1.0. 

1. AWET(I) 

2. FORMF(I) 

Form factors may be obtained from LWP - 1120. 

Fuselage: FORMF(I) = 1.0 + l.SCd/l) 1 -® + 7(d/l)3 (Eq. C-l) 

Empennage: FORMF(I) = 1.0 + 1.44(t/c) + 20/c) 2 (Eq. C-2) 

3. The model reference lengths (FL(I)), are given in the nomenclature 
section. 

4. The model reference area (AREF) is used for coefficients, sq. in. If jet 
exhaust measurements are used, this constant is already specified. 


2 Frankl, F.; and Voishei, V.: Friction in the Turbulent Boundary Layer of a Compressible Gas c? 
High Speeds. TM NACA No. 1032, 1942. 
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5. The constants (KFF, KFAFT, KFNOZ) used to include proper terms in total 
skin friction drag term, CDF*, must equal 0 or 1. 


B. 


Teat for Skin Fri 




Calculation 


If AWET(1) = 0, skip the calculations for the skin friction drag in this module. 


C. Fuselage Skin Friction Drag 

1. The constants required from the project engineer are AWET(1), AWET(2), 
AWET(3), FL(1), FL(2), FL(3), AREF, and FORMF(l). 

J = 3 

If AWET(2) = 0 and AWET(3) = 0, J = 1 
If AWET(2) * 0 and AWET(3) = 0, J = 2 
Calculate CDFR(I) for I = 1, J 


CDFR(I) = 


.472 * AWET(I) * FORMF(l) 


1 + .2MACH 


258 

*{log 10 [(RN/FT)*FL(I)]} • 


• AREF 


(Eq. C-3) 


If MACH < .1, CDFR(I) = 0.0 

CDFF * CDFRd) (Eq. C-4) 

If AWET(2) * 0, 

CDFAFT = CDFR(1) - CDFR(2) (Eq. C-5) 

If AWET(3) * 0, 

CDFNOZ * CDFRd) - CDFR(3) (Eq. C-6) 
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w 


D - Emnfinna ? a Skin FricHft» p rfir 

1. The constants required from the project engineer are AWET(4), AWET(5), 
AWETX6), FL(4), FL(5), FL(6), AREF, FORMF(2), FORMF(3), FORMF(4), 
KFF, KFAFT, and KFNOZ. 



Calculate CDFR(I) for I = 4, 7 


J = I - 2 


CDFR(I) = 


.472* AWET(I) * FORMF(J) 


( 2V i67 

|^1+.2MACH I * {log 10 [(RN / FT) * FLd)]} 2,68 ♦ AREF 

(Eq. C-7) 


IF MACH < .1, CDFR(I) = 0 
CDFW = CDFR(4) 
CDFVT = CDFR(5) 
CDFHT = CDFR(6) 


Total Skin Friction 

1. Skin friction drag of the entire model is computed. 

CDF = (KFF * CDFF) + (KFAFT * CDFAFT) + (KFNOZ * CDFNOZ) + 

CDFW + CDFVT + CDFHT + CDFR(7) 

(Eq. C-8) 
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APPENDIX D 

Balance Loads and Model Attitudes Calculations 


Nomenclature ® 

Required Constants 72 

Test for Balance Loads and Model Attitudes 77 

Balance Component Naming System.... 78 

Uncorrected Balance Quantities 78 

Tunnel Support Pitch Angle ® 

Balance Quantities Corrected for Interactions, Weight Tares 

and Momentum Tares ® 

Balan ce Quantities Corrected for Method of Attachment 91 

Angle of Attack and Sideslip Angle 92 

Body Axis Components; Rotation and Translation from 

Balance-to -Model Axis 96 

Pressure Corrections to Body Axis Components 97 

Stability Axis Components 99 

Wind Axis Components 

Alternate Reference Axis Components 101 

Base Force and Moment Tare Coefficients 103 

Base Pressure Coefficients 104 

Model (Body) Axis Coefficients 106 

Stability Axis Coefficients **106 

Wind Axis Coefficients 10® 

Alternate Reference Axis Coefficients 107 

Miscellaneous Equations 107 

Calculation of Initial Weight Tares and Attitude Load Constanta 108 
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Calculation ox Initial Weight Tares and Attitude Load Tares 

(Strain Gage Balance) 

Calculation of Attitude Load Constants by Method I 

Calculation of Attitude Load Constants by Method II 118 

Balances Without Six Components 119 

Initial Weight Tare Calculations 

New Values of Initial Weight Tares 120 

Point Calculations 120 

Figures 

D-l Definition of gravity and balance axes showing positive 

directions and rotation angles for gravity-to-b alance 

transformations 122 

D-l(a) Gravity* to- tunnel support axes 122 

D-l(b) Tunnel support to undeflected balance axes 19a 

D-l(c) Illustrations of knuckle angles 124 

D-l(d) Undeflected balance- to-bal ance ax es iss 

D-l(e) Final balance orientation; gravity-to-balance axes 126 

D-l(f) Illustration of primary balance to undeflected secondary 

balance rotations 127 

D-l(g) Definition of initial or wind-off balance attitude 128 

D-2 Model-balance orientation stings 129 

D-2(a) Case 1, normal balance arrangement 129 

D-2(b) Case IA, case 1 held by opposite end 130 

D-2(c) Case 2, balance rolled 90* (clockwise) 131 

D-2(d) Case 2A, same as case 2 held by opposite end 132 

D*2(e) Case 3, balance rolled 180* (inverted) 133 
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D-2(f) Case 3A, case 3 held by opposite end 134 

D-2(g) Case 4, balance rolled 90° (counterclockwise) 135 

D-2(h) Case 4, held by opposite end 136 

D-2(i) Case 5, balance yawed 180“ or (pitched 180*and rolled 

180“) - reversed 137 

D-2(j) Case 5A, case 5 held by opposite end 138 

D-2(k) Case 6, balance yawed 180“ and rolled 90* clockwise 

(reversed and rolled 90' clockwise), or pitched 180“ 

and rolled 90“ counterclockwise 139 

D-2( 1 ) Case 6A, case 6 held by opposite end 140 

D-2(m) Case 7, balance yawed 180’ and rolled 180“ (reversed 

and inverted), or pitched 180“ 141 

D-2(n) Case 7A, case 7 held by opposite end 142 

D-2(o) Case 8, balance yawed 180“ and rolled 90* counter- 

clockwise (reversed and rolled 90* counter-clockwise), 

or pitched 180* and rolled 90* clockwise 143 

D-2(p) Case 8A, case 8 held by opposite end — 144 

D-3 Definition of gravity and wind axes showing positive 

directions and rotation angles for wind-to-gravity 

transformations 145 

D-4 Definition of balance and body axes showing positive 

directions and rotation angles for balance-to-model 

(body) transformations 145 

D-4(a) Rotation 145 

D-4(b) Translation 145 

D-5 Definition of angle of attack and angle of sideslip 148 


MODULE D BALANCE LOADS AND MODEL ATTITUDES 


S ^ M£QL NOMENCI.ATTrR.TT 

The arrays FO through F20 are forces and moments whose units are 
lbs and in. lbs. 


AF(I,J) 

AFO(I) 

AFT(I) 

AFTARE(I) 

ALPHA 

AMOM(I) 

ARB(II,K) 


ARP(II.K) 


Axial force, lbs., where I = balance number and 
J = correction number. 

Initial axial load, lbs., where I = balance number. 

Total axial load, lbs., where I = balance number. 

Axial weight tares, lbs., where I = balance number. 

Model angle of attack, degrees. 

Axial force momentum correction, lbs., where I = balance 
number. 

Areas or moment arms * areas used with PBASE(II) for 
computing base force and base moment tares. Care 
should be used to insure proper tare force signs. Area 
and arm units must be consistent with units of base 
pressures and balance components. Second balance, sq. 
in., where K = component number and II = orifice 
number. 

Areas or moment arms * areas used with PBASE(II) for 
computing base force and base moment tares. Care 
should be used to insure proper tare force signs. Area 
and arm units must be consistent with units of base 
pressures and balance components. Third balance, sq. 
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SYMBOL 



; 

I 

v 

? 


| 

1 

i 

I 

i. 

» 

f 

f 

i 

f 

r 

r 

i 





m 



ARPB(II.K) 


^0 

^3 


A 4 

BETA 

BSPAN(I) 

CA(1) 

CABASE(I) 

CAREF(I) 

CC(1) 


i 

NOMENCLATURE 

in., where K = component number and II = orifice 
number. 

Areas or moment arms * areas used with PBASE(II) for 
computing base force and base moment tares. Care 
should be used to insure proper tare force signs. Area 
and arm units must be consistent with units of base 
pressures and balance components. First balance sq. 
in., where K = component number and II = orifice 
number. 

Initial balance loads, axial force, lbs. (Weight Tares) 

Balance component quantity corrected for high 
interactions coupled with high model restraints, axial 
force, lbs. (Weight Tares) 

Balance component quantities corrected for balance 
•orientation to gravity axis, axial force, lbs. (Weight Tares) 

Angle of sideslip, degrees. 

Roll and yaw moments reference length. Normally wing 
span, inches, where I = balance number. 

Axial force coefficient in the body axis, where I = balance 
number. 

Base axial force coefficient, where I = balance number. 

Axial force coefficient in the reference axis, where 
I = balance number. 

Crosswind coefficient in the wind axis, where I = balance 
number. 
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SY MBOL 

CD(I) 

CDBASE(I) 

CDS(I) 

CHORD(I) 

CUI) 

CLS(D 

CLSQR(D 

CMX(D 

CMXREF(D 

CMXS(D 

CMXW(D 

CMY(D 

CMYREF(D 

CMYS(D 


NOMENCLATURE 

Drag coefficient in the wind axis, where I = balance 
number. 

Base drag coefficient, where I = balance number. 

Drag coefficient in the stability axis, where I = balance 
number. 

Pitching moment reference length. Normally wing mean 
aerodynamic chord, inches, where I = balance number. 
Lift coefficient in the wind axis, where I = balance 
number. 

Lift coefficient in the stability axis, where I = balance 
number. 

Lift coefficient squared, where I = balance number. 
Rolling moment coefficient in the body axis, where 
I = balance number. 

Rolling moment coefficient in the reference axis, where 
I = balance number. 

Rolling moment coefficient in the stability axis, where 
I = balance number. 

Rolling moment coefficient in the wind axis, where 
I = balance number. 

Pitching moment coefficient in the body axis, where 
I = balance number. 

Pitching moment coefficient in the reference axis, where 
I = balance number. 

Pitching moment coefficient in the stability axis, where 
I = balance number. 
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SYMBOL 

CMYW(I) 

CMZ(I) 

CMZREF(I) 

CMZS(I) 

CMZW(I) 

CN(I) 

CNBASE(I) 

CNREF(I) 

CPBASE(II) 

CPMBASEI) 

CRMBASE(I) 

CY(I) 

CYBASE(I) 

CYMBASE(I) 


NOMENCLATURE 

Pitching moment coefficient in the wind axis, where 
I = balance number. 

Yawing moment coefficient in the body axis, where 
I = balance number. 

Rolling moment coefficient in t’le reference axis, where 
I = balance number. 

Yawing moment coefficient in the stability axis, where 
I = balance number. 

Yawing moment coefficient in the wind axis, where 
I = balance number. 

Normal force coefficient in the body axis, where 
I = balance number. 

Base normal force coefficient, where I = balance number. 
Normal force coefficient in the reference axis, where 
I = balance number. 

Base pressure coefficient, where II = orifice number. 
Base pitching moment coefficient, where I = balance 
number. 

Base rolling moment coefficient, where I = balance 
number. 

Side force coefficient in the body axis, where I = balance 
number. 

Base side force coefficient, where I = balance number. 
Base yawing moment coefficient, where I = balance 
number. 


58 


SYMBOL 

CYREF(I) 

CYS(I) 

Cl 

C2 

AA 

A*i 

A/ 2 

Amx 
Am 2 
AN 

Anj 

An2 

DPBASE(II) 

AW(I) 


AY 

FA 

FA(I) 

FA(I,L) 

FABASE(I) 


NOMENCLATURE 

Side force coefficient in the reference axis, where 
I = balance number. 

Side force coefficient in the wind axis, where I = balance 
number. 

Linear balance interactions. 

Nonlinear balance interactions. 

W(AF), axial force weight tares, lbs. 

WY(RM), rolling moment weight tares, in. lb. 

WZ(RM), rolling moment weight tares, in. lb. 

WX(PM), pitching moment weight tares, in. lb. 

WZ(PM), pitching moment weight tares, in. lb. 

W(NF), normal force weight tares, lbs. 

WX(YM), yawing moment weight tares, in. lb. 

WY(YM), yawing moment weight tares, in. lb. 

Differential base pressures, where II = orifice number. 
Half weight of balance, lbs., where I = balance number. 
Used in weight tares program. If AW(I) is zero and 

DELW is non- zero from the balance interaction deck, then 
AW(I) is set equal to DELW. 

W(SF), side force weight tares, lbs. 

Axial force, lbs. 

Final body axis axial force, lbs., where I * balance 
number. 

Balance axial force rotated (L = 1) and translated (L = 2) to 
body axis, lbs., where I = balance number. 

Base axial force, lbs., where I = balance number. 
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SYMBOL 

FAMAX 


FC(I) 


FD(1) 


FDS(I) 


FLO) 


FLS(I) 


FNO) 


FN(L,I) 


FNBASE(I) 

FNMAX 


Maximum absolute value of axial force, lbs. 


FAMOM(I) Axial force due to momentum of flow, lbs., where 
I = balance number. 

FAREF'(I) Axial force rotated to reference axis, lbs., where 
I = balance number. 

FAREF(I) Axial force translated to reference axis, lbs., where 


I = balance n amber. 

Crosswind foice in the wind axis, lbs., where I = balance 
number. 

Drag force in the wind axis, lbs., where I = balance 
number. 

Drag force in the stability axis, lbs., where I = balance 
number. 

Lift force in the wind axis, lbs., whore I = balance 
number. 

Lift force in the stability axis, lbs., where I = balance 
number. 

Normal force, lbs. 

Final body axis normal force, lbs., where I = balance 
number. 

Balance normal force rotated (L = 1) and translated (L = 2) 
to body axis, lbs., where I = balance number. 

Base normal force, lbs., where I = balance number. 
Maximum absolute value of normal force, lbs. 


FNREF’(I) Normal force rotated to reference axis, lbs., where 

1 = balance numlu . 


Si 


SYMBOL 

NOMENOI-ATtmE 

FNREF(I) 

Normal force translated to reference axis, lbs., where 
I = balance number. 

FP 

All product combinations of vector FT. 

FT 

Corrected total loads. 

FT ARE 

Tare loads. 

FUT 

Uncorrected total loads. 

FY 

Side force, lbs. 

FY(I) 

Final body axis side force, lbs., where I = balance number. 

FY(I,L) 

Balance side force rotated (L = 1) and translated (L * 2) to 
body axis, lbs., where I = balance number. 

FYBASE(I) 

Base side force, lbs., where I = balance number. 

FYMAX 

Maximum absolute value of side force, lbs. 

FYREF’(I) 

Side force rotated to reference axis, lbs., where I = balance 
number. 

FYREF(I) 

Side force translated to reference axis, lbs., where 
I = balance number. 

FYS(I) 

Side force in the stability axis, lbs., where I = balance 
number. 

FO 

Initial loads. 

FI 

Uncorrected balance quantities. 

F2 

Balance component quantities corrected for interactions. 

F3 

Balance component quantities corrected for high 
interactions coupled with high model restraints. 

F4 

Balance quantities corrected for balance orientation to 
gravity axis, attitude loads, and weight tares. 

F6 

Balance quantities corrected for method of attachment. 
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SYMBOL 

F6 

F7 

F8 

F9 

F10 

Fll 

F12 

F13 

F14 

F15 

F16 

F17 

F18 

F19 

F20 

HIRXX(I) 


IGRND(I) 

KMOM 


NOMENCLATURE 

Balance components rotated to the model (body) axis. 
Balance components rotated and translated to the model 
(body) axis. 

Differential base pressure forces. 

Base force and moment tares. 

Final body axis components. 

Stability axis components. 

Wind axis components. 

Rotation from body axis to reference axis. 

Alternate reference axis coefficients. 

Reference axis coefficients. 

Base force and moment tare coefficients. 

Base pressure coefficients. 

Model (body) .ixis coefficients. 

Stability axis coefficients. 

Wind axis coefficients. 

Corrections for the effect of having a model with high 
restraints coupled with high interactions, where XX is the 
balance component (AF, SF, NF, RM, PM, YM) and 
I = balance number. 

Grounding of balance, where I = balance number. 

Axial momentum correction term. 

= 0, no correction. 

- 1> applies nonblowing correction only and automatically 
computes FAMOM(I) 

* 2, applies nonblowing and blowing corrections 
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SYMBOL 

KPP 


KSIGN(I) 


Ka.1 

Ka.2 

Ka.3 

L/D(I) 

LS/DSG) 

h 


U 

METHOD 

MX(I) 


NOMENCLATURE 

A units conversion factor, initialized at 1. If PBASE is in 
PSF and PO is in PSI, KPP = 144.0 If PBASE is in PSI and 
PO is in PSF, KPP = 0.00694 If PBASE is differential 
(PBASE-PO), KPP = 0.0 If PBASE is absolute, KPP = 1.0 
(Standard). 

Constant for correcting balance quantities for gro undin g 
by wrong end, where I = balance number. KSIGN = 1 for 
normal balance attachment. KSIGN = -1 for grounding 
balance by wrong end. 

COS(THETA0) * COS(PHI0) 

SIN(THETAO) 

COS(THETAO) * SIN(PHIO) 

Lift over drag ratio in the wind axis, where I = balance 
number. 

Lift over drag ratio in stability axis, where I = balance 
number. 

Initial balance loads, roll moment, in. lb. 

Balance component quantity corrected for high 
interactions coupled with high model restraints, roll 
moment, in. lb. 

Balance component quantities corrected for balance 
orientation to gravity axis, roll moment, in. lb. 

Method to be used in the weight tares program. 

Final body axis rolling moment, in. lb., where I = balance 
number. 
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SYMBOL 

MX(I,L) 

MXREF'(I) 

MXREF(I) 

MXS(I) 


NOMENCLATURE 

Balance rolling moment rotated (L = 1) and translated 
(L = 2) to body axis, in. lb., where I = balance number. 
Rolling moment rotated to reference axis, in. lb., where 
I = balance number. 

Rolling moment translated to reference axis, in. lb., 
where I = balance number. 

Rolling moment in the stability axis, in. lb., where 
I = balance number. 


MXW(I) 

MY(I) 

MY(I,L) 

MYREF'(I) 

MYREF(I) 

MYS(I) 

MYW(I) 

MZ(I) 


Rolling moment in the wind axis, in. lb., where 
I = balance number. 

Final body axis pitching moment, in. lb., where 
I = balance number. 

Balance pitching moment rotated (L = 1) and translated 
(L = 2) to body axis, in. lb., where I = balance number. 
Pitching moment rotated to reference axis, in. lb., where 
I = balance number. 

Pitching moment translated to reference axis, in. lb., 
where I = balance number. 

Pitching moment in the stability axis, in. lb., where 
I = balance number. 

Pitching moment in the wind axis, in. lb., where 
I = balance number. 

Final body axis yawing moment, in. lb., where I = balance 
number. 
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SYMBOL 

MZ(I,L) 

MZREF'(I) 

MZREF(I) 

MZS(I) 

MZW(I) 

m o 

m 3 


m 4 


NF(I,J) 

NFO(I) 

NFT(I) 

NFTARE(I) 

NUBAL 

n 0 

n 3 


NOMENCLATURE 

Balance yawing moment rotated (L = 1) and translated 
(L = 2) to body axis, in. lb., where I = balance number. 
Yawing moment rotated to reference axis, in. lb., where 
I - balance number. 

Yawing moment translated to reference axis, in. lb., 
where I = balance number. 

Yawing moment in the stability axis, in. lb., where 
I = balance number. 

Yawing moment in the wind axis, in. lb., where 
I = balance number. 

Initial balance loads, pitch moment, in. lb. 

Balance component quantity corrected for high 
interactions coupled with high model restraints, pitch 
moment, in. lb. 

Balance component quantities corrected for balance 
orientation to gravity axis, pitch moment, in. lb. 

Normal force, lbs., where I = balance number and 
J = correction number. 

Initial normal load, lbs., where I = balance number. 
Total normal load, lbs., where I = balance number. 
Normal weight tares, lbs., where I = balance number. 

Number of balances in the model, (max 5). 

Initial balance loads, yaw moment, in. lb. 

Balance component quantity corrected for high 
interactions coupled with high model restraints, yaw 
moment, in. lb. 
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SYMBOL 

n 4 


N 0 

N 3 


n 4 


PBASE(II) 

PHI 

PHIB 

PHIB2 

PHIB3 

PHID 

PHID2 

PHID3 

PHIDX(I) 

PHIK 

PHIK2 


NOMENCLATURE 

Balance component quantities corrected for balance 

orientation to gravity axis, yaw moment, in. lb. 

Initial balance loads, normal force, lbs. 

Balance component quantity corrected for high 
interactions coupled with high model restraints, normal 
force, lbs. 

Balance component quantities corrected for balance 
orientation to gravity axis, normal force, lbs. 

Base pressure, lbs/sq. in., where II = orifice number. 
Model Euler roll angle, degrees. 

Euler roll rotation angle between primary balance and 
model, degrees. 

Euler roll rotation angle between secondary balance and 
model, degrees. 

Euler roll rotation angle between tertiary balance and 
model, degrees. 

Roll deflection of primary balance, degrees. 

Roll deflection of secondary balance, degrees. 

Roll deflection of tertiary balance, degrees. 

Deflection roll angle constants, where X is balance 
component (A, S, N, R, P, Y) and I = balance number. 
Euler roll angle to account for knuckle and/or primary 
balance angles in relation to tunnel support, degrees. 
Euler roll angle to account for orientation of undeflected 
secondary balance in relation to primary balance, 
degrees. 
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SYMBOL 

PHIK3 

PHIR 


PHIS 

PHI0,I 

PM 

PM(I,J) 

PMBASE(I) 

PMMAX 

PMO(I) 

PMT(I) 

PMTAREI 

PSI 

PSIB 

PSIB2 

PSIB3 

PSID 

PSID2 


NOMENCLATUR E 

Euler roll angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 
Euler roll rotation angle between model (body) axis and 
reference axis, positive in same direction as PHIB, 
degrees. 

Strut roll angle, degrees. 

Wind off zero attitude of each balance, degrees, where 
I = balance number. 

Pitching moment, in. lb. 

Pitching moment, in. lb., where I = balance number and 
J = correction number. 

Base pitching moment, in. lb., where I = balance number. 
Maximum absolute value of pitch moment, in. lb. 

Initial pitching moment, in. lb., where I = balance 
number. 

Total pitching moment, in. ib., where I = balance number. 
Pitching weight tares, in. lb., where I = balance number. 
Model yaw angle, degrees. 

Euler yaw rotation angle between primary balance and 
model, degrees. 

Euler yaw rotation angle between secondary balance and 
model, degrees. 

Euler yaw rotation angle between tertiary balance and 
model, degrees. 

Yaw deflection of primary balance, degrees. 

Yaw deflection of secondary balance, degrees. 
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SYMBOL 

PSID3 

PSIDX(I) 

PSIK 


PS [K2 


PSIK3 

PSIR 


PSIS 

PSIU 

R(I,J) 

RGB 

RM 

RM(I,J) 

RMBASE(I) 

RMMAX 

RMO(I) 

RMT(I) 

RMTARE(I) 


NOMENCLATURE 

Yaw deflection of tertiary balance, degrees. 

Deflection yaw angle constants, where X is the balance 
component (A,S,N,R,P,Y) and I = balance number. 

Euler yaw angle to account for knuckle and/or 
primary balance angles in relation to tunnel support, 
degrees. 

Euler yaw angle to account for orientation of undeflected 
secondary balance in relation to primary balance, 
degrees. 

Euler yaw angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 

Eule yaw rotation angle between model (body) axis and 
reference axis, positive in same direction as PSIB, 
degrees. 

Strut yaw angle, degrees 
Tunnel sideflow angle, degrees. 

I'th row and J’th column in rotation matrix. 

Gravity to balance rotation matrix. 

Rolling moment, in. lb. 

Rolling moment, in. lb., where I = balance number and 
J = correction number. 

Base rolling moment, lbs., where I * balance number. 
Maximum absolute value of roll moment, in. lb. 

Initial rolling moment, in. lb., where I = balance number. 
Total rolling moment, in. lb., where I = balance number. 
Rolling -ve» .*t tares, in. lb., where I = balance number. 
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SYMBOI. 

NOMENCLATURE 

SAREA(I) 

Reference area for balance coefficients. Normally wing 
area, sq. in., where I = balance number. 

SF(I,J) 

Side force, lbs., where I = balance number and 
J = correction number. 

SFO(I) 

Initial side load, lbs., where I = balance number. 

SFT(I) 

Total side load, lbs., where I = balance number. 

SFTARE(I) 

Side weight tares, lbs., where I = balance number. 

TAREA 

Axial momentum tare correction term. 

TAREN 

Normal momentum tare correction term. 

TAREP 

Pitching momentum tare correction term. 

TARER 

Rolling momentum tare correction term. 

TARES 

Side momentum tare correction term. 

TAREY 

Yawing momentum tare correction term. 

THEDX(I) 

Deflection pitch angle constants, where X is the balance 
component (A,S,N,R,P,Y) and I = balance number. 

THETA 

Model euler pitch angle, degrees. 

THETAB 

Euler pitch rotation angle between primary balance and 
model, degrees. 

THETAB2 

Euler pitch rotation angle between secondary balance and 
model, degrees. 

THETAB3 

Euler pitch rotation angle between tertiary balance and 
model, degrees. 

THETAD 

Pitch deflection of primary balance, degrees. 

THETAD2 

Pitch deflection of secondary balance, degrees. 

THETAD3 

Pitch deflection of tertiary balance, degrees. 
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SYMBOL 

THETAK 




T 



THETAK2 


THETAK3 

THETAR 


THETAS 

THETAU 

THETAO.(I) 

W 

x 

XBAR(I) 


XICH(I) 

XK 

XKCH(I) 


NOMENCLATURE 

Euler pitch angle to account for knuckle and/or p rimar y 
balance angles in relation to tunnel support, degrees. 
Euler pitch angle to account for orientation of undeflected 
secondary balance in relation to primary balance, 
degrees. 

Euler pitch angle to account for knuckle and/or tertiary 
balance angles in relation to tunnel support, degrees. 
Euler pitch rotation angle between model (body) axis and 
reference axis, positive in same direction as THETAB, 
degrees. 

Strut pitch angle, degrees. 

Tunnel upflow angle, degrees. 

Wind off zero attitude of each balance, degrees, where 
I = balance number. 

Weight tares. 

Distance of center of gravity to balance center, inches. 
Moment transfer distance measured in the body force axis 
system from the balance moment center to the desired 
moment center, positive in the direction of positive model 
thrust, side and normal force respectively, inches, where 
I = balance number. 

Intercept for momentum term, where I = balance 
number. 

Constants used in calculating momentum correction 
terms. 

Slope for momentum term, where I = balance number. 


- 
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r 

SYMBOL 

XREF 

; y 

YBAR(I) 

, YM 

YM(I,J) 

if 

YMBASE(I) 

YMMAX 

YMO(I) 

YMT(I) 

YMTARE(I) 

YREF 

Y 0 

y 3 


y 4 


ZBAR(I) 

ZREF 





NOMENCLATURE 

Moment transfer distance. Measured relative to and in 
the same direction as XBAR, inches. 

Distance of center of gravity to balance center, inches. 

See XBAR. 

A 

Yawing moment, in. lb. 

Yawing moment, in. lb., where I = balance number and 
J = correction number. 

Base yawing moment, lbs., where I = balance number. 

Maximum absolute value of yaw moment, in. lb. 

Initial yawing moment, in. lb., where I = balance 
number. 

Total yawing moment, in. lb., where I = balance number. 

Yawing weight tares, in. lb., where I = balance number. 

Moment transfer distance. Measured relative to and in 
the same convention as YBAR, inches. 

Initial balance loads, side force, lbs. 

Balance component quantity corrected for high 

interactions coupled with high model restraints, side 
force, lbs. 

Balance component quantities corrected for balance 
orientation to gravity axis, side force, lbs. 

Distance of center of gravity to balance center, inches. 

See XBAR. 

Moment transfer distance. Measured relative to and in 
the same convention as ZBAR, inches. 
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APPENDIX D 
Module D 


Balance Loads and Model Attitude 
A. Required Consent* 

Required constants are defined in the nomenclatures. 

1. Primary balance deflection constants - Aangle/Aload 


PSID A 1 

= APSID/AAF( 1 ,3) 


THEDA1 

= ATHETAD/AAF( 1,3) 

See related 

PHID A 1 

= APHID/AAF(1,3) 

items 2. and 3. 

PSIDN1 

= APSID/ANF( 1,3) 


THEDN1 

etc. 

= ATHETAD/ANF( 1,3) 



2. Primary balance deflection angle names - PSID, THETAD, PHID. 
These names are optional as shown in item 3. However, they are 
suggested and extreme care should be used if changed since this is 
based on these names. No values are required for these angles since 
they are computed internally from the constants supplied under 
item 1 as follows: 


PSID = (PSIDA1)AF(1,3) + (PSIDN1)NF(1,3) 


+ (PSIDS1)SF(1,3) + (PSIDR1)RM(1,3) 


+ (PSIDP1)PM(1,3) + (PSID Y1)YM( 1,3) 

(Eq. D-l) 

THETAD = (THEDA1)AF(1,3) + 

(Eq. D-2) 

PHID = (PHIDA1)AF(1,3) + 

(Eq. D-3) 


3. Input of items 1 and 2 - Deflection angle names and constants are 
input from C-card images (which may be modified) stored on 
magnetic storage disks. A maximum of six deflections is permitted. 

Therefore, the six values assigned in the yaw plane (PSD for 
example are PSIDA1, PSIDS1, PSIDN1, PSIDR1, PSIDP1, and 
PSIDY1 as defined in item 1. 

4. Input of rotations from gravity to primary balance - Order of 
rotations from gravity to primary balance axis system (see Figure 
D- 1(a) to D-l(e)) are input from the R-card image names stored on 
magnetic storage disks. 

5. Secondary balance deflection constants - Aangle/Aload 


PSIDA2 

APSID2/AAF(2,3) 


THEDA2 

ATHETAD2/AAF(2,3) 

See related 

PHIDA2 

APHID2/AAF(2,3) 

Items 6. and 7 

PSIDN2 

APSID2/ANF(2,3) 


THEDN2 = 
etc. 

ATHETAD2/ANF(2,3) 



6. Secondary balance deflection angle names - PSID2, PSID3, 

THETAD2, PHID2. These names are optional as shown in item 7. 
However, they are suggested and extreme care should be used if 
changed since this description is based on these names. No values 
are required for these angles since they are computed internally 
from the constants supplied under item 5 as follows: 
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PSED2 = (PSIDA2)AF(2,3) + (PSIDN2)NF(2,3) 

+ (PSroS2)SF(2,3) + (PSIDR2)RM(2,3) (Eq. D-4) 
+ (PSIDP2)PM(2,3) + (PSIDY2)YM(2,3) 


THETAD2 = (THEDA2)AF(2,3) + (Eq.D-5) 

PHID2 = (PHIDA2)AF(2,3) + (Eq. D 6) 


7 . Input of items 5. and 6. - Deflection names and constants are input 
from C-card image names stored on magnetic disks. Six is the 
maximum number of deflections permitted. 

8. Tertiary balance deflection angles are handled in a manner similar 
to primary and secondary balance constants. 

9. Input of rotations (THETAK2, PSIK2, THETAD2, etc.) from pr imar y 
balance to secondary balance - Order of rotations from the primary 
balance to the secondary balance are input from R-card images 
stored on magnetic disks. See Figure D-l(f). 

10. Wind-off-zero attitude of each balance - Input PHIO, THETAO, from 
card images stored on magnetic disks for each balance. This 
option is normally used as a result of problems associated with 
option 2. It is also used when data zeros are not used in the force 
data reduction scheme. If data zeros are not taken and values are 
not input from the disk, PHIO = THETAO = 0 is assumed. See 
Figure D-l(g). 
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F 



11. Weight tares and attitude loads - Tares are determined 

automatically from a 700 series weight-shift run made 
immediately before each model configuration tunnel run. Do not 
input W, X, Y, Z, W(AF), W(SF) etc. 

12. HIRAFI, HIRNFI, HIRSFI, HIRPMI, HIRYMI, HIRRMI where I 
= balance number - These constants correct for the effect of having 
a model with high restraints (HIR) coupled with a balance with 
high interactions (AF, NF, etc.). Thus, the name HIRAFI, 
HIRNFI, etc. These constants are obtained for each balance 
component by the following equation. 

HIR xx(I) = T unne l balance xx calibration _ ^ 

xx span check 

(Eq. D-7) 

where xx = balance component 

Note that when this correction is applied, the balance spans should 
be used in the standard program for quantities (EU) and not in- 
tunnel calibration. These constants are input from the C-card 
images stored on the magnetic disks for each balance. 

13. KSIGN(I) - Constant for correcting balance quantities for 
grounding by the wrong end, where I = balance number. As shown 
in Figure D-2, grounding the balance by the wrong end ("A" cases) 
rather than the taper end results in a change of each balance 
component sign. Therefore 



T. 


V • 
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KSIGN(I) s 1 for normal balance attachment 
KSIGN(I) = -1 for grounding balance by wrong end. 

14. THETAU - Tunnel upflow angle, see Figure D-3. 

15. PSIU - Tunnel sideflow angle, see Figure D-3. 

16. Input of items 14. and 15. - THETAU and PSIU are the required 
rotations for the wind-to-gravity transformation and are input from 
the T-card images (tables as function of MACH) stored on magnetic 
disks. 

17. Euler yaw, pitch and roll rotation angles (PSIB, THETAB, PHIB) 
between balance and model, are shown in Figure D-4(a). 

18. Input of - PSIB(I), THETAB(I), and PHIB(I) - Required rotations for 
the balance-to-model transformation are input from C-card ima ges 
stored on magnetic disks. 

19. XBAR(I), YBAR(I), ZBAR(I) - Moment transfer distances are 
measured in the body force axis system from the balance moment 
center to the desired moment center, positive in the direction of 
positive model thrust, side and normal force, respectively (see 
Figure D-4(b)). Input from C-card images stored on magnetic 
disks, where I = balance number. 

20. ARPB(II,K) - Areas or momentum arms * areas used with 
PBASE(II) for computing base force and base moment tares, where 
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II = orifice number. Use care to insure proper tare force signs. 

Area and arm units must be consistent with units of base pressures 
and balance components. ARB(II,K) is the same but for the second 
balance. ARP(II,K) is the same but for the third balance. 

21. Input of item 20. - Areas and arm x areas are input from C-card 
images stored on magnetic disks. A maximum of 20 may be used. 

22. KPP - Units conversion factor, initialized at 1. 

If PBASE is in PSF and PO is in PSI, KPP = 144 
If PBASE is in PSI and PO is in PSF, KPP = .0069444 
If PBASE is differential (PBASE-PO), KPP * 0 
If PBASE is absolute, KPP = 1 (standard) 

Input from C-card images stored on magnetic disks if not equal 
to 1.0. 

23. Input of items - PSIR, THETAR, and PHIR are the required 
rotations for the model (body) to reference axis transformation and 
are input from C-card images stored on magnetic disks. 

24. XREF, YREF, ZREF - Moment transfer distances are measured 
relative to and in the same convention as XBAR, YBAR and ZBAR. 
Input from C-card images stored on magnetic disks. 

B Test for Balance Loada a nd Model Attitudes 
If NUBAL = 0, skip module D. 
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1. In general, the balance component naming system follows the 
format of WX(Y,Z), where 

WX = component name is as follows: 

AF = Axial force 
NF = Normal force 
SF = Side force 
PM •- Pitching moment 
YM = Yawing moment 
RM = Rolling moment 

Y = balance number associated with component 

1 = 1st balance 

2 = 2nd balance 
etc. 

Z = number of corrections applied to component 
(uncorrected quantity = 1). 


D. Uncorrecfced Balanc e Quantities 

1. Signs on component quantities are uncorrected and thus are a strict 
function of model-balance orientation and the manner in which the 
model-balance attachment is made. Figure D-2 provides sketches 
showing the eight most frequent cases of model-balance orientation 






► . 




2. For normal NASA typo balances, the component quantities are 
obtained directly from the standard program for quantities. The 
balance components for this type of balance are always named as 
follows: 


Axial force -AF(I,1) 
Normal force -NF(I,1) 
Sideforce -SF(I,1) 
Pitch moment - PM(1, 1) 
Yawmoment - YM(1, 1) 
Roll moment -RM(I,1) 


= [FI] 


(Eq. D-8) 


where I = balance number 


3. For TASK type balances, the component quantities are also obtained 
directly from the standard program for quantities (EU), but 
additional equations must be supplied since axial force and rolling 
moment are generally the only two components obtained directly 
with TASK type balances. The following equations and names are 
suggested for the engineering units program. The following 
equations assume the axes origin is at the center of the balance. 



NF(I,1) = N 1(1,1) + N2(I,1) 

(Eq. D-9) 


PM(I,1) = N1(I,1) - N2(I,1) 

(Eq. D-10) 


SF(I,1) a S1(I,1) + S2(I,1) 

(Eq. D-ll) 


' M(I,1) ss S1(I,1) - S2(I,1) 

(Eq. D-12) 


The names shown for the final quantities are mandatory. 




* 


H 
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E. Gravity to Balance Transformat ion Angles 

The tunnel support pitch, roll and yaw angles are used in gravity to 

balance transformations. 

Tunnel Support. Pi tch Angle 

1. The tunnel support pitch angle is THETAS. See Figure D-l(a). 

2. THETAS is computed in the engineering units program. It may be 
obtained from the strut encoder or from a "dangle" meter in the 
model. 

Tunnel Support Roll Angl* 

1. The tunnel support roll angle is PHIS. See Figure D-l(a). 

2. PHIS is obtained from the engineering units program. It is obtained 
from the strut encoder. 

Tunnel Support Y aw Angle 

1. The tunnel support yaw angle is PSIS. See Figure D- 1(a) 

2. PSIS is obtained from the engineering units program. 

F * Ba l ance Quantities Corrected for Int eractions. Weight Tares and 

Momentum Tares 

1. Balance component quantities corrected for interactions are named 
as follows: 
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Axial force 

AF(I,2) 

Normal force 

NF(I,2) 

Side force 

SF(I,2) 

Pitch moment 

PM(I,2) 

Yaw moment 

YM(I,2) 

Roll moment 

RM(I,2) 


= [F2] 


(Eq. D-13) 

2. Balance component quantities corrected for high interactions 
coupled with high model restraints are named as follows: 


Axial force 

AF(I,3) 

Normal force 

NF(I,3) 

Side force 

SF(I,3) 

Pitch moment 

PM(I,3) 

Yaw moment 

YM(I,3) 

Roll moment 

RM(I,3) 


= [F3] 


(Eq. D-14) 

3. Balance component quantities corrected for the attitude loads and 
weight tares are named as follows: 


Axial force 

AF(I,4) 

Normal force 

NF(I,4) 

Side force 

SF(I,4) 

Pitch moment 

PM(I,4) 

Yaw moment 

YM(I,4) 

Roll moment 

RM(I,4) 


= [F4] 


(Eq. D-15) 
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4. Initial balance loads or weight tares are named as follows: where 
I s balance number. 


AFO(I), NFO(I), SFO(I) ' 
PMO(I), YMO(I), RMO(I) 


(Eq. D-16) 


5. Total bidance loads (AF(I,1) + AFO(I), NF(I,1) + NFO(I), etc. are 
named as follows: 


AFT(I), NFT(I), SFT(I) 
PMT(I), YMT(I), RMT(I) 


= (FT] 


(Eq. D-17) 


6. First order interactions are represented by a matrix Cl; second 
order interactions are represented by a matrix Cf\ 

7. Attitude weight tares are named as follows: 

AFTARE(l), NFTARE(l), SFTARE(I) 

PMTARE(I), YMTARE(I), RMTARE(I) 

(Eq. D-18) 

8. Constants required from the project engineer are: 

a. For gravity-to-primary-balance rotations, see Figure D-l(e). For 
gravity- to-tunnel-strut rotation, see Figure D-l(a). 

THETAS, PHIS, and PSIS are supplied from Section E. 

For tunnel strut-to-undeflected-primary balance rotations, see 
Figure D-l(b) and D-l(c). 

PSIK, THETAK, PHIK 


[F'tajre]- 
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For undeflected balance-to-deflected-balance rotations, see Figure 
D-l(d). 

PSIDA1, THEDA1, PHIDA1 
PSroSl, THEDS1, PHIDS1 
PSIDN1, THEDN1, PHIDN1 
PSIDR1, THEDR1, PHIDR1 
PSIDP1, THEDP1, PHroPl 
PSIDY1, THEDY1, PHIDY1 

b. Primary -to- secondary -balance rotations 

For primary balance-to-undeflected-secondary balance 
rotations, see Figure D-l(f). 

PSIK2, THETAK2, PHIK2 

Undeflected secondary balance-to-deflected-secondary balance 
rotations (with respect to primary balance). 

PSIDA2, THEDA2, PHIDA2 
PSIDS2, THEDS2, PHEDS2 
PSIDN2, THEDN2, PH1DN2 
PSIDR2, THEDR2, PHIDR2 
PSIDP2, THEDP2, PHIDP2 
PSDDY2, THEDY2, PHIDY2 

The third to fifth balance is similar to the above but with the 
number 3 to 5 replacing the number 2 in the second balance. 
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For wind-off-zero attitude of each balance (See Figure D-l(a)) 


PHIO, I; THETAO, I, 


c. High restraint and interaction constants 
HIRAFI, HIRNFI, HIRSFI 
HIRPMI, HIRYMI, HIRRMI 


9. The following description on correcting balance quantities for 
interactions and weight tares does not provide the exact equations 
for computing corrected balance quantities. The PAB balance check 
point program or the contractor's user manual must be consulted 
for these. However, this does provide the general outline for 
computing corrected balance quantities. 


Determine uncorrected total loads, 


[FUT] = [FI] + [FO] = 


AF(I,1) + AFO(I) ' 
SF(I,1)+SF0(I) 
NF(1, 1) + NFO(I) 
RM(1, 1) + RMO(I) 
PM(1, 1) + PMO(I) 
YM(I,1) + YMO(I) 


(Eq. D-19) 


Correct for interactions 3 


^ Smith, David L.: An Efficient Algorithm using Matrix Methods to Solve Wind-Tunnel Force- 
Balance Equations. NASA TN D- 6860 , 1972 . 
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a. [FUT1 = [C 1 ]*IFT]+[C 2 ]*IFP] 


(Eq.D-20) 



where [Cj] and [C 2 ] are balance interaction constants 


b. Therefore 


[FT1 = [C 1 ]- 1 •[FUT]-[C 1 ]" 1 *[C 2 1*IFP] 


(Eq. D-21) 


Compute corrected delta balance loads, 


[F2] = [FT] - [FO] = 


AF(I,2) 


’AFT(I) 

SF(I,2) 


SFT(I) 

NF(I,2) 


NFT(I) 

RM(I,2) 


RMT(I) 

PM(I,2) 


PMT(I) 

YM(I,2) 


YMT(I) 


- AFO(I) ' 

- SFO(I) 

- NFO(I) 

- RMO(I) 

- PMO(I) 

- YMO(I) 


(Eq. D-22) 


10. Correct forces and moments for high model restraints coupled 
with high balance interactions 

AF(I) + (HIRAF)AF(I,1) ' 

SF(I) + (HIRSF)SF(I.l) 

NF(I) + (HIRNF)NF(U) 

RM(I) + (HIRRM)RM(I,1) 

PM(I) + (HIRPM)PM(I.l) 

YM(1) + (HIRYM)YM(U) 

(Eq. D-23) 


(F31 = [F2] + K[F1) = 


AF(I,3) 


SF(I,3) 


NF(I,3) 


RM(I,3) 


PM(I,3) 


YM(I,3) 
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11. Depending on the value of the constant KMOM, balance 

components are further corrected for balance/bellcws interactions 
and momentum flow effects. 



\ 


* \ 


l 

■ u 


If KMOM = 0, 

no further balance corrections are applied and equations 
D-26 to D-33 are skipped. 


If KMOM > 0 


'AF' 


AF(I,3) 

SF 


SF(I,3) 

NF 


NF(I,3) 

RM 


RM(I.3) 

PM 


PM(I,3) 

YM 


YM(I,3) 


(Eq. D-24) 


► - 


Balance/bellows interactions and momentum flow effects on the 
balance are computed after high restraint corrections. 


PTZERO 

PTZERO 

PTZERO 


PTANKG if MCODE = 1 or 3 
PTANKH if MCODE = 2 or 4 
PTKSON if MCODE = 5 


DELP = PCH(l) - PTZERO 


(Eq. D-25) 


If PTJZPO is less than or equal to 1.2 then DELP = 0.0 
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Where PTJZPO is the weighted average of the nozzle pressure ratios 
for the primary air system. (Normally air system n umbe r 1.) 


DELSQ 

AREA 

ASQ 

FNO 

PMO 

XNSQ 

XPSQ 


DELP*DELP 
XK(8,3V12.0 
AREA * AREA 
XK(6,3V12.0 
XK(7,3V12.0 
FNO * FNO 
PMO * PMO 


TAREN = XK(1,1) + XK(2,1) * FN + XK(3,1) * PM + XK(4,1) * RM + XK(5,1) * YM 
+ XK(6,1) * SF + DELP * (XK(46,1) + XK(47,1) * FN + XK(48,1) * PM 
+ XK(49,1) * RM + XK(50,1) * YM + XK(51,1) * SF + FNO *XK(52,1) 

+ XNSQ * XK(53,1) + PMO * XK(54,1) + XPSQ * XK(55,1) + AREA * (XK(56,1) 
+ XK(57,1) * FN + XK(58,1) * p M + XK(59,1) * RM + XK(60,1) * YM 
+ XK(61,1) * SF) + ASQ(XK(62,1) + XK(63,1) * FN + XK(64,1) * PM 
+ XK(65,1) * RM + XK(66,1) * YM + XK(67,1) * SF)) 

(Eq. D-26) 

TAREA = XK(7,1) + XK(8,1) * FN + XK(9,1) * PM + XK(10,1) * RM + XK(11,1) * YM 
+ XK(12,1) * SF 

(Eq. D-27) 


Then for PTJZPO greater than or equal to 1.2, the value of PTJZPO is the 

weighted average value of the nozzle pressure ratios for the primary air 
system. 
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TAREA = 


TAREP = 


TAKER = 


TAREY = 


TAREA + XK(37,1) + XK(38,1) * DELP + XK(39,1) * DELSQ + AREA * 
(XK(40,1) 4 XK'41,1) * DELP + XKf42,l) * DELSQ) + ASQ * (XK(43,1) 

+ XK(44,1) * DELP + XK(45,1) * DELSQ) 

(Eq. D-28) 

XK(13,1) + XK(141) * FN + XK(15,1) * PM + XK(16,1) * RM + XK(17,1) * YM 
+ XK(18,1) * SF + DELP * (XK(68,1) + XK(69,1) * FN + XK(70,1) * PM 
4 XK(71,1) * RM + XK(72,1) * YM + XK(73,1) * SF + FNO * XK(74,1) 

+ XNSQ * XK(75,1) + PMO * XK(1,2) + XPSQ * XK(2,2) + AREA * (XK(3,2) 

+ XK(4,2) * FN 4 XK(5,2) * PM + XK(6,2) * RM + XK(7,2) * YM + XK(8,2) * SF) 
4 ASQ * (XK(9,2) + XK(10,2) * FN + XK(11,2) * PM + XK(12,2) * RM 
4 XK03,2) * YM + XK(14,2) * SF)) 

(Eq. D-29) 

XK(19,1) + XK(20,1) * FN 4 XK(21,1) * PM + XK(22,1) * RM + XK(23,1) * YM 
+ XK(24,1) * SF + DELP * (XK(15,2) + XK(16,2) * FN 4 XK(17,2) * PM 
4 XK(18,2) * RM 4 XK(19,2) * YM 4 XK(20,2) * SF 4 FNO * XK(21,2) 

4 XNSQ * XK(22,2) 4 PMO * XK(23,2) 4 XPSQ * XK(24,2) 4 AREA * (XK(25,2) 
4 XK(26,2) * FN 4 XK(27,2) * PM 4 XK(28,2) * RM 4 XK(29,2) * YM 
4 XK(30,2) * SF) 4 ASQ * (XK(31,2) 4 XK(32,2) * FN 4 XK(33,2) * PM 
4 XK(34,2) * RM 4 XK(35,2) * YM 4 XK(36,2) 4 SF)) 

(Eq. D-30) 

XK(25,1) 4 XK(26,1) * FN 4 XK(27,1) * PM 4 XK(28,1) * RM 4 XK(29,1) * YM 
4 XK(30,1) * SF 4 DELP * (XK(37,2) 4 XK(38,2) * FN 4 XK(39,2) * PM 
4 XK(40,2) * RM f XK(41,2) * YM 4 XK(42,2) * SF 4 FNO * XK(43,2) 

4 XNSQ * XK(44,2) 4 PMO * XK(45,2) 4 XPSQ * XK(46,2) 4 AREA * (XK(47,2) 
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+ XK(48,2) * FN + XK(49,2) * PM + XK(50,2) * RM + XK(51,2) * YM 
+ XK(52,2) * SP) + ASQ * (XK(53,2) + XK(54,2) * FN + XK(55,2) * PM 
+ XK(66,2) * RM + XK(57,2) * YM + XK(58,2) * SF)) 


(Eq. D-31) 

TARES = XK(31,1) + XK(32,1) * FN + XK(33,1) * PM + XK04,1) * RM + XK(35,1) * YM 

+ XK(36,1) * SF + DELP * (XK(59,2) + XK(60,2) * FN + XK(61,2) * PM 
+ XK(62,2) * RM + XK(63,2) * YM + XK(64,2) * SF + FNO * XK(65,2) 

+ XNSQ * XK(66,2) + PMO * XK(67,2) + XPSQ * XK(68,2) + AREA * (XK(69,2) 
+ XK(70,2) * FN + XK(71,2) * PM + XK(72,2) * RM + XK(73,2) * YM 
+ XK(74,2) * SF) + ASQ * (XK(75,2) + XK(1,3) * FN + XK(2,3) * PM 
+ XK(3,3) * RM + XK(4,3) * YM + XK(5,3) * SF)) 

(Eq. D-32) 


AF(I,3) ' 


AF(I,3) - TAREA' 

SF(I,3) 


SF(I,3) - TARES 

NF(I,3) 


NF(I,3) - TAREN 

RM(I,3) 


RM(I,3) - TAKER 

PM(I,3) 


PM(I,3) - TAREP 

YM(I,3) 


YM(I,3) - TAREY 


(Eq. D-33) 

12. Perform gravity-to-balance transformations. 


Let[R|] denote specific Euler transformation matrixes 


[P'bal] [^strutl^knucklelRdeflecUonsJ^'g] 

^[ R GBjFg] 


(Eq. D-34) 
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where 


[Ftai] = vector representing balance quantities in balance axis. 
[Fg] = vector representing balance quantities in gravity axis. 
[F gb ] = gravity- to-balance axis transformation matrix. 


13. Determine weight tares (attitude loads) 


AFTARE 

SFTARE 

NFTARE 

RMTARE 

PMTARE 

YMTARE 


w(sin 0g - sin 0 o j 

w(cos 0 g sin <j) g - cos 0 O sin <j> 0 ) 

-w(cos 0g cos <|)g - cos 0 O cos 4> 0 ) 

SFTARE(Z) - NFTARE(Y) 
AFTARE(Z) + NFTARE(X) 
SFTARE(X) + AFTARE(Y) 


(Eq.D-35) 


Correct for weight tares (attitude loads) 


[F4] = [F3] - [FT ARE] = 


AF(I,4) 


AF(I,3) - AFTARE(I) ' 

SF(I,4) 


SF(I,3) - SFTARE(I) 

NF(I,4) 


NF(I,3) - NFTARE(I) 

RM(I,4) 


RM(I,3) - RMTARE(I) 

PM(I,4) 


PM(I,3) - PMTARE(I) 

YM(I,4) 


YM(I,3) - YMTARE(I) 


(Eq. D-36) 


IfKMOM = 1, 

nonblowing balance corrections are applied and FAMOM(I) is 
automatically computed along with the values of PTJ/PO and 
FI, the weighted average values of each air system. 
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APCH = 0.0 
AMOM(I) = 0.0 



FJCON/FI - ftPTJ/PO) Table lookup and linear interpolation. 
FAMOM(D=[AF(I,4)]-FI[FJCON/FI] 


* 

A 


t 


i 



(Eq. D-37) 


The values of FJCON/FI are obtained from an input table which 
results from averaged Stratford choke nozzle data obtained over 
many years. Typical table values are given below: 


PTJ/PO FJCON/FI PTJ/PO FJCON/FI 


1.0 

1.3 

1.5 

2.0 

3.0 

4.0 

4.5 


0.0 

5.0 

0.9700 

0.9820 

6.0 

0.9600 

0.9905 

7.0 

0.9500 

0.9960 

8.0 

0.9*25 

0.9920 

10.0 

0.9300 

0.9815 

0.9760 

14.0 

0.9125 


A maximum of 15 values can be input to the computer as a T table. 

G * Balance Quantities Corrected for Metho d of Attachment 

1. Balance component quantities corrected for method of attachment 
are named as follows: 

[F5] = 


(Eq. D-38) 
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AF(1,5) 

SF(I,5) 

NF(I,5) 

RM(I,5) 

PM(I,5) 

YM(I,5) 


Where I = balance number. 


2. The constant required from the project engineer is KSIGN(I). 


[F5] = KSIGN * [F4] = 


AF(I,5) 


’KSIGN(I) * AF(I,4) 

SF(I,5) 


KSIGN(I) * SF(I,4) 

NF(I,5) 


KSIGN(I) * NF(I,4) 

RM(I,5) 


KSIGN(I) * RM(I,4) 

PM(I,5) 


KSIGN(I) * PM(I,4) 

YM(I,5) 


KSIGN(I) * YM(I,4) 


(Eq. D-39) 


H. Angle of Attack and Sideslip Angle 

1. The following defin itions denote various transformation matrixes 
which are obtained from given orders of Euler rotation angles. 

[Rwg] = v,ind-axis-togravity-axis transformation matrix 


[R G b] = gravity-axis-to-balance-axis transformation 

matrix. This matrix is established from rotation 
angles supplied in section F, therefore 

[R G b] * [Rstrut] [Rknuckle] [Rdeflection] 


[R bm ] = balance-axis-to-model axis transformation 
matrix 


2. The constants required from the project engineer are THETAU, 
PSIU, PSIBI, THETABI and PHIBI. 
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For wind-to-gravity rotation angles, see Figure D-3. 

For balance-to-model rotation angles, see Figure D-4(a). 

The matrix [Rqb]» which transforms a vector in the gravity axis system to 

the balance axis system, may now be computed by a yaw, pitch, and roll 
rotation. 


[ r gb] - 


R ii R12 R 13 
R21 R 22 R 23 

R 31 R32 R33 


[Rgb ] - [Rz (<b)lR y (0)jRx ( V )] 

(Eq. D-40) 


[ r gb] = 


1 

0 0 

COS0 

0 

- sin 0* 

"cosy -siny 

o' 

0 

co8<j> — sin. <{> 

0 

1 

0 

sin y cos \|f 

0 

0 

sin<j> cos<|> 

sin0 

0 

cos 0 

0 0 

1 


cos 0 

0 

-sin 0 

cos y 

-siny 

O' 

-sin<|>sin0 

cos (}> 

- sin 4> sin 0 

siny 

cosy 

0 

cos <J> sin 0 

sin<|> 

cos <|> cos 0 

0 

0 

1 


cos 6 cos \|/ 

-sin <J> sin 0 cos y + cos <{> sin \p 
cos (|> sin 0 cos v + sin <J> sin y 


- sin v cos 0 - sin 0 

sin 4» sin 0 sin y + cos <|> cosy - sin 4> cos 0 
- cos <(> sin 0 sin y + sin $ cos y - cos $ cos 0 


(Eq.D-41) 


where 0 is pitch angle, <J> is roll angle, and y is yaw angle. 
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For wind-to-gravity rotation angles, see Figure D-3. 


For balance-to-model rotation angles, see Figure D-4(a). 

The matrix [Rgb]» which transforms a vector in the gravity axis system to 

the balance axis system, may now be computed by a roll, yaw, and pitch 
rotation. The result is the final rotation matrix from the wind axis to model 
axis. 


[ r wm] = [ r bm][ r gbI r wg] 


W n W 12 W 13 " 
W 21 w 22 w 23 
W 31 W 32 W33 


( R Wm 1 = [ R y(S)l R z(V)][ R x(*)] 

(Eq. D-42) 


[ r wm] 


COS0 

0 

- sin 0 

cosy -siny 

0 ‘ 

'1 

0 

0 ' 

0 

1 

0 

sin y cos y 

0 

0 

cos<> 

-- sin<|> 

sin0 

0 

cos0 

0 0 

1 

0 

sin<)> 

cos<|> 


Performing the matrix multiplcations 



cos 0 cos y 
siny 

sin 0 cos y 


cos 0 sin y cos <j> - sin 0 sin <J> 
cos v cos <j> 

sin 0 sin y cos <J) + cos 0 sin <J> 


cos 9 siny sin<(>-sin 0 cos<() 
-cosy sin <J> 

sin 9 sin y sin$ cos 0 cos<J> 


(Eq. D-43) 


A discussion of these matrices is in Gainer, Thomas G. and Hoffman, 
Sherwood, Summary of Transformation Equations and Equations of Motion 
Used in Free-flight and Wind-tunnel Data Reduction and Analysis, NASA 
SP-3070. 
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Using the definitions shown in Figure D-5 and the above information 


ALPHA = TAN 1 



Note that for <J> = 0\ a = 0 


PSI = SIN (W 21 ) 
BETA = -PSI 

-l 

THETA = SIN (-R 13 ) 


-1 

PHI = TAN 


-Bail 

R33 J 


(Eq. D-44) 


(Eq.D-45) 

(Eq.D-46) 

(Eq. D-47) 
(Eq.D-48) 


I- Body Axis Components; Rotation and Translation fro m Balance-to- 
Model Axis 

1. Balance components rotated to the model (body) axis are named as 
follows: 


Axial - FA(1,I) 

Side - FY(1,I) 

Normal - FN(1,I) 

Roll - MX(1,I) 

Pitch - MY(i,I) 

Yaw - MZ(U) 
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2. Balance components rotated and translated to the model (body) axis 
are named as follows: 


Axial 

- FA(2,I) 

Side 

- FY(2,I) 

Normal 

- FN(2,I) 

Roll 

- MX(2,I) 

Pitch 

- MY(2,I) 

Yaw 

- MZ(2,I) 


3. The constants required from the project engineer are XBAR, YBAR 
and ZBAR. (See Figure D-4.(b).) 

The matrix is used to transform the components in the balance axis to the 
model (body) axis system as follows: 


'FA(il)' 


AF(I,5)" 

FY(1,I) 

= [ r bm] 

SF(I,5) 

FN(1,I) 


NF(I,5) 


and 


-MX(1,I)' 


-RM(I,5)" 

MY(U) 

s 

BQ 

II 

PM(I,5) 

-MZ(U)_ 


-YM(I,5) 


(Eq. D-49) 


or 
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o a 


T 


'FA(l.I) * 


bj 1 AF(I.5)+ b 12 SF(I,5) + b 13 NF(I,5) ' 

FY(1.I) 


b 21 AF(I,5) + b 22 SF(I.5) + b 23 NF(I.5) 

FN(1,I) 


b 31 AF(I,5) + b 32 SF(I,5) + b^NFd.S) 

MX(l.I) 


hi ]RM(1,5) - b 12 pm(1.5) + b 13 YM(I,5) 

MY(1,I) 


-b 21 RM(I,5)+ b 22 PM(I,5)- b 23 YM(I,5) 

MZ(l.I) 


b 31 RM(I,5) - b 32 PM(I,5) + b 33 YM(I,5) 


(Eq. D-50) 


The components are then translated as follows: 


FA(2,I) ' 


FA(l.I) 

FY(2,I) 


FY(l.I) 

FN(2,I) 


FN(1,I) 

MX(2,I) 


MX(l.I) + FN(l.I) * YBAR - FY(1. 1) * ZBAR 

MY(2,I) 


MY(1.I) -FN(1,I)* XBAR - FA(1, 1) * ZBAR 

MZ(2,I) 


MZ(1,I) - FY(1,I)* XBAR - FA(1. 1) * YBAR 


(Eq. D-51) 


J- Pressure Corrections to Body Axis Co mponents 

1. Base and/or cavity pressures are obtained from the standard 
program for quantities and are named PBASE(II). Where 
II = orifice number. 

2. Tunnel static pressure it computed in module A and is named PO. 

3. Base force and moment tares are named as follows: 


Axial - PABASE(I) 

Side - FYBASE(I) 

Normal • FNBASE(I) 

Roll - RMBASE(I) 
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Pitch - PMBASE(I) 

Yaw - YMBASE(I) 

4. Final body axis components, corrected for base tares, are named as 
follows: 

Axial - FA(I) 

Side - FY(I) 

Normal - FN(I) 

Roll - MX(I) 

Pitch - MY(I) 

Yaw - MZ(I) 

Note that axial force is not corrected for internal (duct) axial force. 

5. The constants required from the project engineer are ARPB(II,K) 
and KPP. 

To determine differential base and cavity pressures 

APBASE(II) = PBASE ((II) - [(PO(KPP))], 

(Eq. D-52) 

Noting that a positive differential pressure acting on the base of a model 
causes a thrust, then base pressure force and moment tares are defined as 
follows: 


FABASE(I) = - I [APBASE(II)] * [ARPB(II, U] 

11=1 

(Eq. B-53) 
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FYBASE(I) = - I [APBASEGD] * [ARPB(II,2)] 
H=r 


FNBASE(I) = - I [APBASE(II)] * [ARPB(II,3)] 
ll=i 


RMBASE(I) = - I [APBASE(II)] * [ ARPB(II,4)] 
ll=i 


PMBASEd) = - 2 [APBASE(II)] * [ARPB(II,5)] 
ll=i 


(Eq.B-54) 


(Eq. B-55) 


(Eq. B-56) 


(Eq. B-57) 



1 


YMBASEd) - - I [APBASE(II)]*[ARPB(II,6)] 
ll=i 

(Eq. B-58) 

I 

i 


'FA(I) ' 


FA(2,I) ‘ 


'FABASE(I) ‘ 

FY(I) 

FN(I) 


FY(2,l) 

FN(2,I) 


FYBASE(I) 

FNBASE(I) 

MX(I) 

MY(I) 


MX(2.1) 

MY(2,I) 


RMBASE(I) 

PMBASEil) 

MZ(I) 


MZ(2.I) 


YMBASE(l) 


(Eq. D-59) 


K. Stability Aria Components 

1. Force and moment components in the stability axis are called 
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Drag-FDS(I) 
Side -FYS(I) 
Lift -FLS(I) 
RoU -MXS(I) 
Pitch - MYS(I) 
Yaw -MZS(I) 


where I = balance number. 

Note that drag is not corrected for internal (duct) drag. 

FDS(I) = [FA(I)1 * [COS( ALPHA)] + [FN(I)] * [SIN(ALHPHA)] (Eq. D-60) 

FYS(I) = FY(I) (Eq- D-6D 

FLS(I) = [FN(I)] * [COS( ALPHA)] - [FA(I)] * [SIN(ALPA)] (Eq. D-62) 
MXS(I) = [MX(I)] * [COS( ALPHA)) + [MZ(I>] * [SIN(ALPHA)] (Eq. D-63) 

MYS(I) = MY(I) (Eq. D-64) 

MZS(I) = [MZ(I)] ♦ [COS( ALPHA)] - [MX(I)] * [SIN(ALPHA)] (Eq. D-65) 

L. Wind Aria Comnonenta 

1. Force and moment components in the wind axis are called 


Drag - FD(I) 
Crosswind - FC(I) 


100 


Lift 

- FIXI) 

Roll 

- MXW(I) 

Pitch 

- MYW(I) 

Yaw 

- MZW(I) 


Note that drag is not correct for internal (duct) drag. 

FD(I) = [FDS(D] * [COS(BETA)] - [FYS(I)] * [SIN(BETA)] (Eq. D-66) 

FC(I) = [FYS(I)] * [COSKBETA)] + [FDS(I)] * [SIN(BETA)] (Eq. D-67) 

FL(I) = FIS(I) (Eq. O’ 68 ) 

MXW(I) = [MXS(I)] * [COS(BETA)] + [MYS(I)] * [SIN(BETA)] (Eq. D-69) 

MYW(I) = [MYS(I)] * [COS(BETA)] - [MXS(I)] * [SIN(BETA)] (Eq. D-70) 

MZW(I) = MZS(I) (Eq- D-71) 

M. Alternate Reference Axis Componen ts 

1. Body axis components rotated and translated to an arbitrary 
reference axis system are called 

Axial - FAREF(I) 

Side - FYREF(I) 

Normal - FNREF(I) 

Roll - MXREF(I) 
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Pitch - MYREF(I) 

Yaw - MZREF(I) 

where I = balance number. 

Note that axial force ia corrected for internal (duct) axial force. 

2. The transformation matrix for model axis to reference axis rotations 
is defined as [RmrI- 

3. The constants required from the project engineer are PSIR, 
THETAR, PHIR, XREF, YREF, ZREF and SAREAI where 

I = balance number for model-(body)-to-reference axis rotations. 

4. CAI is from module E. 

The matrix [Rmr] used to transform the components in the model 
(body) axis to a reference axis system as follows: 


FA(I)' = FA(I) - CAI * QO ♦ SAREA(I) (Eq. I>72) 


‘FAREF(I)'" 

JFA(I)'1 

FYREF(I)' 

= [^mr! fy(I) 

FNREF(I)' 

l.FN(l). 


(Eq.D-73) 


and 


‘-MXREF(I)'' 

J-MX(I)' 

MYREF(I)' 

= [Rmr]MY(I) 

-MXREF(I)' 

1-MZ(I). 


(Eq.D-74) 


or 


'FAREFd)' ' 


‘m n FA(iy + m 12 FY(I)+ m 13 FN(I) 


FYREFd)' 


m 2l FAd)' + + m 23 FN(I) 


FNREFd)' 


m 31 FA(iy + m 32 FY(I) + m 33 FN(I) 

(Eq.D-75) 

MXREF(I)' 


m n MX(I) - m 12 MY(l) + m 13 MZ(I) 

MYREFdV 


--m 21 MX(I) + m 22 MY(D - m^MZd) 


MZREF(I)' 


m 31 MX(I) - m 32 MY(I) + m 33 MZ(I) 



The components are now translated as follows: 


’FAREF(I) 


FAREFd)' 

FYREF(l) 


FYREF(I)' 

FNREF(I) 


FNREF(I)' 

MXREF(I) 

— 

MXREF(I)' + FNREF(I)' * YREF - FYREF(I)' * ZREF 

MYREF(I) 


MYREFd)' - FNREF(I)' * XREF - FAREFd)' * ZREF 

MZREF(I) 


MZREF(I)' - FYREF(I)' * XREF - FAREFd)' * YREF 


(Eq. D-76) 


N. Rasa Force a nd Moment, Tare Coefficients 

1. Base force and moment tare coefficients are called 

Axial - CABASE(I) 

Side - CYBASE(I) 

Normal -CNBASE(I) 

Roll - CRMBASE(I) 

Pitch - CPMBASE(l) 

Yaw - CYMBASE(I) 

where I = balance number. 
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2. Free-stream dynamic pressure is defined in module A and is called 
QO. 

3. The constants required from the project engineer are SAREA(I), 
CHORD! I), and BSPAN(I). 


CABASE(I) ' 


’FABASE(I) 

CYBASE(I) 


FYBASE(I) 

CNBASE(I) 

1 

FNBASE(I) 

CRMBASE(I) 

[QO * SAREA(I)] 

[RMBASE(I) / BSPAN(I)] 

CPMBASE(I) 


[PMBASE(I) / CHORD(I) 

CYMBASK(I) 


[YMBASE(I) / BSPAN(I) _ 


(Eq. D-77) 

O. Base Pressure Coefficients 

1. Base pressure coefficients are called CPBASE(II) 


(Eq. D-78) 


CPBASE(II) = -Ua PBASE(II) 
QO 


where II = orifice number. 


P. Model (Bo dy) Axis Coefficients 

1. Model (body) axis coefficients are called 


Axial 

- CA(I) 

Side 

- CY(I) 

Normal 

- CN(I) 

F ” 

- CMX(I) 

Pit- 

- CMY(I) 

Yaw 

- CMZ(I) 
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where I = balance number. 


2. CAI is from module E. 


rcAd) i 


'FA(I) 


*CAI" 

CY(I) 


FY(I) 


0 

CN(I) 

1 

FN(I) 


0 

CMX(i) 

“ [QO*SAREA(I)l 

[MX(I) / BSPAN(I)] 


0 

CMY(I) 


[MY(D / CHORD(I)] 


0 

CMZ(I) 


[MZ(I> / BSPAN(I) _ 


0 


(Eq.D-79) 

Q. Stability A xis Coefficients 

1. Stability axis coefficients are called 

Drag -CDS(I) 

Side - CYS(I) 

Lift - CLS(I) 

Roll -CMXS(I) 

Pitch -CMYS{I) 

Yaw -CMZSU) 

where I = balance number. 

2. CDIS is from module E. 


'CDS(I) 


'FDS(I) 


CAI' 

CYS(I) 


FYS(I) 


0 

CLS(I) 

1 

FLS(I) 


0 

CMXSG) 

~ [QO*SAREA(I)3 

[MXS(I) / BSPAN(I)] 


0 

CMYS(I) 


[MYS(I) / CHOREKD1 


0 

CMZS(I) 


[MZS(I) / BSPAN(I) . 


.0 


(Eq.D-80) 

k 
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R. Wind Axis Coefficients 


1. Wind axis coefficients are named 


Drag 

- CD(I) 

Crosswind • CC(I) 

Lift 

- CUI) 

Roll 

- CMXW(I) 

Pitch 

- CMYW(I) 

Yaw 

- CMZW(I) 

where I = 

balance numb 


2. CDI is from module E. 


’CD(I) 


'FD(I) 


CDI' 

CC(I) 


FC(I) 


0 

CL(I) 

1 

FL(I) 


0 

CMXW(I) 

[QO*SAREA(Dl 

[MXW(I) / BSPAN(I)] 


0 

CMYW(I) 


[MYW(I) / CHORD(I)] 


0 

CMZW(I) 


[MZW(I) / BSPAN(I) . 


0 


(Eq. D-81) 


S. Alternate R eference Axis Coefficients 
1. Reference axis coefficients are named 


Axial 

- CAREF(I) 

Side 

- CYREF(I) 

Normal 

- CNREF(I) 

Roll 

- CMXREF(I) 

Pitch 

- CMYREF(I) 

Yaw 

- CMZREF(I) 
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where 1 = balance number. 


CAREF(I) ' 


’FAREF(I) 

CYREF(I) 


FYREF(I) 

CNREF(I) 

1 

FNREF(I) 

CMXREF(I) 

[QO*SAREA(I)] 

[MXREF(I) / BSPAN(I)] 

CMYREF(I) 


[MYREF(I) / CHORD(I)] 

CMZREF(I) 


[MZREF(I) / BSPAN(I) 


(Eq. D-82) 


T. Miscellaneous Equations 

1. Base drag coefficient is called CDBASE(I). Where I = balance 
number. 


CDBASE(I) = [CABASE(I)] * [COS(ALPHA)] + [CNBASE(I)] * 
[SINXALPHA)] 


(Eq. D-83) 


2. Lift-over-drag ratio in the stability axis is called LS/DS(I). 


LS/DS(I) = CLS(I)/CDS(I) (Eq. D-84) 

3. Lift-over-drag ratio in the wind axis is called L/D(I). 


IVD(I) = CL(I)/CD(I) (Eq. D-85) 

4. Lift coefficient squared is called CLSQR(I). 

CLSQR(I) = [CLS(I)] * [CLS(I)] (Eq. D-86) 
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U. Calculation of Initial Weight Tares and Atti tude Load Constants 

1. The initial weight tares and attitude load constants may be obtained 

by either of three methods for each strain gage balance. 

a. Method I - Data obtained at an arbitrary series of pitch angles (2 £ 
number of pitch angles ^ 30). This method cannot be used with a 
balance without an axial force component. 

b. Method II - Data obtained at an arbitrary series of roll angles (4 £ 
number of roll angles < 30). Normally, the roll angles will be 0*, 
90°, 180", and 270°. The roll angle must be specified in a digital 
channel with name PHIK or with PHIS. (Note that this method 
must be used for balances without an axial force component). 
This method cannot be used with a balance that does not have a 
rolling moment coefficient. 

c. Method III - Data obtained at an arbitrary series of roll and pitch 
angles, (number of angles £ 30) 

V. Calculation of Initial Weight Tares a nd Attitude Load Tares (Strain 

Gage Balance) 

1. Calculate 

a. Ka, 1 = cos(J>.cos0. (Eq. D-87) 

b. K a>2 = 3in0. (Eq. D-88) 

c • K A) 3 = sin0.co80. (Eq. D-8S) 
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2. Determine from balance deck number of components and what these 
components are. 

3. Determine maximum value of each equipment over entire tare run. 


a. 

FNMAX(I) = ABS(NF(I,l))max 

(Eq. D-90) 

b. 

FAMAX(I) = ABS(AF(I,l))max 

(Eq. D-91) 

c. 

FYMAX(I) = ABS(SF(I,l))max 

(Eq. D-92) 

d. 

PMMAX(I) = ABS(PM(I,l))max 

(Eq. D-93) 

e. 

RMMAX(I) = ABS(RM(1 , 1 ) )max 

(Eq. D-94) 

f. 

YMMAX(I) = ABS(YM(I,l))max 

(Eq. D-95) 


4. Initialize initial weight tares and attitude load constants. 

a . Set AA = AN = AY = 0 

Amj = A m 2 = Anj = An 2 = A i\= A ^2 = 0 
x — y = z = 0 

(Eq. D-96) 

b. Assume NFO(I) = AFO(I) = PMO(I) = RMO(I) = YMO(I) = 

SFC(I) = 0. (Eq. D-97) 


5. For each data point correct balance quantities for interactions. 
Determine uncorrected total loads, [FUT] 


|FUT1 = |F1)+(F0] = 


AF(I,1) + AFO(I) 
SF(I,l) + SFO<l) 
NF(I,1) + NFO(I) 
RM(I,1) + RM0(I) 
PM(I.1) + PM0(1) 
YM(I,1) + YMO(I) 


(Same as Eq. D-19) 
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Correct for interactions (see footnote 3) 


a. [FUT] = [C,]*[FT] + [C 2 ]*[FP] (Same as Eq. D-20) 


where [Cj] and [C 2 ] are balance interaction constants 


b. Therefore 


[FT] ^Cj' 1 *[FUT]-[C 1 ] -1 *[C 2 ]*[FP] 


(Same as Eq. D-21) 


Compute corrected delta balance loads, [F2] 


[F2] = [FT]-[F0] = 


[AF(I,2)' 


’AFT(I) 

SF(I,2) 


SFT(I) 

NF(I,2) 


NFT(I) 

RM(I,2) 


RMT(I) 

PM(I,2) 


PMT(I) 

YM(I,2) 


YMT(I) 


- AFO(I) ' 

- SFO(I) 

- NFO(I) 

- RMO(I) 

- PM0(P 

- YMO(I) 


(Same as Eq. D-22) 


Correct forces and moments for high model restraints coupled with 
high balance interactions 


(F3| = |F2] + KIF11 = 


[AF(I.3)' 


AF(I,2) 

SF(I,3) 


SF(I,2) 

NF(I,3) 


NF(I,2) 

RM(I,3) 


RM(I,2) 

PM(I,3) 


PM(I,2) 

YM(I,3) 


YM(I.2) 


+ (HIRAF)AF(I.l) ' 
+ (HIRSF)SF(U) 
+ (HIRNF)NF(I.l) 
+ (HIRRM)RM(U) 
+ (HIRPM)PM(U) 
+ (HIRYM)YM(U) 


(Same as Eq. D-23) 
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6. Determine balance rotation from gravity axis. 

a. Determine rotation matrix for each matrix. See first part of this 
module. 

b. Determine [Rg B ] = product of each individual rotation 



c. Then: 


Rob - 


COS 6 COS y 

-sin <(> sin 0 cos y + cos <j> sin y 
cos $ sin 0 cos y + sin 0 sin y 


- sin y cos 0 ~ sin 0 

sin <j> sin 0 sin y + cos $ cos y - sin cos 0 
+ cos $ sin 0 sin y + sin cos y - cos <(» cos 0_ 


(Eq. D-98) 


d. 


[Rcb]- 


R(l.l) R(1.2) R(1.3) 
R(2.1) R(2.2) R(2,3) 
R(3,l) R(3.2) R(3,3) 


(Eq.D-99) 


e. calculate 


THETA = SIN 
PHI = TAN '( 


(-R(l,3)) 

R(2,3) ^ 

R(3,3)J 


(Same as Eq. D-46 thru Eq. D-47) 
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Calculation of Attitu de Load Constanta bv Method I 

1. Solve following matrix equation using a least squares technique 
(MINFIT 4 routine) for AA. 


(-R(1,3)-K a , 2 ) 1 




( A 3)i 

(-Rd,3)-K A2 ) 2 




( A 3) 2 

• 


AA 


• 

• 




o 

• 




• 

(-R(U)-K A , 2 ) k 


1 




where k is the number of data points £ 30 


(Eq.D-100) 


2. AN = AA + Aw 

AY = AN (Eq. D-101) 

where Aw is obtained from balance interaction deck 


3. If PMMAX(I) > YMMAX(I) and > RMMAX(I) 


a . Solve following matrix equation using least squares technique for 
Ami and Am 2 - 




4 Golub, G. H. and Reinsch, C. Numer: Singular Value Decomposition and Least Squares 
Solutions. Math 14 403-420 (1970) Reprinted in Wilkinson, J. H. and Reinsch, C., Linear Algebra, 
Springer Verlag, Berlin, (1971). 
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J 


(Eq. D-102) 


(K A , 1 -R(3,3)) l (-R(l,3)-K A , 2 ) l 

(K a ,i - R(3,3)) 2 (-RU.3) - K a2 ) 2 

• • 


Air*! 


( m 3)! 

(“ 3)2 

• 

• • 

(K A , 1 -R(3,3)) k (-R(l,3)-K A , 2 ) k 


Am 2 


• 

(”3)k 


, Arrij 

b. x = — -f- 
AN 


(Eq. D-103) 


Amo 

c. z = - 

AA 

d. A / 2 = Am 2 

e. An j = Amj 


(Eq. D-104) 

(Eq. D-105) 
(Eq. D-106) 


f. If YMMAXd) > RMMAX(I) solve the following equation for An 2 
and A^i. 


(-R(l,3)-K A . 2 ) l 




(— Rd.3) — K A 2 ) 




• 


An 2 


• 

(-R(l,3)-K A , 2 ) k 

1 

1 



I n 3 + An!(R(2,3) + K A>3 ))^ 


(Eq. D-107) 


and A/i = An 2 
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g. If RMMAX(I) > YMMAX(I), solve the following equations for 
and An 2 


(R(3.3)-K a ,,) 1 




[+/ 3 + A/ 2 (R(2,3) + 

(R(3.3)-K a ,,) 2 

• 




+/ 3 + A/ 2 (R(2,3) + K 

• 

• 




• 

(R<3,3)-K M ) k 




+/ 3 + A* 2 (R(2,3) + K A.3)] k 


(Eq. D-108) 


and An 2 = A/^ 


h v = +^ (Eq.D-109) 

n * y AA 

4. If YMMAX(I) > PMMAX(I) and > RMMAXd) 

a. Solve following matrix equation using a least square technique 
for Anj and An 2 


j(-R(2,3)-K A 3) ( 

(-R(L3)-K a2 ) 1 




(n 3 ) 1 

(-R(2,3)-K Ai3 ) 2 

(-R(1,3)-K a . 2 ) 2 




( n 3) 2 

• 

# 


An t 



• 

• 

• 




• 

• 

• 


An 2 


• 

(-Rtta-K^ 

(-RO,3)-K A , 2 ) k 




(»3) V 


(Eq. D-110) 
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(Eq.D-111) 


(Eq. D-112) 

(Eq.D-113) 
(Eq.D-114) 

f. If PMMAX(I) > RMMAX(I), solve the following equation for Am 2 
and A/ 2 

|(-Ra3)-K A , 2 ) l 
(-R(1,3)-K A2 ) 2 

(-R(l,3)-K A2 ) k 

(Eq.D-115) 


km 2 


m,-Am 1 p A>1 -R(3.3^ l 

m 3 -Am 1 (K A>1 -R(3,3))] 2 


lm 3 - - RO.S))^ 


b. x = + 


A Til 
AY 


An 2 

c - y=+ M 

d. A/j = An 2 

e. Am 2 = Ani 


and A / 1 = Am 2 

g. If RMMAX(I) > PMMAXd), solve following equations for A/ 2 and 
Am 2 
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(Eq.D-116) 


to2,3)-K A , 3 ) t 



[ + < 3 -/ 1 (R<3,3)-K a , 1 )] i 

(-R(1,3)-K a , 2 ) 2 




+*3 - / 1 (R(3,3) - K Atl )] 2 

• 


A^2 



• 

• 




• 

(-R(l*3)-K A>2 ) k 




• 

[+/ a -/ l (R(3.3)-K A , 1 )] k 


and Am 2 = A 


h. 


Z — + 


Am 2 

AA 


(Eq.D-117) 


5. If RMMAX(I) > PMMAX(I) and > YMMAX(I) 

a . Solve the following matrix equation using a least squares 
technique for Al 1 and Al 2 - 


(R(3.3)-K a .,) 2 

• 

-(RtiSl+K^), 

-(R(2.3)+K a . 3 ) 2 

• 


A/i 


+ ( / 3)i 

+ (^) 2 

• 

• 

• 

(R(3.3)-K A .,) k 

• 

• 

-(R(2.3l+K A , 3 ) k 


m 2 


• 

• 

+(^3) k 


(Eq. D-118) 

b. y = ^ (Eq.D-119) 

AN 


c. 


z- 



AY 


(Eq. D-120) 
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d. An 2 * A/j 


(Eq.D-121) 


6 . Am 2 = At 2 


(Eq.D-122) 


f. If PMMAX(I) > YMMAX(I), solve the following equation for Amj 
and An^ 


|(K a , 1 -R(3,3 )) 1 

(K a , 1 -R(3,3 )) 2 





1 m 3 - Am 2 (-R(l,3) - K A 2 )^ 
m 3 -Am 2 (-R(l,3)-K A>2 )] 2 

[m 3 - Am 2 (-R(l,3) - K A 2 )] k 


(Eq.D-123) 


and An! = Amj 


g. If YMMAX(I) > PMMAXd), solve the following equations for An^ 
and Amj. 


-^( 2.31 + ^), 
-{R(2,3) + K Ai 3 J 1 


|An 2 


WR(2.3) + K A3 ) k 


[n 3 -An 2 (-R(l,3)-K A _ 2 )] 1 " 
[n 3 -An 2 (-R(1.3)-K A . 2 )] 2 

[n 3 -An,(-R(1.3)-K A . 2 )] k 

(Gq.D-124) 


and Ami = Ani 
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(Eq.D-125) 


h. 


x = 



AN 


X. Calculation of Attitude Lo ad Constants hv Method II 

1. If FNMAX(I) > FYMAXd), solve the following equations for AN and 

AY. 


k-KM), 

1 

(Ns), 

(K M -R<3,3)) 2 




(Ns), 

• 


AN 


• 

• 



• 

• 




• 

• 




• 

i - R(3,3)) k 




(Ns), 


(Eq.D-126) 


and AY = AN 


2. If FYMAXd) > FNMAX, solve the following equations for AY and 
AN. 


-(R(2,3) + K*, 3 ) 1 

. I 

'(*•),' 

-{R<2,3) + K a , 3 ) 1 




(Ys), 

• 


AY 


• 

• 


= 

• 

• 




• 

• 




• 

-{RM + K*^ 




(Ys), 


(Eq. D-127) 


and AN = AY 
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3. AA = AN - Aw 


(Eq.D-128) 


4. Determine Amj, Am 2 , Anj, An 2 , A l j, A^ 2 » x> y> and z by calculation 
procedure given in Subsection W., item 3. 

Y. Balances Without Six Components 

1. For balances that do not have six components, set the appropriate 
attitude tare constant to zero as indicated below. 

a. If balance does not have a normal-force component: AN = 0 

b. If balance does not have an axial-force component: AA = 0 

c . If balance does not have a side-force component: AY = 0 

d. If balance does not have a pitching-moment component: 

Am} = Am 2 = 0 

e. If balance does not have a rolling-moment component: 

a* 1 = a/ 2 =o 

f. If balance does not have a yawing moment component: 

An 3 = An 2 = 0 

Z. Initial Weight Tare Calculations 


1. Calculate initial weight tares 



a. N 0 = -ANK Atl 

NFO 

(Eq. D-129) 

b. A 0 = -AA K A 2 

AFO 

(Eq. D-130) 
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c. 

m 0 = -AmjKAi + Am 2 K A 3 

PMO 

(Eq.D-131) 

d. 

/ 0 a Af jKa,1 + 

RMO 

(Eq.D-132) 

e. 

no = AnxK^s + An 2 K A ,2 

YMO 

(Eq.D-133) 

f. 

y 0 * ayk A(3 

SFO 

(Eq. D-134) 


AA. New Values of Initial Weight Tares 

1. Go to Subsection V., item 5. and repeat the calculation using new 
values of initial weight tares. Repeat iteration procedure until 
initial weight tares repeat to following accuracy. 


New -Old 
New 


<0.005 


(Eq. D-135) 


BB. Point Calculations 

1. For each point, calculate: 

a. N 4 = N s -[aN(K A j-R( 3,3))] (Eq.D-I36) 

b. A 4 = A s -[aA(-R( 1,3)-K a , 2 )] (Eq.D-137) 

c. m 4 = m, - {[Am,(K A-1 -R(3.3))]-[Am s (Ra3) + K A ^)]} 

(Eq.D-138) 
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d. l,=t 3 - {[ A*,(R<3.3) - K A1 )] - [a/ 2 (R<2,3) + K A3 )| 


e. n 4 — a 3 


f. Y 4 = Y 3 


(Eq. D-139) 


{[An, (Ri2.^)+ K a>3 )] -[An 2 (R(1.3) + K A>2 )]} 

(Eq. D-140) 


[ AY(R(2,3) + K A 3 )] 


(Eq. D-141) 
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THETAS (© s ) is measured in the 
tunnel or gravity X-2 plane. 

(*s - 0). 

©sis generally termed tunnel 
pitch angle. 

PHIS (4> s ) is measured in 
the tunnel or gravity Y~Z plane. 
4>s is generally termed tunnel 
roll angle. 

PSIS (Ys) is measured in 
the tunnel or gravity X-Y plane. 
X s is generally termed tunnel 
yaw angle. 


(a) Gravity to tunnel support axes. 


figure 0-1. Definition of gravity and balance axes showing 

positive directions and rotation angles for gravity 
to balance transformations. 
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•jr.nel support to undeflected balance axes. 
Figure 0-1. Continued. 








(') Illustration of primary balance to undeflected 
secondary balance rotations. 

figure D--1. Continued. 







Aero mode Jet (ofterbody) mode! 



(a) Caee 1, Normal balance arrangement. 
Figure D-2. Model-balance orientation stings. 
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fb) Case 1 A, Case 1 held by opposite end. 




(c) Coee 2. Balance rolled 90° (clockwise). 
Fiaur* D-2. Continued. 




(ofUrbody) mod#! 



(<j) Casa 2A. Same at cast 2 held by opposite end. 
Figure 0—2. Continued. 





Jet (efterbody) mod«i 



(f) Case 3A* Cass 3 held by opposite end. 
Figure D-2. Continued. 



J«t (cfUrbody) mod# I 



(g) Com 4, Balance rolled 90 ° (counterclockwlee). 
Figure D-2. Continued. 
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Bolance yawed 180° or (pitched 180° ond rolled 180°) - reversed. 
Figure D— 2. Continued. 








f k) Case 6 Balance yawed 180° and rolled 90° clockwise (reversed and rolled 90° 
clockwise) or pitched 180° and rolled 90° counterclockwise. 

Figure D-2. Continued. 



Jet (afterbody) model 



(I) Coee &K Cose 6 held by opposite end. 
Figure D-2. Continued. 









Aero mo<M , Jet <<*•*<><*) model 



(n) Coee 7A, Same oe C ate 7 hold by opposite end. 
Figure D-2. Continued. 
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(ofterbody) model 



(p) Cose 8A, Cose 8 held by opposite end. 
Fiour* D— 9 
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MODULE E INTERNAL DRAG (OR EXIT-FLOW DISTRIBUTIONS) 



I 


I 


SYMBOL 

AEXITl 

AEXIT2 

ARAKE(I) 

CAI 

CAI1 

CAI2 

CDI 

CDIS 

CLI 

CNI 

CYI 

FTMDOT1 

FTMDOT2 

FTPR1 

Frpii2 

INDX(I.J) 


IRAKE 


NOMENCLATURE 

Exit areas for duct 1, sq. in. Not required for IRAKE = 2 or 3. 

Exit areas for duct 2, sq. in. Not required for IRAKE = 2 or 3. 

Exit area assigned to each rake total pressure PROBE(I), sq. in. 
Not required for IRAKE = 2 or 3. 

Total internal axial force coefficient. 

Internal axial force coefficient for duct 1. 

Internal axial force coefficient for duct 2. 

Total internal drag coefficient in the wind axis. 

Total internal drag coefficient in the stability axis. 

Total internal lift coefficient. 

Total internal normal force coefficient. 

Total internal side force coefficient. 

Mass flow rate at exit of duct 1, slugs/sec. 

Mass flow rate at exit of duct 2, slugs/sec. 

Ratio of nozzle exit total pressure to free stream static pressure 
for duct 1. 

Ratio oi nozzle exit total pressure to free stream static pressure 
for duct 2. 

Table of values used to assign rake total pressures to specific 
static pressures, where 1 = static pressure probes assigned to 
J = table position. Not required for IRAKE = 2 or 3. 

RAKE code. 

= 0, set CAI=CDIS=CDI=0.0 and skip module 5. 

= 1, computes internal drag. 


1 


:« 
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SYMBOL 


KPR(I) 


MEXIT1 

MEXIT2 

MODOT1 

MODOT2 

M/M01 
M/MO 2 
NPR1 

NPR2 

NPTR1 


NOMENCLATIVE 

= 2, measures exit flow distribution only. 

= 3, obtains internal drag from a given table. 

* 4, obtains internal drag from a row of internal static 
pressures. 

= 5, inlet distortion with rotating rake. 

= 6, inlet distortion with nonrotating rake. 

Needed to correct for bad rake static pressure probes, where 
I ® static pressure probe. Not required for IRAKE = 2 or 3. If no 
correction is to be made to the pressure probe, then its value is set 
to 1.0. If the probe is faulty or does not exist, then its value is set to 
0 . 0 . 

Average exit mach number for duct 1. 

Average exit mach number for duct 2. 

Mass flow rate based on free-stream conditions for duct 1, 
slugs/sec. 

Mass flow rate based on free-stream conditions for duct 2, 
slugs/sec. 

Mass flow ratio for duct 1 . 

Mass flow ratio for duct 2. 

Number of static pressure probes on the rake for duct 1. 

Maximum of 10. Not required for IRAKE = 3. 

Number of static pressure probes on the rake for duct 2. 

Maximum of 10-NPR1. Not required for IRAKE = 3. 

Number of total pressure probes on the rake for duct 1. Maximum 
of 50. Not required for IRAKE = 3. 
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SYMBOL 

NPTR2 

pm/m) 

pi)2/i y ro 

PRAKE(I) 

1>R/PT0(I) 

PS1N1 

PSIN2 

prm/i y ro 

i y ri)2/i y ro 

PTRAKE(I) 

PTR/PTOU) 

KHOE 
SC API 

SCAP2 


NOMENCLATURE 

Numbe r of total pressure probes on the rake for duct 2. Maximum 
ofSO-NPTRl. Not required for IRAKE = 3. 

Ratio of the average duct static pressure to free-stream total 
pressure for duct 1. 

Ratio of the average duct static pressure to free-stream total 
pressure for duct 2. 

Rake static pressure, where I = probe number. 

Ratio of rake static pressure to free-stream total pressure, where 
- probe number. 

Thrust axis yaw angle (degrees) for duct 1. Not required for 
IRAKE = 2 or 3. 

Thrust axis yaw angle (degrees) for duct 2. Not required for 
IRAKE = 2 or 3. 

Ratio of the average duct total pressure to free-stream total 
pressure for duct 1. 

Ratio of the average duct total pressure to free-stream total 
pressure for duct 2. 

Rake total pressure, where I = probe number. 

Ratio of rake total pressure to free-stream total pressure, where 
1 = probe number. 

F roe- stream density. 

Inlet capture area for duct 1, sq. in. Not required for IRAKE = 2 or 
3. 

Inlet capture area for duct 2, sq. in. Not required for IRAKE = 2 or 
3 
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SYMBOL 

THETAN1 

THETAN2 


NOMENCLATURE 

Thrust axis Euler pitch angle (degrees), with respect to body axis 
for duct 1. Not required for IRAKE = 2 or 3. 

Thrust axis Euler pitch angle (degrees), with respect to body 
axis for duct 2. Not required for IRAKE = 2 or 3. 
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APPENPIX-E 

Module E 


Internal Drag (or Exit Flow Distributions) 


A. Required C onstants 

The constants for internal drag calculations are given in the nomenclatures. 
All constants are initialized to a value of 0.0. 


1. IRAKE - 
where 


Rake code 
IRAKE = 

IRAKE = 
IRAKE = 
IRAKE = 
IRAKE = 

IRAKE = 
IRAKE = 


0, Set CAI = CDIS = CDI = 0.0 and skip this 

module. 

1, compute internal drag. 

2, measure exit flow distribution only. 

3, obtain internal drag from a given table. 

4, obtain internal drag from a row of 

internal static pressures. 

5, inlet distortion with rotating rake 

6, inlet distortion with nonrotating rake 


NPTRl 

I 


1=1 


ARAKE(I) = total exit area for duct 1 


(Eq. E-l) 


NPTR2 

I 


ARAKE(I) = total exit area for duct 2 


I=N1 J TR1+1 


(Eq.E-2) 


2. SCAP1, SCAP2 - inlet capture area, where SCAP1 is for duct 1 and 

SCAP2 is for duct 2. Not required for IRAKE = 2 or 

3. 

3. AEXIT1, AEXIT2 • exit areas for ducts 1 and 2, respectively. Not 

required for IRAKE = 2 or 3. 




4. PSIN1, PSEN2 - 


W 






Thrust axis yaw angle, with respect to body axis, 
for ducts 1 and 2, respectively. Positive direction is 
shown on Figure E-l. Not required for IRAKE = 2 
or 3. 


5. THETAN1, THETAN2 - Thrust axis Euler pitch angles, with respect to 

body axis, for ducts 1 and 2, respectively, deg. 
Positive direction is shown on Figure E-l. Figure 
E-l also gives relations to obtain the Euler angle if 
not known directly. Not required for IRAKE = 2 or 
3. 


6. AREF • Model reference area used for coefficients, sq. in. 

If Module B or C is used, this constant is already 
specified. Not required for IRAKE = 2 or 3. 


B. Teat for Module E Computations 

IF IRAKE * 0, skip module E. 

IF IRAKE = 3, do section T only. 
IF IRAKE = 4, do section U only. 
IF IRAKE = 5 or 6, do Module I 


C. Rake Total Pressure 

1. Rake total pressures are called PTRAKE(I). Note that provisions are made 
to survey two exits at one time; however probes are numbered consecutiv sly 
(max. of 50). For example, probes in the first exit may be numbered 1 


P 




-r 

-m . 


157 



through 30; probes in the second exit must start with number 31. 

Where I = probe numbei . 

2. The ratio of rake total pressure to free-stream total pressure is called 
PTR/PTCKI), where PTO is from module A. 

3. The constants required from the project engineer are NPTR1, and NPTR2. 
Calculate PTR/PTO(I) for I = 1, NPTR1 + NPTR2 

PTR / PTOfl) = F1ME| " 

PTO 

(Eq. E-3) 


D. Rake Stati c Pressures 

1. Rake static pressures are called PRAKE(I). Comments C.l. above apply 
except that the maximum number of probes is 10. 

2. The ratio of rake static pressure to free-stream total pressure is called 
PR/PTOG). PTO is from module A. 

3. The constants required from the project engineer are NPR1, and NPR2. 
If NPR1 = 0, skip this part. 

Calculate PR/PTOG) for I = 1, NPR1 + NPR2 


PR / PTO(I) = 


PTRAKEG) 

PTO 


(Eq. E-4) 
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E. Rake Total Pressure/Static Pressure Assignments 

1. If internal drag is to be computed, the project engineer must assign 
specific total pressure measurement to each static pressure 
measurement. 

INDX(I,J) This is the index of the total rake pressure PTRAKE that goes 
with the static rake pressure, PRAKE. The array INDEX has 
two indices, I,J. The value of INDEX(I,J) contains the index of 
the total pressure that goes with the Ith static pressure. The 
value of J starts at one (1) and increments by 1 until the value 
of INDEX(I,J) is zero or if other values of I or J exceed limits. 

Example: To assign PRAKE3 to PTRAKE 24 

INDEX(3,J) = 24 where J is table position 
starting at 1. 

2. For example: 

1 = 1, 2,3,4 
1 = 5,9, 11 
I = 6, 7, 8, 10 

I(Max) = NPTR1 + NPTR2 

3. The constants required from the project engineer are from the I, J table. 

If IRAKE = 2, skip this section. 


, 1=1 
J = 2 
J = 3 


b * 


J = NPR1 + NPR2 


) 
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F. Pact Flow Static-to-Total-Pressure Ratio 

1. The ratio of duct flow static pressure to duct flow total pressure is called 
PR/PTR(J,I), where J and I are the combinations supplied in section E 
above. For the example shown in E., values of PR/PTR(J,1) are obtained 
for: 

PR/PTRl.l 

PR/PTR1.2 

PR/PTR1,3 

PR/PTR1.4 

PR/PTR2.5 

PR/PTR2,9 

PR/PTR2.11 

PR/PTR3,6 

PR/PTR3.7 

etc. 


If IRAKE = 2, skip this section. 

Do the following calculation for J = 1, NPR1 + NPR2 

Do the following calculation for I = those values assigned 


PR/PTR(J,I) 


PR / PTO( J) 
PTR / PTO(I) 


(Eq. E-5) 


G. Correct for Supersonic Duct Mach Numbers 

1. Local duct Mach number is called MD(I). Where I =■ total pressui'e probe 
number on which local Mach number is based. 

If IRAKE = 2, skip this section. 
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Do the following calculation for J = 1, NPR1 + NPR2 
Do the following calculation for I = those values assigned 

If PR/PTR(J,I) < 0.5283, calculate MD(I) using the Newton Raphson method 
with an initial assumption of MD(I) = 1.0001 and correct the total pressure 
ratio for normal shock. 


MD(I) = 



r 2 i 

5/7 

5 . 

7MD(I) -1 

JPTR/PTO(I)f /7 

I 6 

6 

L PR / PTO( J) J 


(Eq. E-6) 


PTR / PTO(I) = PR / PTO(J) 


1 + 


MD(I) 


2 x7/2 


(Eq. E-7) 


H. Compute Subsonic Duct Mach Numbers 

1. This calculation is made for those I, J combinations for which PR/PTR 
(J, I) > 0.5283. 


If IRAKE = 2, skip this section. 


Do the following calculation for J = 1, NPR1 + NPR2 
Do the following calculation for I = those values assigned 

If PR/PTR(J,I) > 0.5283, calculate MD(I) 


MD(I) = ,5* 


PTR/PTCKI) 
PR / PTO( J) J 


i2/7 


-5 


(Eq. E-8) 


161 


I. Compute A verage Duct Pressure Ratios 

1. The ratio of the average duct total pressure to free-stream total pressure is 
called PTD1/PTO for duct 1 and PTD2/PTO for duct 2. 

2. The ratio of the average duct static pressure to free-stream total pressure 
is called PD1/PTO for duct 1 and PD2/PTO for duct 2. 


3. 


The c mstanls required from the project ngineer are ARAKE(I), KPRd), 
NPTR1, NPTR2, NPR1, and NPR2 


PTDl/PTO = 


NPTRl 

I 

1=1 


ARAKE(I)[PTR / PTO(I)] 

NPTRl 


I ARAKE(I) 


i=i 


NPTRl 

If I ARAKE(I) = 0.0, then PD1/ PTO = 10 

i=i 


(Eq. E-9) 


PDl/PTO = 


npri _ 

I KPRdKPR / PTO(I)] 

i=i 

RPin 

I KPR(I) 

i=i 


(Eq. E-10) 


If 


NPTRl 

I 


KPR(I) = 0.0, then PD1/ PTO = L0 


i=i 


If NPTR2 = 0.0, skip equations E-ll and E-12. 

NPTR2 

I ARAKE(I*PTR/PTO(I)l 

'TD2 / PTO l =NPTRltl 

I ARAKE(l) 

I=NPTR1+1 

(Eq. E-ll) 


162 


If 


NPTR2 

z 


ARAKE(I) = 0.0, then PD2 / PTO = 10 


I-NPTRl+1 


NPR2 

I KPR(I)[PR/PTO(I)l 
PD2 / PTO = i=npri+i 

ru^/riKj- NpR2 

I KPR(I) 

I=NPR1+1 

(Eq. E-12) 


If 


NPR2 

I 


KPR(I) = 0, then PD2 / PTO = L0 


I=NPRl+l 


J. Mass-Flow Rates 

1. Mass-flow rate at the duct exit is called FTMDOT1 for duct 1 and 
FTMDOT2 for duct 2. 


2. Mass-flow rate based on free-stream conditions is called MODOT1 for duct 
1 and MODOT2 for duct 2. 


3. TTO, MACH, and PO come from the tunnel parameters, module A. 

4. The constants required from the project engineer are ARAKE(I), SCAP1, 
SCAP2, NPTR1, and NPTR2. 


If IRAKE = 2, skip equations E-13, E-14, E-15 and E-16. 


FTMDOTl = 


0.028563 .nftri 
VTTO + 459.67 Ui 


ARAKE(I) * PRAKE(J) * 


PTR/PTOd) ] 177 
. PR/PTO(J) . 


*MD(I) 


(Eq. E-13) 


where J corresponds to I from E.l. above. 
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(0.028563) * (SC API) * (MACH) * (PO) . 
MUL)U 11= >. ■ - , - 

VTTO+ 459.67 


2 

1+0.2(MACH) 



(Eq.E-14) 


If NPTR2 = 0, skip equations E-15 and E-16. 


(0.028563) 


NPTR2 


FTMODOT2 = * I 

VTTO + 459.67 i= 


NPTR1+1 


ARAKE(I).PRAKE(J).[^^^ 


ll/7 


*MD(I) 


(Eq. E-15) 


where J corresponds to I from E.l. above. 


MODOT2 = 


(0.028563) * (SCAP2) ♦ (MACH) * (PO) , 
VTTO + 459.67 ~~~ 


2 

1+0.2(MACH) 



(Eq. E-16) 


K. Mass-Flow Ratio 

1. M ass -flow ratio for duct 1 is called M/M01. Mass-flow ratio for duct 2 is 
called M/M02. 


If IRAKE = 2, skip the remainder of this section. 


M/M01= 


FTMDQT1 

MODOT1 


(Eq. E-17) 


If MODOTl = 0.0, M/M01 = 0.0 


If NPTR2 = 0, skip equation E-18. 


M / M02 = 


FTMD0T2 

M0D0T2 


(Eq. El- 18) 


1. Free-stream velocity is called VO. 


2. TTO and MACH are from module A. 


VO = 


49.021179VTTOT459.67 
Vl + 0.2(MACH) 2 


(Eq. E-19) 


1. The average exit Mach number for duct 1 is always called MEXIT1. The 
average exit Mach number for duct 2 is always called MEXIT2. 

If IRAKE = 2, skip the remainder of this section. 

f77PTDl/PTOf' 7 c 


MEXIT1=J5* 


PD1/PT0 


If NPTH2 = 0, skip equation E-21. 


(Eq. E-20) 



c PTD2 / PTOl _ 

MEXIT2 = * 5 * : - 5 

1 PD2 / PTO 


\*\ 


N. Internal Axial Force 


1. The internal axial force is called All and AI2 for ducts 1 and 2, 
respectively. 

2. The internal axial force coefficient is called CAI1 and CAI2 for ducts 1 and 
2, respectively. 

3. The total internal axial force coefficient is called CAI. 

4. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI and THETA are from the balance and weight tare calculations, module 
D. Positive directions for PSI and THETA are shown on Figure 

E- 2. 

6. The constants required from the project engineer are AEXIT1, AEXIT2, 
PSIN1, PSIN2, THETAN1, THETAN2, AREF, and NPTR2. 


If IRAKE = 2, skip the remainder of this section. 


All = [(FTDOT1) * VO * COS(PSI) * COS(THETA)] 

- {[1.4 * (PD1/PTO) * PTO * (MEXIT1) 2 ] + [((PD1/PTO) * PTO) - PO]} 
* (AEXIT1) * COS(PSINl) * COS(THETANl) 


CAI1 = 


All 

(QO) * (AREF) 


(Eq. E-22) 


(Eq. E-23) 


CAI = CAI1 


(Eq. E-24) 
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If NPTR2 = 0, skip equations E-25, E-26 and E-27. 



i 





AI2 = [(FTMDOT2) * VO * COS(PSI) * COS(THETA)] 

- {[1.4 * (PD2/PTO) ♦ FrO * (MEXIT2) 2 ] + [((PD2/PTO) * PTO) - PO]} 

+ (AEXIT2) * COS(PSIN2) * COS(THETAN2) 

(Eq. E-25) 


CAI2 = 


AI2 

(QO)*(AREF) 


(Eq. E-26) 


CAI = CAI1 + CAI2 


(Eq. F-27) 


O. Internal Normal Force 

1. The internal normal force is called Nil and NI2 for ducts 1 and 2, 
respectively. 

2. The internal norma! force coefficient is called CNI1 and CNI2 for ducts 1 
and 2, respectively. 

3. The total internal normal force coefficient is called CNI. 

4. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI, THETA, and PHI are from the balance and weight tare calculations, 
module D. Positive directions are shown on Figure E-2. 

6. The constants required from the project engineer are AEXIT1, AEXIT2, 
THETAN1, THETAN2, AREF, and NPTR2. 
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If IRAKE = 2, skip the remainder of this section. 


Nil = {(FTMDOTi) * VO * |COS(PHI) * SIN(THETA) * COS(PSI) + SIN(PHI) + SIN(PSI)]} 
+{[1.4 * (PD1/PTO) * PTO * (MEXIT1) 2 ] + [((PD1/PTO) * PTO) - PO]} 

* (AEXIT1) * SIN(THETANl) 

(Eq. E-28) 


CNI1 = 


Nil 

(QO)*(AREF) 


(Eq. E-29) 


CNI = CNI1 


(Eq. E-30) 


If N?TR2 = 0, skip equations E-31, E-32 and E-33. 


NI2 = {(FTMDOT2) * VO * [COSd'HI) * SIN(THETA) * COS(PSI) + SIN(PHI) + SIN(PSI)]) 
+{[1.4 * (PD2/PTO) * PTO * (MEXIT2) 2 | + [((PD2/PTO) * PTO) - PO]} 

* (AEXIT2) * SIN(THETAN2) 

(Eq. E-31) 


CNI2 = 


NI2 

(QOWAREF) 


(Eq. E-32) 


CNI = CNI1 + CNI2 


(Eq. E-33) 
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P. Internal Side Fores 


1. The internal side force is called YI1 and YI2 for ducts 1 and 2, 
respectively. 

2. The internal side force coefficient is called CYI1 and CYI2 for ducts 1 and 
2, respectively. 

3. The total internal side force coefficient is called CYI. 

4. PTO, PO, and QO are from the tunnel parameters, module A. 

5. PSI, THETA, and PHI are from the balance and weight tares calculations, 
module D. 

6. The constants required from the project engineer are AEXIT1, AEXTT2, 
THETAN1, THETAN2, PSIN1, PSIN2, AREF, and NPTR2. 

If IRAKE = 2, skip the remainder of this section. 

YI1 = {(FTMDOT1) * VO * [SIN(PHI) * SIN(THETA) * COS(PSI) - COS(PHI) * SIN(PSI)]} 
+{(1.4 * (PD1/PTO) * PTO * (MEXIT1)2] + [((PD1/PTO) * PTO) - PO]} 

* (AEXIT1) * COS(THETANl) * SIN(PSINl) 

(Eq. E-34) 


CYI1 = 


YI1 

(QO) * (AREF) 


(Eq. E-35) 


CYI = CYI1 


(Eq. E-36) 


If NPTR2 = 0, skip equations E-37, E-38 and E-39. 


YI2 = {(FTMDOT2) * VO * ISIN(PH1) * SIN(THETA) * COS(PSI) - COS(PHI) * SIN(PSI)]) 
+{[1.4 * (PD2/PTO) * PTO * (MEXIT2) 2 ] + [((PD2/PTO) * PTO) - PO]} 

* (AEXIT2) * COS(THETAN2) * SIN(PSIN2) 

(Eq. E-37) 


CYI2 = 


YI2 

(QO) * ( AREF) 


(Eq. E-38) 


CYI = CYI1 + CYI2 


(Eq. E-39) 


Q. Flow-Throu gh Pressure Ratio 

1. The nozzle exit (flow-through) total pressure in ratio to free- stream static 
pressure is called PTD1/PO and PTD2/PO for ducts 1 and 2, respectively. 

2. PTO and PO are from tunnel parameters, module A. 

If IRAKE = 2, skip the remainder of this section. 

^/PO.IPTDI^MPTO) 

(Eq. EMO) 


If NPTR2 = 0, skip equation E-41. 


^ 2 / poJpto™wpto) 


(Eq. E-41) 
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|JF *p 

i 

E. Internal Drap ' " 

1. The internal drag coefficient based on the stability axes is called CDIS1 
and CDIS2 for ducts 1 and 2, respectively. 

2. The total internal drag coefficient in the stability axis is called CDIS. 

3. ALPHA is from the balance and weight tares computations, module D. 

4. The internal drag coefficient based on the wind axis is called CDIi and 
CDI2 for ducts 1 and 2, respectively. 

5. The total internal drag coefficient in the wind axis is called CDI. 


6. BETA is from module D. 

If IRAKE = 2, skip the remainder of this section. 

CDIS1 = (CNI1) * SIN(ALPHA) + (CAI1) * COS(ALPHA) (Eq. E-42) 

CDII = (CDIS1) * COS(BETA) - (CYI1) * SIN(BETA) (Eq. E-43) 

CDIS = CDIS1 (Eq. E-44) 

CDI = CDII (Eq. E-45) 

If NPTR2 = 0, skip equations E-46, E-47, E-48 and E-49. 

CDIS2 = (CNI2) * SIN(ALPHA) + (CAI2) * COS(ALPHA) (Eq. E-46) 

CDI2 = (CDIS2) * COS(BETA) - (CYI2) * SIN(BETA) (Eq. E-47) 

CDIS = CDIS1 + CDIS2 (Eq. E-48) 

CDI = CDII + CDI2 (Eq. E-49) 


E> 


171 


1. The internal lift coefficient based on stability axis (also wind axis) is called 
CLI1 and CLI2 for ducts 1 and 2, respectively. 


2. The total internal lift coefficient is called CLI. 


3. ALPHA is from the balance and weight tares computations in module D. 


If IRAKE = 2, skip the remainder of this section. 


CLI1 = (CNI1) * CQSK ALPHA) - (CAI1) * SIN(ALPHA) 


(Eq.E-50) 


CLI = CLI1 


(Eq.E-51) 


If NPTR2 = 0, skip equations E-52 and E-53. 


CLI2 = (CNI2) * COS( ALPHA) - (CAI2) * SIN(ALPHA) 


CLI = CLI). + CLI2 


If IRAKE * 3, skip this section. 


CAI, CDSI, and CDI are supplied in tables as functions of MACH, ALPHA 
and PSI. 


If IRAKE is 4, the following formulas are used. When the Mach number, 
MACH, is less than 0.1, then 
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CAI1 = CNI1 = CAI = CNI = CDIS = CD1 = CLI = CYI 1 = CY1 = 

MPM01 = MEXIT = MODOTI = FTMDOT1 = 0.0 

(Eq.E-54) 

For the Mach number, MACH, greater than 0.1, the following formulas are 
used. 

The area, AEXIT1, is used to interpolate R vs. Z to obtain YM. 

Then the average pressure, PL, is the average of the static pressures 

^Vrakeu) 

PL= 1 1 NPR1 

(Eq.E-55) 


When PL is greater than PTO, then go to equation (54). 
The Mach number, ML, is computed as 



(Eq. E-56) 


The ratio, BRATIO, becomes 


BRATIO s* 


L728 ML 


1+0.2ML 


(Eq. E-57) 
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Then the area, ASTAR, is 


ASTAR = n * YM 2 * B RATIO 


The Mach number ratio, MFR, is 


MFR = 


ML 


MACH 


(f 2^ 3 

1+0.2 MACH 

. i 

2 

1+0.2 ML 


YMf 

R(l)J 


The drag coefficients are now computed. 


The ratio value, BRATIO, is computed 

BRATIQ = ASTA ft 

riRd ) 2 


The value of the Mach number, M, is computed which satisfies 


1.728M/U + 0.2M 2 ) 3 - BRATIO = 0 


The static pressure, PL, is then computed 

pl= — pto — 
(1 + 0.2 M 2 )™ 


(Eq.E-58) 


(Eq. E-59) 


(Eq.E-60) 


(Eq. E-61) 


(Eq. E-62) 
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The pressure coefficient, CP(I), is 


CP(I) = (PL - PO)/QO 


The dynamic pressure, Q(I), is 


Q(I) = 0.7 * PL * M 2 


The slope, SLOPE, is 

SLOPE = TAN* 1 ((R(I) * R(I - 1))/(Z(I) - Z(I - 1))) 

The value of SL becomes 

SL = n(R'I) + R(I - 1)>/(R(I) - Rd - 1)) 2 + (Z<I> “ Z(I - 1)) 2 (M) 


The value of SWET is 


SWET = 


NPRI 

1 

i=l 


SL COS (SLOPE) / AREF 


The value of QRATIO is 


NPRI 


QRATIO = 1 

I=i 


( Q(I) + Q(I - 1)) SL » COS(SLOPE) 
2QS 


(Eq.E-63) 


(Eq.E-64) 


(Eq. E-65) 


(Eq.E-66) 


(Eq. E-67) 


(Eq.E-68) 
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The value of CDPA is 


NPRi n 

CDPA = - I (R(I) CP(I) + R(I - 1) CP(I - 1)) (R(I) - R(I - D) 

1=1 AKitir 

(Eq. E-69) 


The Mach number, MEFF, is computed to satisfy 


( MEFF Y 

f 2 1 

(1+0.2 MACrt ) 

ImachJ 

^ (1+0.2 MEFF 2 ) j 


- QRATIO = 0 


The static temperature, TE, is 


(Eq. E-70) 


TTO + 459.67 
1+0.2 MEFF 2 


(Eq. E-71) 


The density, RHOE, is 


RHOE = 


144.0*PQ 
17 16. 4829 TE 


(Eq. E-72) 


The viscosity, VISE, is 


_ 2.27*10 8 »TE*VTE 
TE + 198.6 


(Eq. E-73) 


The free-stream velocity is 


UE = 49.021179*MEFF*VTE 


The local Reynolds number is 


RNA = 


RHOE * UE * Z(NPRT1) 
12 * VISE 


When the Reynolds number is positive and nonzero 

0.472 


CFO = 


(LOG 10 (RNA)) 


2.58 


( 


1 + 0.2 MEFF 


2 \ 0.467 


V 


(Eq. E-74) 


(Eq. E-75) 


(Eq. E-76) 


Then 


COF = CFO * QRATIO * SWET 


(Eq.E-77) 


The pre-entry drag, CDPRE, is 


CDPRE = - 


f 

f 

2MFR 

M'l) 

MACH1 


» 


1+0.2 MACH 
1 + 0.2 M 2 


-COS( ALPHA) 


+ CP(I) 


HR(1) 

AREF 


(Eq. E-78) 


177 


CAI1 = CDPRE + (CDPA + COF) COS (PSIN) COS (THETAN) 



The internal normal force coefficient is 

2 

CAI1 = 2 * MFR * — — r SIN( ALPHA) + (CDPA + COF) SIN(THETAN) 
AREF 


The total axial force coefficient is 

CAI = NPRT2 * CAI1 


The total normal force coefficient is 

CNI = NPRT2 * CNI1 


The total internal drag force coefficient in the stability axis is 
CDIS = CNI * SIN(ALPHA) + CAI * COS(ALPHA) 



The total internal drag force coefficient in the wind axis is 

CDI = CDIS * COS(BETA) 


The total internal lift coefficient is 

CLI = CNI * COSCALPHA - CAI * SIN(ALPHA) 
The total internal side force coefficient is 


(Eq. E-79) 


(Eq.E-80) 


(Eq. E-81) 


(Eq. E-82) 


(Eq. E-83) 


(Eq. E-84) 


(Eq. E-85) 
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The total internal side force coefficient is 


CYI1 = CYI = 0.0 


The mass flow ratio is 


M/M01 = MFR 


The average exit Mach number is 

MEX1T1 = ML(NPRTl) 


The mass flow rate based on free-stream conditions is 


MODOT1 = 0.028563 * SCAP1 * MACH * PO * 


The mass flow ratio at exit is 

FTMDOT1 = MFR * MODOT1 


(Eq. E-86) 


(Eq. E-87) 


(Eq. E-88) 


2 

1 + 0.2 MACH 

(Eq. E-89) 


(Eq. E-90) 
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MODULE P PRESSURE COEFFICIENT AND INTEGRATED FORCES 


SYMBOL 


NOMENCLATURE 


ARAAUCIfi 

ARABU(I) 

ARADU(I) 

ARAFU(I) 

ARAGU(I) 

ARAHU(I) 

ARANU(I) 

ARASU(I)^ 


Axial 

Force 


Areas to be used with pressure groups to compute 
integrated forces, where I = orifice number, sq. in. 


AREF 


ARHAU(I) 

ARHBU(I) 

ARHDU(I) 

ARHFU(I) 

ARHGU(I) 

ARHHU(I) 

ARHNU(I) 

ARHSU(I) 


A 




J 


Hinge 

Moment 


Model reference area from module B, sq. in. 


Areas times moment arm to be used with pressure 
groups to compute integrated forces, where 
I = orifice number, sq. in. 


ARNAUdfi 

ARNBU(I) 

ARNDUd) 

ARNFU(I) 

ARNGUCI) 

ARNHUd) 

ARNNU(I) 

ARNSUd). 


Normal 

Force 


Areas to be used with pressure groups to compute 
integrated forces, where I = orifice number, sq. in. 


ARPAUdfi 

ARPBU(I) 

ARPDU(I) 

ARPFU(I) 

ARPGUd) 

ARPHU(I) 

ARPNU(I) 

ARPSUd) J 


Pitch 

Moment 


Area times moment arm to be used with pressure 
groups to compute integrated moments, I = orifice 
number, sq. in. 
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SYMBOL NOMENCLATURE 


CBAR 


Pitching moment reference length module D, in. 


CDAUN(in 

CDBUN(I) 

CDDUN(I) 

CDFUN(I) L 

CDGUN(I) 

CDHUN(I) 

CDNUN(I) 

CDSUN(I)J 


Integrated pressure drag coefficients. 


CDPR 


Total integrated drag coefficient. 


CFAUNdfi 

CFBUN(I) 

CFDUN(I) 

CFFUN(I) 

CFGUN(I) 

CFHUN(I) 

CFNUN(I) 

CFSUN(I) 


Integrated pressure axial force coefficients. 


CHMAUN(I) A 

CHMBUN(I) 

CHMDUN(I) 

CHMFUN(I) 

CHMGUN(I) 

CHMHUN(I) 

CHMNUN(I) 

CHMSUN(I)J 


Integrated pressure hinge moment coefficients. 


* 


CLAUN(If| 

CLBUN(I) 

CLDUN(I) 

CLFUN(I) 

CLGUN(I) 

CLHUN(I) 

CLNUN(I) 

CLSUN(I) 


Integrated pressure lift coefficients. 
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CLPR 


Total integrated lift coefficient. 


CNAUNdfl 
CNBUN(I) 
CNDUN(I) 
CNPUN(I) L 
CNGUN(I) 
CNHUN(I) 
CNNUN(I) 
CNSUN(I) J 


CPMAUNdfi 

CPMBUNd) 

CPMDUNd) 

CPMFUNd) 

CPMGUNd) 

CPMHUNd) 

CPMNUNd) 

CPMSUN(I)^ 


CPMPR 

KCDA 

KCDB 

KCDD 

KCDF 

KCDG 

KCDH 

KCDN 

KCDS 


PAUN(I)') 
PBUNd) 
PDUNd) 
PFUNd) 
PGUNd) f 
PHUNd) ! 
PNUNd) 
PSUN(I)J 


Integrated pressure normal force coefficients. 


Integrated pressure pitching moment coefficients. 


Total integrated pitching 
Constants provided by the 
Constants provided by the 
Constants provided by the 
Constants provided by the 
Constants provided by the 
Constants provided by the 
Constants provided by the 
Constants provided by the 


moment coefficient, 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 
engineer. (See note) 


Individual pressures to be used with each type of 
pressure coefficient for computation of integrated 
forces and moments. Maximum number of each 
type is 125. 
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PRATI(II) 


TCI?! 


# 


Ratio of nozzle internal static pressure to nozzle total 
pressure, where II = orifice number (PGUN only). 


Enter 1.0 if term is to be included in the total pressure drag 
coefficient, or 0.0 if it is to be excluded. 
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APPENDIX F 
Module F 

Pressur e Coefficients and Integrated Forces 

Eight groups of pressure coefficients may be computed under this module. 
Names assigned to each group are arbitrary. Final names may be inserted with 
finalized data printout headers. These groups are PAUN, PBUN, PDUN, PFUN, 
PGUN, PHUN, PNUN, and PSUN. 

A. Required Constants 

The required constants for module F are given in the nomenclatures. 

1. All constants are initialized to a value of zero. The project engineer 
need only supply those constants which are required for those 
quantities to be computed. 

2. KAUN, KBUN, KDUN, KFUN, KGUN, KHUN, KNTJN, KSUN - 
number of individual pressures to be used with each type of pressure 
coefficient for computation of integrated forces and moments. (125 
maximum for each type) 

B. Test for Module F Computations 
If KPRESS = 0, skip this module. 

If KPRESS = 1 performs calculations as listed in module along with ratio of 
engine static pressure to engine total pressure. 

If KPRESS = 2 performs calculations as listed in module along with ratio of 
engine static pressure to free stream total pressure. 

C. Frec^Stream Static and Dvnamip P ressures 

Free-stream static and dynamic pressures to be used for computing pressure 
coefficients are obtained from module A; however, for individual pressure 
transducers, an average value is used. 


D. Coefficient Calculations 


CPSUN(I) = (PSUN(I) - PO)/QO 


KSUN 

CFSUN = I CPSUN(I)*ARASU(I)/ AREF 
1=1 


(Eq. F-l) 


(Eq. F-2) 


KSUN 

CNSUN = I CPSUN(I) * ARASU(I) / AREF 
1=1 


KSUN 

CPMSUN = X CPSUN(I)* ARPSU(I) / AREF*CBAR 
1=1 


(Eq. F-3) 


(Eq.F-4) 


KSUN 

CHMSUN = X CPSUN(I)* ARHSU(I) / AREF*CBAR 

l=i 


(Eq. F-5) 


CDSUN = CFSUN * COS (ALPHA) + CNSUN * SIN (ALPHA) 


(Eq. F-6) 


CLSUN = CNSUN * COS (ALPHA) - CFSUN * SIN (ALPHA) 


(Eq. F-7) 


These equations are the same for all pressure groups. 
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E. Total Pressure Drag Coefficient 


CDPRessure: = KCDS * CDSUN + KCDA * CDAUN + KCDB * CDBUN + KCDN * 
CDNUN + KCDD * CDDUN + KCDF * CDFUN + KCDG * CDGUN + KCDH * 
CDHUN 

(Eq. F-8) 

where KCDS, KCDB, KCDA, KCDN, KCDD, KCDF, KCDG, KCDH are 
constant inputs, 1.0 if term is to be included in the total pressure drag 
coefficient, or 0.0 if it is to be excluded. 

F. Internal Static Press ure Ratio 
If KPRESS = 1 

PRATIdl) = PGUN(n)/PTENGl (Eq. F-9) 

If KPRESS = 2 

PRATKII) = PGUN (II)/PTO 

(Eq. F-10) 

NOTE: In addition to the pressure coefficients, this ratio is for PGUN 

measurements only. 
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APPENDIX G 

Thrust Removal Options 


Nomenclature 193 

General Information 198 

Required Constants 198 

Quantities Required 198 

Compute Thrust and Static Thrust Terms 

IF = 1 199 

Single Balance/All Metric IF1 = 1 206 

Single Balance/Afterbody Metric IF2 = 1 209 

Two-Balance/Afterbody Metric IFAF1 = 1 212 

Two-Balance/Afterbody Metric IFAF2 =1 215 

Two-Balance/ Afterbody Metric IFAFN1 = 1 218 

Two-Balance/Afterbody Metric IFAFN2 = 1 221 

Single Balance, Thrust Removal All 

Components IFAF = 1 224 

When IFAFN = 1.. 229 

Bifurcate Support Mode Two Balance/Afterbody 

Metric IDN = 1 231 

Other Options 236 


w 

i 

! MODULE G THRUST REMOVAL OPTIONS 


SYMBOL 

NOMENCLATURE 

AEX 

Nozzle exit area. 

CAAERO 

Thrust removed axial force coefficient. 

CASCADE 

Resultant angle of jet exhaust, degrees. 

CDAERO 

Thrust removed drag coefficient. 

CDNOZ 

Nozzle drag. 

C(F-ANOZ) 

Thrust minus nozzle axial force coefficient. 

C(F-DNOZ) 

Thrust minus nozzle drag coefficient. 

CDSAER 

Thrust removed stability axis drag coefficient. 

CDWAER 

Thrust removed wind axis drag coefficient. 

CF 

Jet axial force coefficient (from balance and pressures). 

CF/CFI 

Ratio of thrust (from balance and pressures) to ideal thrust. 

CFJ 

Jet axial force coefficient. 

CFJC 

Computed jet axial force coefficient. 

CFJET 

Jet reaction axial force coefficient. 

C(F-A) 

Thrust minus axial force coefficient. 

C(F-D) 

Thrust minus drag coefficient. 

CLAERO 

Thrust removed lift coefficient. 

CIJET 

Jet reaction lift coefficient. 

CLNOZ 

Thrust removed nozzle lift. 

CLNQZT 

Nozzle lift plus thrust. 

CLSAER 

Thrust removed stability axis lift coefficient. 

CLWAER 

Thrust removed wind axis lift coefficient. 

CMAERO 

Thrust removed pitching moment coefficient. 

CMJ 

Jet pitching moment coefficient. 







193 


SYMBOL 




CMJC 

CMJET 

CMNOZ 

CMNOZT 

CMSAER 

CMWAER 

CNAERO 

CNJ 

CNJC 

CRAERO 

CRJC 

CRJET 

CRSAER 

CRWAER 

CSAERO 

CSJC 

CSJET 

CSSAER 

CSWAER 

CT 

CTS 

CTST 

CTSY 


Computed jet pitching moment coefficient. 

Jet reaction pitching moment coefficient. 

Thrust removed nozzle pitching moment. 

Nozzle pitching moment plus lift. 

Thrust removed stability axis pitching moment coefficient. 

Thrust removed wind axis pitching moment coefficient. 

Thrust removed normal force coefficient. 

Jet normal force coefficient. 

Computed jet normal force coefficient. 

Thrust removed rolling moment coefficient. 

Computed jet rolling moment coefficient. 

Jet reaction rolling moment coefficient. 

Thrust removed stability axis rolling moment coefficient. 

Thrust removed wind axis rolling moment coefficient. 

Thrust removed side force coefficient. 

Computed jet side force coefficient. 

Jet reaction side force coefficient. 

Thrust removed stability axis side force coefficient. 

Thrust removed wind axis side force coefficient. 

Computed resultant thrust coefficient about pitch axis. 

Resultant static thrust coefficient, main balance, about pitch 
axis. 

Resultant static thrust coefficient, main balance. 

Resultant static thrust coefficient, main balance, about yaw 
axis. 


i 




; 


J! 




H 
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SYMBOL 




CTS2 


CTS2T 

CTS2Y 


CTT 

CTY 

CYAERO 

CYJC 

CYJET 

CYSAER 

CYWAER 

DELTA 

DELTAY 

DELTA1 


DELTA2 


DELT1Y 


DELT2Y 


ETAABS 

(F-A)/FI 

(F-ANOZ)/FI 

FGT/FI 

FGT2/FI 

FGY/FI 



Resultant static thrust coefficient, second balance, about pitch 
axis. 

Resultant static thrust coefficient, second balance. 

Resultant static thrust coefficient, second balance, about yaw 
axis. 

Computed resultant thrust coefficient. 

Computed resultant thrust coefficient about yaw axis. 

Thrust removed yawing moment coefficient. 

Computed jet yawing moment coefficient. 

Jet reaction yawing moment coefficient. 

Thrust removed stability axis yawing moment coefficient. 

Thrust removed wind axis yawing moment coefficient. 

Computed thrust vector angle about pitch axis, degrees. 

Computed thrust vector angle about yaw axis, degrees. i 

| 

Static thrust vector angle, main balance, about pitch axis, ' 

degrees. 

Static thrust vector angle, second balance, about pitch axis, ■ 

degrees. ! 

i 

Static thrust vector angle, main balance, about yaw axis, 

degrees. ( 

Static thrust vector angle, second balance, about yaw axis, 
degrees. 

Isentropic vacuum or stream thrust coefficient. 

Ratio of thrust minus axial force to ideal thrust. 

Ratio of thrust minus nozzle axial force to ideal thrust. 

Total static resultant thrust ratio, main balance. 

Total static resultant thrust ratio, second balance. 

Static resultant thrust ratio, main balance, about yaw axis. 

195 

f>i 
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SYMBOL 

FG/FI 

FG2/FI 

FG2Y/FI 

FJ1/FI 

FJ2/FI 

(F-D)/FI 

(F-DNOZVFI 

FN/FI 

FT/FI 

F/FI 

IDA 

IDN 

IF 

IF1 

IF2 

IFAF 

IFAF1 

IFAF2 

IFAFN 

IFAFN1 

IFAFN2 

LENGTH(I) 

PM/FI 

RM/FI 

SF/FI 


NOMENCLATURE 

Static resultant thrust ratio, main balance, about pitch axis. 
Static resultant thrust ratio, second balance, about pitch axis. 
Static resultant thrust ratio, second balance, about yaw axis. 
Static thrust ratio, main balance. 

Static thrust ratio, second balance. 

Ratio of thrust minus drag to ideal thrust. 

Ratio of thrust minus nozzle drag to ideal thrust. 

Ratio of normal force to ideal thrust. 

Total resultant thrust ratio. 

Ratio of thrust to ideal thrust. 

Engineer's option. 

Future option. 

Computes thrust and static thrust terms when IF=1. 
Computes single balance/all metric when IF1=1. 

Computes single balance/afterbody metric when IF2=1. 
Single balance, thrust removal from all components. 
Computes two balances/afterbody metric when IFAF1=1. 
Computes two balances/afterbody metric when IFAF2=1. 
Future option. 

Computes two balances/afterbody metric when IFAFN1=1. 
Computes two balances/afterbody metric when IFAFN2=1. 
Lengths for transferring moments to relative station. 

Ratio of pitching moment to ideal thrust. 

Ratio of rolling moment to ideal thrust. 

Ratio of side force to ideal thrust. 
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SYMBOL 


NOMENCLATURE 


SPLAY 

SPLAY1 

YM/FI 


Projected roll angle of jet exhaust, degrees. 
Projected roll angle of jet exhaust, degrees. 
Ratio of yawing moment to ideal thrust. 


APPENDIX G 
Module G 

Thrust Removal Options 


r 


i A. General Information 

The following options are used to remove thrust and to obtain various 
. aerodynamic and aeropropulsion parameters usually required for most 16-Ft. 

i Transonic Tunnel investigations. The various constants are keyed to typical 

balance arrangements used and may be used for most test setups. This 
‘ section requires computed inputs from modules A, B, C, D and F. The 

t en P neer should refer to each module for exact definition of the computed 

quantity. These options will work for both fully and partially metric models for 
j both longitudinal and lateral data. 

I 

f 

B. Required Constants 

1. IF, IF1, IF2, IFAF, IFAF1, IFAF2, IFAFN, IFAFN1, IFAFN2, ID and 
IDN. 

C. Quantities Required 

MODULE A 

a. PO&QO 

MODULE B 

a. NPR 

b. CFI 

MODULE C 

a. CDFAFT * afterbody + nozzle skin friction 

b. CDFNOZ * nozzle skin friction 
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4. MODULE D 


1. ALPHA 

2. CN1, CA1, CMY1 

CY1, CMX1, CMZ1 MAIN BALANCE 

3. CDS1.CLS1 

4. CN2, CA2, CMY2 

CY2, CMX2, CMZ2 SECOND BALANCE (2) 

5. CDS2.CLS2 

5. MODULE F 

1. CFSUN, CNSUN, CPMSUN 

2. CDSUN, CLSUN AFTERBODY PRESSURE FORCES 

3. CFBUN, CNBUN, CPMBUN 

4. CDBUN, CLBUN NOZZLE PRESSURE FORCES 

D. Compute T hrust and Static Thrust Terms IF = 1 
1. Compute Thrust 

The value of NPR is the average nozzle pressure ratio from each air system, 

a. IF NPR <; 1.2, CFJ = CMJ = 0 = CRJ = CYJ = CSJ 


b. The computed jet axial force coefficient is 

CFJC = ^ * (KCFJ(NPR) + ICFJ] 


c. The computed jet normal force coefficient is 


(Eq. G-l) 
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(Eq. G-2) 


CNJC = ^ *f R CNJ(NPR) + ^CNJ] 


d. The computed jet pitching moment coefficient is 

PO 

CMJC = * [KCMJ(NPR) + ICMJ] 


e. The computed jet rolling moment coefficient is 

PO 

CRJC = Sr [KCRJ(NPR) + ICRJl 

lq|U 


f. The computed jet yawing moment coefficient is 

PO 

CYJC = ♦ [KCYJ(NPR) + IC YJ] 


g. The computed jet side force coefficient is 

PO 

CSJC = ^ * [ RCSJ(NPR > + ICSJ ] 


(Eq. G-3) 


(Eq. G-4) 


(Eq. G-5) 


(Eq. G-6) 


h. Table input is as follows: (Need six tables) Up to five 
values per table 
may be used. 

NPR Range Slope 
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Intercept 


2. Compute Static Thrust Terms 

a. The resultant static thrust coefficient about the pitch axis for the main 
balance is 


CTS = 


i 


2 2 
CNl +CA1 


(Eq.G-7) 


b. The resultant static thrust ratio about the pitch axis for the main balance 
is, where CPI is the average of each air system 


FG/FI = CTS/CFI 


(Eq. G-8) 


c. The static thrust ratio for the main balance is 

FJ1/FI = -CA1/CFI 


(Eq. G-9) 


d. The static thrust vector angle about the pitch axis for the main balance is 


DELTA1 = TAN* 1 (-CN1/CA1) 


(Eq. G-10) 


e. The resultant static thrust coefficient about the yaw axis for the main 
balance is 


I 2 

CTSY = VCYl +CA1 


(Eq. G-ll) 
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f. The resultant static thrust ratio about the yaw axis for the main balance is 


FGY/FI = CTSY/CFI 

(Eq. G-12) 


g. The static thrust vector angle about the yaw axis for the main balance is 


DELTA1Y = TAN 1 (-CY1/CA1) 

(Eq.G-13) 


h. The resultant static thrust coefficient for the main balance is (Eq. G-14) 


222 
CTST = vCNl + CA1 + CY1 


i. The total resultant static thrust ratio for the main balance is 


FGT/FI = CTST/CFI 

(Eq. G-15) 


j _ The isentropic vacuum thrust or stream thrust coefficient. 

Specific equations to be supplied when required by test hardware. 


k. The nozzle exit Mach number ME is computed from 


AS/ AEX = 


216 

125 


ME(1 + 0.2ME 2 )' 3 


(Eq. G-16) 


where AS = WPWITO(M) * AT(1). 
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and the nozzle exit pressure ratio is from 

PE/PTJ = (1 + 0.2 ME 2 r 7/2 


1. The resultant static thrust coefficient about the pitch axis for the second 
balance is 


CTS2 = vCN2 2 + CA2 2 


m . The resultant static thrust ratio about the pitch axis for the second 


balance is 


FG2/FI = CTS2/CFI 


n . The static thrust ratio for the second balance is 


FJ2/FI = -CA2/CFI 


(Eq.G-20) 


o. The static thrust vector angle about the pitch axis for the second balance is 


DELTA2 = TAN 1 (-CN2/CA2) 


(Eq. G-21) 


p. The resultant static thrust coefficient about the yaw axis for the second 
balance is 


CTS2Y = VCY2 2 +CA2 2 


(Eq. G-22) 
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q. The resultant static thrust ratio about the yaw axis for the second balance 
is 


FG2Y/FI = CTS2Y/CFI 


(Eq. G-23) 

r. The static thrust vector angle about the yaw axis for the second balance is 

DELTA2Y = TAN' 1 (-CY2/CA2) 


(Eq. Gr- 24 ) 


s. The resultant static thrust coefficient for the second balance is 


CTS2T 


J 


CY2 2 + CA2 2 + CY2 2 


(Eq. G-25) 


t. The total resultant static thrust ratio for the second balance is 


FGT2 = CTS2T/CFI 


(Eq. G-26) 


u . The splay angle is 


SPLAY = ATAN(CY1/CN1) 


(Eq. G-27) 


v. The cascade angle is 


CASCADE = TAN 


-l 


£ 


TAN (DELTA1) 


TAN 2 (SPLAY)*TAN 2 (DELTA1) 


(Eq. G-28) 
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w . The ratio of normal force to ideal thr ust is 


FN/FI = CN1/CFI 


(Eq. G-29) 


x . The ratio of side force to ideal thrust is 

SF/FI = CY1/CFI 


(Eq. G-30) 


y. The ratio of rolling moment to ideal thrust is 

RM/FI = CMX1/(CFI * LENGTH1) 


(Eq. G-31) 


z. The ratio of pitching moment to ideal thrust is 

PM/FI = CMY1/(CFI * LENGTH2) 


(Eq. G-32) 


aa. The ratio of yawing moment to ideal thrust is 

YM/FI = CMZ1/(CFI ♦ LENGTH3) 


(Eq.G-33) 



1. The resultant thrust coefficient about the pitch axis is 


ct = \cnjc 2 +cfjc 2 




(Eq.G-34) 


2. The thrust vector about the pitch axis is 

DELTA = TAN 1 (CNJC/CFJC) 


3. The jet reaction lift coefficient is 

CLJET = CT [SIN (ALPHA + DELTA)] 


(Eq. G-35) 


(Eq. G-36) 
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4. The jet reaction axial force coefficient is 

CFJET « CT [COS (ALPHA + DELTA)] 

5. The thrust removed lift coefficient is 

CLAERQ = CLS1 - CUET 

6. The thrust removed drag coefficient is 

CDAERO = CDS1 + CFJET 

7. The thrust removed pitching moment coefficient is 

CMAEEO = OMYS1 - CMJC 

8. The thrust minus axial force coefficient is 

C(F - A) = -CA1 

9. The thrust minus drag coefficient is 

C(F - D) = -CDS1 

10. The ratio of thrust minus axial force to ideal thrust is 

(F - AVFI * C(F - A)/CFI 



(Eq. G-37) 


(Eq. G-38) 


(Eq. G-39) 


(Eq.G-40) 


(Eq. G-41) 


(Eq. G-42) 


(Eq. G-43) 
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11. The ratio of thrust minus drag to ideal thrust is 


(F - DyFI * C(F - D)/CFI 

(Eq.G-44) 

12. The ratio of thrust to ideal thrust is 

F/FI = CFJC/CF1 

(Eq. G-45) 

13. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) - CFJET - CDBUN 

(Eq. G-46) 


14. The ratio of thrust minus nozzle drag to ideal thrust is 

(F - DNOZVFI = C(F - DNOZ)/CFI 


(Eq. G-47) 
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1. The resultant thrust coefficient about the pitch axis is 


ct = )/cnjc 2 +cfjc 2 


(Same as Eq. G-34) 


2. The thrust vector angle about the pitch axis is 


DELTA = TAN 1 (CNJC/CFJC) 


(Same as Eq. G-35) 


3. The jet reaction lift coefficient is 


CLJET = CT [SIN (ALPHA + DELTA)] 


(Same as Eq. G-36) 
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4. The jet reaction axial force coefficient 


CFJET = CT [COS (ALPHA + DELTA)] 


(Same as Eq. G-37) 


5. The thrust removed lift coefficient is 


CLAERO = CLS1 - CLJET 


(Same as Eq. G-38) 


6. The thrust removed drag coefficient 


CDAERO = CDS1 + CFJET 


7. The thrust removed pitching moment coefficient 

CMAERO = CMYS1 - CMJC 


8. The thrust minus axial force coefficient 


C(F - A) = -CA1 


9. The thrust minus drag coefficient 


(Same as Eq. G-39) 


(Same as Eq. G-40) 


(Same as Eq. G-41) 


C(F - D) = -CDSl 


(Same as Eq. G-42) 
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10. The ratio of thrust minus axial force to ideal thrust is 


(F - AVFI « C(F - AyCFI 

(Same as Eq. G-43) 


11. The ratio of thrust minus drag to ideal thrust is 

(F - DyFI = C(F - DVCFI 


(Same as Eq.G-44) 


12. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) = C(F - D) + (CDFAFT - CDFNOZ) 

(Eq.G-48) 


13. The ratio of thrust minus nozzle drag to ideal thrust is 

(F - DNOZyFI = C(F - DNOZVCFI 


(Eq.G-49) 


14. The ratio of thrust to ideal thrust is 

F/FI = [C(F - DNOZ) + CDFNOZ + CDBUN]/CFI 

(Eq. G-60) 


1 • 
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^ 777777 / 77777 ? 


Balance l-(D-T) 

Balance 2 -(DAFT + DNOZ) 


1. The jet axial force coefficient is 

CFJ = CA2 - CA1 


(Eq. G-51) 


2. The jet normal force coefficient is 

CNJ * CN1 - CN2 


(Eq. G-52) 


3. The jet pitching moment coefficient is 

CMJ = CMY1 - CMY2 


(Eq. G-53) 


4. The resultant thrust coefficient about the pitch axis is 



CT = 




cfj 2 +cnj 2 


5. The thrust vector angle about the pitch axis is 

DELTA = TAN 1 (CNJ/CFJ) 


6. The thrust removed lift coefficient is 

CLAERO = CLS2 


7. The thrust removed drag coefficient is 

CDAERO = CDAS2 


(Eq.G-54) 


(Eq. G-55) 



if 


i 

f 1 


i 


(Eq. G-56) 


(Eq. G-57) 


8. The thrust removed pitching moment coefficient is 

CMAERO = CYMS2 


I 

i 


(Eq. Gr-58) 


9. The thrust minus axial force coefficient is 


^ ♦ 


C(F - A) = -CA1 


(Same as Eq. G-41) 
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10. The thrust minus drag coefficient is 


C(F - D) = -CDSl 


11. The ratio of thrust minus axial force to ideal thrust is 


(F - AVFI = C(F - AVCFI 


12. The ratio of thrust minus drag to ideal thrust is 


(F - D)/FI = C(F - D)/CFI 


13. The thrust minus nozzle drag coefficient is 


(Same as Eq. G-42) 


(Same as Eq. G-43) 


(Same as Eq. G-44) 


C(F - DNOZ) = C(F - D) + (CDFAFT - CDFNOZ) + CDSUN 


14. The ratio of thrust minus nozzle drag to ideal thrust is 


(Eq. G-59) 


(F - DNOZVFI = C(F - DNOZVCFI 


(Same as Eq. G-49) 



16. The ratio of thrust to ideal thrust is 


F/FI = CFJ/CFI 


(Eq. G-60) 
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16. The jet reaction lift coefficient is 

CLJET = CT [SIN(ALPHA + DELTA)] 

17. The jet reaction axial force coefficient is 

CFJET = CT [COS( ALPHA + DELTA)] 



(Eq. G-61) 


(Eq. G-62) 



Note: Wings or tails usually attached 
to balanco 2 


1. The thrust removed lift coefficient is 

CLAERO = CLS2 


(Same as Eq. G-56) 
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2. The thrust removed drag coefficient is 


CD AERO = CDS2 


(Same as Eq. G-57) 


3. The thrust removed pitching moment coefficient is 

CMAERO = CYMS2 


(Same as Eq. G-58) 


4. The thrust minus axial force coefficient is 


C(F - A) = -CA1 


(Same as Eq. G-41) 


5. The thrust minus drag coefficient is 


C(F - D) = -CDS1 


(Same as Eq. G-42) 


6. The ratio of thrust minus axial force to ideal thrust is 

(F - A)/FI = C(F - A)/CFI 


(Same as Eq. G-43) 


7. The ratio of thrust minus drag to ideal thrust is 

(F - D)/FI = C(F - D)/CFI 


(Same as Eq. G-44) 
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8. The thrust minus nozzle drag coefficient is 


W 


L 



OF - DNOZ) . C(F - D) + CDS2 


(Eq.G-63) 


9. The jet axial force coefficient is 

CFJ = C(F - DNOZ) + (CDFNOZ + CDBUN) 


(Eq. G-64) 


10. The jet normal force coefficient is 


CNJ = CN1 - CN2 - CNBUN 


(Eq. G-65) 


11. The jet pitching moment coefficient is 

CMJ = CMY1 - CMY2 - CPMBUN 


(Eq. G-66) 


12. The resultant thrust coefficient about the pitch axis is 


CT 


-J 


CFJ 2 + CNJ 2 


(Same as Eq. G-54) 


! i: 


i 

-l 

4 


t 


4 




13. The thrust vector angle about the pitch axis is 


DELTA = TAN 1 (CNJ/CFJ) 


(Same as Eq. G-55) 
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14. The jet reaction lift coefficient is 


CLJET = CT [SIN(ALPHA + DELTA)] 

(Same as Eq. G-61) 


15. The jet reaction axial force coefficient is 

CFJET = CT [COS( ALPHA + DELTA)] 

(Same as Eq. G-62) 


I. Two-Balance/ Afterbo dy Metric IFAFN1 = 1 



1. The jet axial force coefficient is 


CFJ = -CA2 


(Eq. G-67) 


2. The jet normal force coefficient is 

CNJ = CN2 


(Eq. G-68) 
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3. The jet pitching moment coefficient is 


CMJ = CMY2 


(Eq.G-69) 


4. The resultant thrust coefficient about the pitch axis is 


CT 


4 


cfj 2 +cnj 2 


(Same as Eq. G-54) 


5. The thrust vector angle about the pitch axis is 

DELTA = TAN 1 (CNJ/CFJ) 


(Same as Eq. G-55) 


6. The jet reaction lift coefficient is 

CLJET = CT [SIN(ALPHA + DELTA)] 


(Same as Eq. G-61) 


7. The jet reaction axial force coefficient is 

CFJET = CT [COS( ALPHA + DELTA)] 

(Same as Eq. G-62) 


8. The thrust minus axial force coefficient is 


C(F - A) » -CA1 


(Same as Eq. G-41) 
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9. The thrust minus drag coefficient is 



i 

\ 




I 


C(F - D) * -CDS1 


(Same as Eq. G-42) 


10. The ratio of thrust minus axial force to ideal thrust is 


(F -AVFI = C(F - AVCFI 


(Same as Eq. G-43) 


11. The ratio of thrust minus drag to ideal thrust is 


(F - DYFI = C(F - D)/CFI 


(Same as Eq. G-44) 


12. The thrust removed lift coefficient is 

CLAERO = CLS1 - CLS2 


(Eq. G-70) 


13. The thrust removed drag coefficient is 

CDAERO = CDS1 - CDS2 


(Eq. G-71) 


14. The thrust removed pitching moment coefficient is 

CMAERO = CMYS1 - CMYS2 


(Eq. G-72) 
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15. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) a C(F -D) + (CDFAFT - CDFNOZ) + CDSUN 

(Same as Eq. G-59) 


lyo-Balance/ Afterbody Met ric IFAFN2 = 1 



1. The jet axial force coefficient is 

CFJ = (CDFNOZ + CFBUN) - CA2 


(Eq. G-73) 


2. The jet normal force coefficient is 


CNJ = CN2 - CNBUN 


(Eq. G-74) 
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3. The jet pitching moment coefficient is 


CMJ = CMY2 - CPMBUN 


(Eq. G-75) 


4. The resultant thrust coefficient about the pitch axis is 


ct=\cfj 2 +cnj 2 


(Same as Eq. G-54) 


5. The thrust vector about the pitch axis is 

DELTA = TAN 1 (CNJ/CFJ) 


(Same as Eq. G-55) 


6. The jet reaction lift coefficient is 

CLJET = CT * [SIN(ALPHA + DELTA)] 


(Same as Eq. G-61) 


7. The jet reaction axial force coefficient is 

CFJET = CT * [COS( ALPHA + DELTA)] 


(Same as Eq. G-62) 


8. The thrust minus axial force coefficient is 


C(F - A) = -CA1 


(Same as Eq. G-41) 
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9. The thrust minus drag coefficient is 


C(F - D) = -CDS1 


(Same as Eq. G-42) 


10. The ratio of thrust minus axial force to ideal thrust is 


(F - A)/FI = C(F - A)/CFI 


(Same as Eq. G-43) 


11. The ratio of thrust minus drag to ideal thrust is 


(F - D)/FI = C(F - D)/CFI 


(Same as Eq. G-44) 


12. The thrust removed lift coefficient is 

CLAERO = CLS1 - CLS2 


(Same as Eq. G-70) 


13. The thrust removed drag coefficient is 

CDAERO = CDS1 - CDS2 


(Same as Eq. G-71) 


14. The thrust removed pitching moment coefficient is 

CMAERO = CMYS1 - CMYS2 

(Same as Eq. G-72) 
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(Eq. G-76) 

l 



1. The thrust removed normal force coefficient is 


CNAERO = CN1 - CNJC 


(Eq. G-77) 



2. The thrust removed axial force coefficient is 

CAAERO = CA1 + CFJC 

3. The thrust removed pitching moment coefficient is 

CMAERO = CMY1 - CMJC 

4. The thrust removed rolling moment coefficient is 

CRAERO = CMX1 - CRJC 

5. The thrust removed yawing moment coefficient is 

CYAERO = CMZ1 - CYJC 


(Eq. G-78) 


(Eq. G-79) 


(Eq. G-80) 


(Eq. G-81) 
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6. The thrust removed side force coefficient 


is 


CSAERO = CY1 - CSJC 


(Eq. G-82) 


7. The thrust removed lift coefficient is 


CLAERO = CNAERO * COS(ALPHA) - CAAERO ♦ SIN(ALPHA) 

(Eq. G-83) 

8. The thrust removed drag coefficient is 


CDAERO = CAAERO * COS(ALPHA) + CNAERO * SIN(ALPHA) 

(Eq. G-84) 


9. The resultant thrust coefficient about the pitch axis 

CT = a/cNJC 2 + CFJC 2 


18 


(Same as Eq. G-34) 


10. The thrust vector angle about the pitch axis is 


DELTA = TAN 1 (CNJC/CFJC) 


(Same as Eq. G-35) 


11. The jet reaction lift coefficient is 


CLJET = CNJC * COS( ALPHA) + CFJC * SIN(ALPHA) 


(Eq. G-85) 
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12. The jet reaction axial force coefficient is 


CFJET = CFJC * COS( ALPHA) - CNJC * (ALPHA) 


(Eq. G-86) 


13. The jet reaction side force coefficient is 


CSJET = CSJC 


(Eq. G-87) 


14. The jet reaction pitching moment coefficient 


is 


CMJET = CMJC 


(Eq. G-88) 


15. The jet reaction rolling moment coefficient is 


CRJET * CRJC * COS(ALPHA) + CYJC * SIN(ALPHA) 


(Eq. G-89) 


16. The jet reaction yawing moment coefficient 


is 


CYJET = CYJC * COS( ALPHA) - CRJC * SIN(ALPHA) 


(Eq. G-90) 


17. The splay angle is 


SPLAY1 = TAN 1 (CSJC/CNJC) 


(Eq.G-91) 
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18. The thrust minus axial force coefficient is 


C(F - A) = -CA1 


(Same as Eq. G-41) 


19. The thrust minus drag coefficient is 


C(F - D) = -CDS1 


(Same as Eq. G-42) 


20. The ratio of thrust minus axial force to ideal thrust is 

(F - A)/FI + C(F - A)/CFI 


(Eq. G-92) 


21. The ratio of thrust minus drag to the ideal thrust is 

(F - DVFI = C(F - D)/CFI 


(Eq. G-93) 


22. The ratio of thrust to the ideal thrust is 

F/FI = CFJC/CFI 


(Same as Eq. G-45) 


23. The resultant thrust coefficient about the yaw axis is 


CTY 


-fc 


csjc 2 +cfjc 2 


(Eq. G-94) 
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24. The thrust vector angle about the yaw axis is 


DELTAY = TAN 1 (CSJC/CFJC) 


(Eq. G-95) 


25. The total resultant thrust coefficient is 


CTT = \CNJC 2 + CFJC 2 + CSJC 2 


(Eq. G-96) 


26. The total resultant thrust ratio is 

FT/FI = CTT/CFI 


(Eq. G-97) 


27. The thrust minus nozzle axial force coefficient is 

C(F - ANOZ) = C(F - A) + (CDFAFT - CDFNOZ) + CFSUN) 

(Eq. G-98) 


28. The ratio of thrust minus nozzle axial force to the ideal thrust is 

F - ANOZ/FI = C(F - ANOZ)/CFI 


(Eq. G-99) 


29. The thrust minus nozzle drag coefficient is 

C(F - DNOZ) = C(F - D) + (CDFAFT - CDFNOZ) + CDSUN 


(Same as Eq. G-59) 


30. The nozzle drag coefficient is 




\ 


1 


V 
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CDNOZ = CDAERO - (CDFAFT - CDFNOZ) - CDSUN 

(Eq. G-100) 


31. The ratio of thrust minus nozzle drag to the ideal thrust is 

(F - DNOZyFI = C(F - DNOZyCFI 

(Same as Eq. G-49) 


32. The thrust coefficient (from balance and pressures) is 

CF = (F - ANOZ) + CDFNOZ + CFBUN 

(Eq. G-101) 

33. The ratio of thrust (from balance and pressures) to the ideal thrust is 


CF/CFI = CF/CFI 


(Eq.G-102) 


L. When IFAFN.= 1 

1. The thrust removed stability axis lift coefficient is 

CLSAER = CLAERO 


(Eq.G-103) 


t 






4 


2. The thrust removed stability axis drag coefficient is 

CDSAER = CDAERO 


(Eq. G-104) 
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3. The thrust removed stability axis side force coefficient is 


CSSAER = CSAERO 


(Eq.G-105) 


4. The thrust removed stability axis pitching moment coefficient is 


CMSAER = CMAERO 


(Eq.G-106) 


5. The thrust removed stability axis rolling moment coefficient is 

CRSAER = CRAERO * [COS(ALPHA)] + CYAERO * [SIN (ALPHA)] 

(Eq. G-107) 


6. The thrust removed stability axis yawing moment coefficient is 

CYSAER = CYAERO ♦ [COS(ALPHA)] - CRAERO * [SIN (ALPHA)] 

(Eq. G-108) 


7. The thrust removed wind axis drag coefficient is 

CDWAER = CDSAER * COS(BETA) - CSSAER * SEN (BETA) 

(Eq. G-109) 


8. The thrust removed wind axis side force coefficient is 

CSWAER - CSSAER * COS(BETA) + CDSAER ♦ SIN (BETA) 

(Eq. G-UO) 


230 


9. The thrust removed wind axis lift coefficient is 


CLWAER = CLSAER 


10. The thrust removed wind axis rolling moment coefficient is 


(Eq. G-lll) 


CRWAER = CRSAER * COS(BETA) + CMSAER * SIN (BETA) 


(Eq. G-112) 


11. The thrust removed wind axis pitching moment coefficient is 


CMWAER = CMSAER * COS(BETA) - CRSAER * SIN (BETA) 


(Eq. G-113) 


12. The thrust removed wind axis yawing moment coefficient is 


CYWAER = CYSAER 


1. The axial force coefficient is modified 


CA1 = CA1 - 0.0004 


2. The drag coefficient in the stability axis is modified 


CDS1 = CDS1 - 0.0004 


(Eq. G-114) 


(Eq. G-115) 


(Eq.G-116) 


h 
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3. The tlirust removed normal force coefficient is 

CNAERO = CN1 - CNJC 


(Same as Eq. G-77) 


4. The thrust removed axial force coefficient is 

CAAERO = CA1 + CFJC 


(Same as Eq. G-78) 


5. The thrust removed pitching moment coefficient is 

CMAERO = CMY1 - CMJC 


(Same as Eq. G-79) 


6. The thrust removed rolling moment coefficient is 

CRAERO = CMX1 - CRJC 


(Same as Eq. G-80) 


7. The thrust removed yawing moment coefficient is 


CYAERO = CMZ1 - CYJC 


(Same as Eq. G-81) 


8. The thrust removed side force coefficient is 

CSAERO = CY1 - CSJC 


(Same as Eq. G-82) 
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9. The thrust removed lift coefficient is 

CLAERO = CNAERO * COS(ALPHA) - CAAERO * SIN(ALPHA) 

(Same as Eq. G-83) 


10. The thrust removed drag coefficient is 

CDAERO = CAAERO * COS(ALPHA) + CNAERO ♦ SIN(ALPHA) 

(Same as Eq. G-84) 

11. The computed resultant thrust about the pitch axis is 


ct=Vcnjc 2 +cfjc 2 


(Same as Eq. G-34) 

12. The computed thrust vector angle about the pitch axis is 


DELTA = TAN 1 (CNJC/CFJC) 


(Same as Eq. G-35) 


13. The jet reaction lift coefficient is 


CLJET = CT [SIN (ALPHA + DELTA)] 


(Same as Eq. G-36) 


The jet reaction axial force coefficient is 

CFJET = CT [COS (ALPHA + DELTA)] 


(Same as Eq. G-37) 
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15. The thrust minus axial force coefficient is 


C (F - A) » -CA1 


(Same as Eq. G-41) 


16. The thrust minus drag coefficient is 

C(F - D) = -CDS1 


(Same as Eq. G-42) 


17. The ratio of thrust minus axial force to ideal thrust is 


(F - A)/FI = C(F - AyCFI 


(Same as Eq. G-43) 


18. The ratio of thrust minus drag to ideal thrust is 


(F - DyFI = C(F - DVCFI 


(Same as Eq.G-44) 


19. The ratio of thrust to ideal thrust is 

F/FI = CFJC/CFI 


(Same as Eq. G-45) 


20. The computed resultant thrust coefficient about the yaw axis is 


ct=^csjc 2 +cfjc 2 


(Same as Eq. G-94) 
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21. The computed thrust vector angle about the yaw axis is 


DELTAY = TAN 1 (CSJC/CFJC) 


22. The computed resultant thrust coefficient is 


CTT = yCNJC 2 + CFJC 2 + CSJC 2 


23. The total resultant thrust ratio is 


FT/FI = CTT/CFI 


24. The thrust minus nozzle drag coefficient is 

C (D - FNOZyFI = CDS2 - CDS1 


25. The ratio of thrust minus nozzle drag to ideal thrust is 

(F - FNOZyFI = C (D - FNOZ v CFI 


26. The nozzle lift plus thrust coefficient is 


(Same as Eq. G-95) 


(Same as Eq. G-96) 


(Same as Eq. G-97) 


(Eq. G-117) 


(Eq.G-118) 


CLNOZT = CLS1 - CLS2 


(Eq. G-119) 
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27. The nozzle pitching moment plus lift coefficient is 


CMNOZT * CMYS1 - CMYS2 


(Eq.G-120) 


28. The nozzle drag coefficient is 

CDNOZ = CDAERO - CDS2 


(Eq. G-121) 


29. The thrust removed nozzle lift coefficient is 

CLNOZ = CLAERO - CLS2 


(Eq. G-122) 


30. The thrust removed nozzle pitching moment coefficient is 


CMNOZ = CMAERO - CMYS2 


(Eq. G-123) 


N. Other Options 

1. ID - Engineer's option 

a. If ID s 1, the engineer may write his own option with the following 
restrictions: 

(1) Names must be identical to those already used. 


(2) No more terms may be added to the output. 
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MODULE H TURBOPROP OPTIONS 


SYMBOL 

NOMENCLATURE 

AD(J) 

Area at rake in exhaust duct J, sq. in. 

AE(J) 

Exhaust area for exhaust duct J, sq. in. 

ALPHAP(J) 

Angle of attack at propeller J, degrees. 

ARATIO(J) 

Area ratio for motor J. 

AT(J) 

Throat area of exhaust duct J, sq. in. 

CDP(J) 

Propeller drag coefficient for motor J. 

CDTP 

Total propeller drag coefficient. 

CE 

Exhaust sonic velocity, feet per second. 

CO 

Sonic velocity, feet per second. 

CHPROP(J) 

Chord length at 75 percent radius of propeller J. 

CMPROP(J) 

Pitching moment coefficient of propeller J. 

CNPROPCJ) 

Normal force coefficient of propeller J. 

CPPROP(J) 

Power coefficient of propeller J. 

CTPROP(J) 

Thrust coefficient of propeller J. 

DIAP(J) 

Diameter of propeller J, feet. 

ETA(J) 

Efficiency for motor J, per cent. 

ETAP(J) 

Efficiency for propeller J, per cent. 

FAPT 

Total system thrust in streamwise direction, lbs. 

FNPT 

Total propeller normal force coefficient. 

FTE(J) 

Propeller thrust plus jet thrust due to exhaust flow for motor 
J, lbs. 

FTGE(J) 

Total system thrust of motor J, lbs. 

JP(J) 

Advance ratio of propeller J. 
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SYMBOL NOMENCLATURE 


KPINM(I,J) 

KPOUTM(I,J) 

KPST(I,J) 

KPW 

KTINM(I.J) 

KTOUTM(I.J) 

MD(J) 

MEW) 

MTIP(J) 

MYPT 

NPROP(J) 

NSAME(J) 

PDRIVE 

PE(J) 

PHIANG 

PINM(I,J) 

PITCH(J) 

POUTM(I.J) 

PR/PTR 

PST(I.J) 


Constant for input drive pressure tap I and motor J, (must be 
0.0 or 1.0). 

Constant for output drive pressure tap I and motor J, (must be 
0.0 or 1.0). 

Constant for rake static pressure tap I and motor J. (must be 
0.0 or 1.0). 

Power coefficient constant. 

Constant for input motor temperature tap I and motor J, (must 
be 0.0 or 1.0). 

Constant for output motor temperature tap I and motor J, 
(must be 0.0 or 1.0). 

Rake mach number for motor J. 

Exhaust mach number for motor J. 

Mach number of propeller tip J. 

Total propeller pitching moment. 

Revolutions per second of propeller J. 

Constant of propeller J set equal to 0.0 or 1.0. 

Pressure drop through air turbine motor, lbs/sq. in. 

Exhaust static pressure for motor J, lbs/sq. in. 

Angle between forward and rotational velocities, degrees. 

Motor input static pressure for motor J and pressure tap I, 
lbs/sq. in. 

Measured value of geometric pitch of propeller J. 

Motor output static pressure for motor J and pressure tap I, 
lbs/sq. in. 

Ratio of static to total pressures. 

Static pressure for motor J and pressure tap I at rake, 
lbs/sq. in. 



e 

t 
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SYMBOL 


PSTATC(J) 

PW1(J) 

PW2(J) 

RHO 

RPS 

TDRIVE(J) 

TE(J) 

TINM(I.J) 

TOUTM(I.J) 

TSPROPCJ) 

TO 

VE(J) 

VIS 

VO 

VRES(J) 

VRN 


NOMENCLATURE 

Average static pressure at rake for motor J, lbs/sq. in. 

Horsepower output by motor J with ideal gas calculations, HP. 

Horsepower calculated using Isentropic equation multiplied by 
efficiency for motor J, HP. 

Density of free-stream air, slugs per cubic feet. 

Engine's revolutions per second. 

Temperature differential across the air turbine motor J, *F. 

Exhaust temperature for motor J, *F. 

Input temperature for motor J and temperature tap I, *F. 

Output temperature for motor J and temperature tap I, *F. 

Rotational tip speed of propeller J, feet per second. 

Tunnel static temperature, *R. 

Exhau st velocity in motor J, feet per second. 

Free-stream air viscosity based on tunnel air static 
temperature, lbs., secVsq. ft. 

Free-stream velocity, feet per second. 

Total velocity of propeller tip J, feet per second. 

Total velocity at 75 percent of propeller radius (for Reynolds 
number), feet per second. 
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APPENDIX H 
Module H 
Turboprop Options 


A. Introduction 

1. Module B with its constants must be run first. All constants are to be 
initialized to a value of zero. The project engineer must supply only those 
constants which are required for those quantities to be computed. In 
addition, by logical use of combinations of these constants, several options 
are available to the project engineer. 

2. Set NSAME(J) = 1 if POUTM(U) = PST(U), and TOUTM(I,J) = TTJ(U) 

where J = engine number 
I = probe number 

3. Set the constant, KTINM(I,J), equal to 1.0 for the temperature measuring 
probe. If the temperature probe is defective or does not exist, set the 
constant equal to 0.0. Use only a maximum of six probes per engine. 

4. The meaning of the values of KTOUTM(I,J) is the same as KTINM(I^J). 

5. Set the constant, KPINM(I,J) equal to 1.0 for the pressure measuring 
probe. If the probe is defective or does not exist, set the constant equal to 
0.0. The pressure probes may be weighted using this constant, if desired. 
Use only a maximum of 12 probes per engine. 

6. The meaning of the values of KPOUTM(I,J) is the same as KPINM(I^J). 
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7. Set the constant, KPST(I^J), equal to 1.0 for the pressure measuring probe. 
If the probe is defective or does not exist, set the constant equal to 0.0. The 
pressure probes may be weighted using this constant, if desired. Use only 
a m a xim u m of 12 probes per engine. 

8. AE(J) is equal to AT(J) for a converging nozzle. Both constants are 
required. Values of AT(J) come from Module B. 

B. Test for Air Turbine Simulator 

1. The constant required from the project engineer input at Module B is 
NUMENG (0 to 4). 

If NUMENG = 0, skip module H. 


C. Compute Common Consent 

1. The constants required from the proj ct engineer input at Module B are 
GAMJ and RJ. 


GAMJ+l 

f | ( 2 ^2(GAMJ-1) jGAMJ ^ 

IU1= l GAMJ + l J V RJ 


32.174 


(Same as Eq. B-l) 


KJ2 = 


GAMJ *64.348 
(GAMJ-l)RJ 


(Same as Eq. B-2) 


l 


KJ3 


■f, 


2(GAMJ) (RJ) 
(GAMJ -1)32.174 


(Same as Eq. B-3) 


„ 4 

I , 


4 


f 




D 
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KJ4 = 


GAMJ-1 

GAMJ 


(Same as Eq. B-4) 


KJ5 = 


1 

GAMJ 


(Same as Eq. B-5) 

2. To continue, the equations are given to show calculations for other 
constants. 

a. The static temperature is 


TO = (TTO + 459.67 V( 1.0 + 0.2 * MACH 2 ) 


(Eq. H-l) 


b. The free-stream density is 

RHO = PO * 144.0/(1716.4829 * TO) 


(Eq. H-2) 


c. The viscosity is 

VIS = 2.270 * 10~ 8 * TO * VTO/ (TO + 198.6) 


d. The free-stream velocity of sound is 

CO = 49.02 1 1 79 * VTO 


(Eq.H-3) 


(Eq.H-4) 
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e. The free- stream velocity is 


VO * CO * MACH 


(Eq. H-5) 

D. Individual Engine Mea surements 

1. This module provides the computations for four separate engines with the 
following instrumentation in each engine, 
a. Input pressure to engine 
*b. Output pressure of engine 
c. Input temperature to engine 
*d . Output temperature of engine 

e. Static exhaust pressure at rake 

f. Revolutions per second indicator 

g. Geometric pitch of propeller 


E. Propeller Coefficient Calculations 
1. The tip speed of propeller is 

TSPROP(J) = 3.14159 * DIAP(J) ♦ NPROP(J) 


(Eq. H-6) 


2. The advance ratio of propeller is 

JP(J) * VO/(NPROP(J) * DIAP(J)) 


(Eq. H-7) 


* May be replaced with rake measurements. 
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3. The angle of attack of the propeller ia the geometric pitch of the propeller at 
the 3/4 chord, in degrees, minus the resultant angle between the free- 
stream velocity and rotational velocity. 


ALPHAP(J) = PITCH(J) - PHIANG 


(Eq. H-8) 


where 


PHIANG = TAN' 1 (VO/VROT) 


(Eq. H-9) 


VROT = (3/4) TSPROP(J) 


(Eq. H-10) 


4. The Reynolds number for the propeller is calculated at the 3/4 chord. 


RNPROP(J) = VRN ♦ RHO ♦ CHPROP(JWIS 


(Eq. H-ll) 


where 


VRN = Resultant velocity at the 3/4 chord 


= J VROT 2 + VO 2 


VIS = Free-stream air viscosity calculated by Ames table equation, based on 
tunnel air static temperature 
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5. The Mach number of the propeller tip is 


MTIP(J) = VRES/CO 


where 


VRES = \V0 2 + TSPROP(J) 2 


6. Calculate the thrust coefficient of the propeller and hub using 
SCALE = RHO ♦ NPROP(J) 2 * DIAP(J) 4 


CTPROP(J) = FAREF1/SCALE 


where 


FAREF1 comes from Equation D-76 for J + 1 balances 


7. Calculate the normal force coefficient of the propeller and hub 


using 


CNPROP(J) = FNREF 1/SCALE 


(Eq. H-16) 


where 


FNREF 1 comes from Equation D-76 for J + 1 balances 


8. Calculated the pitching moment coefficient of the propeller and hub using 


S> 


CMPRORJ) = MYREF1ASCALE * DIAP(J) * 12.0) 


(Eq. H-17) 


where 


MYREF1 comes from Equation D-76 for J + 1 balances 


9. If NSAME(J) equals 1, then 

POUTM(I,J) = PST(I,J) 


(Eq. H-18) 


and 


TOUTM(I.J) = TTJ(I,J) 


(Eq. H-19) 


10. Calculations for the power coefficient of the propeller and hub are: 

a. Turbine inlet temperature 

TIN{ J) = I TTNM(I.J) ♦ KTINM(I, J) 

I KTINM(I.J) 


(Eq. H-20) 


b. Turbine outlet temperature 


TOUT(J) = 


I TOUTM(I, J) * KTOUTM(I,J) 
I KTOUTM(I.J) 


(Eq. H-21) 


c. Turbine inlet pressure 
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PIN(J) = £ PINMd J) ♦ KPINM(I, J) 
I KPINM(I, J) 


(Eq. H-22) 


d. Turbine outlet pressure 

POUT(J) = I POUTM(I,J) » KPOUTMd, J) 
X KOUTM(I,J) 


(Eq. H-23) 


e. The drive pressure across the air turbine engine is 


PDRIVE(J) = PIN(J) - POU1XJ) 


(Eq. H-24) 


f. The drive temperature across the air turbine engine is 


TDRIVE(J) = TIN(J) - TOUT(J) 


(Eq. H-26) 


g. The engine's revolutions per second are 


RPS = NPROP(J) / V(TIN(J) + 459.67)/ 518.7 

(Eq. H-26) 

h . Calculate the horsepower output from the air turbine engine u sing 

PW1(J) = (6006.0 * (WPENG(J)/32.174) * TDRIVE(J)/550) * 

((KPW13 * PIN(J) + KPW12 * RPS + KPW11) * RPS + KPW10) 

(Eq. H-27) 


PW2(J) * (6006.0 * (WPENG(jy32. 174) 
* (TIN(J) + 459.67) 
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* (1.0 - (POUT(jyPIN(J)) 2/7 ) 

* ETA(jy550 

(Eq. H- 28) 


where 

ETA(J) is determined by linear interpolation from a table. 

i. The power coefficient of the propeller and hub is 

CPPROP(J) = PW/(RHO * NPROP(J) 3 * DIAP(J) 5 ) 

PW = PW2(J) IF KPW = 0 
PW = PW1(J) IF KPW = 1 

(Eq. H-29) 


j. The propeller efficiency is 

ETAP(J) = CTPROP(J) * JP(J)/CPROP(J) 


(Eq.H-30) 


Exhaust Calculations 

1. Calculate the exhaust duct Mach number (rake position) using 


a. Duct static pressure 


PSTATKX J) = 


I PST(I,J)»KPST(I,J) 
I KPST(I,J) 


(Eq. H-31) 


b. The pressure ratio at the duct rake is 

PR/PTR(J) = PSTATIC(jyPTENG(J) 


(Eq. H-32) 
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c. If PR/PTRCJ) » PSTATIC(J) < 0.5283, use the Newton Raphson method 


for MD(J). 


MD(J) = . 



2 > 

5/7 

5 

7*MD(J) -1 

J PR T 

\6* 

6 

Ipto(j)J 

1 

J 



(Eq. H-33) 


d. If PR/PTR(J) > 0.5283, use this calculation of subsonic duct Mach 
numbers for MD(J). 

MD(J) = \5 • (PR / PTR(J))" 2/7 - 5 


(Eq. H-34) 


2. The ratio of A* to Area at the rake position and at the exit is 


a. Calculate A*/A at the rake station using 


ASTR / A = (1.728 * MD( J)) * 


-3 


1 + 


MD(J) 


(Eq. H-35) 


b. Calculate AVA of the exhaust exit using 

ARATIO(J) = ASTR/A ♦ AD(J)/AE(J) 

(Eq. H-36) 

3. Calculate the exhaust Mach number at the exit using an iteration 
technique on the formula 
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ME(J) = -^1 * (ARATIO(J)) * 1 + ^£L 
216 5 


4. The exhaust static temperature calculation is 


TE = (TTENG(J) + 459.67) • 1.0 + 


ME(J) 


where 


TTENCKJ) comes from Equation B-9. 


5. The exhaust sonic velocity is 


CE= 49.0211 79 *VTE 


6. The exhaust velocity is 


VE(J) = ME(J) * CE 


7. The exhaust static pressure is 


(Eq. H-37) 


(Eq. H-38) 


(Eq. H-39) 


(Eq. H-40) 


PE( J) = PTENG(J) • 1 + 


ME(J) 


(Eq. H-41) 
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8. The total propeller pitching moment is 


MYPT * I MYj 


(Eq. H-42) 


9. The total propeller normal force is 


NUMENG 

FNPT= I NF, 

1=1 

(Eq.H-43) 


10. The total propeller axial force is 


FAPT = I AF, 
1=1 


(Eq.H-44) 


11. The axial force coefficient in the body axis with propeller and jet thrust 
removed is 


CAPRS = CAAERO + FAPT 


(Eq.H-45) 


12. The drag coefficient in the stability axis with propeller and jet thrust 
removed is 


CDPRS = CAPRS [COS ( ALPHA)] + CNAERO {SIN (ALPHA)] 

(Eq.H-46) 

13. The side force coefficient in the stability axis with propeller and jet thrust 
removed is 
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CSPRS = CSAERO 


(Eq.H-47) 


14. The lift coefficient in the stability axis with propeller and jet thrust 
removed is 


CLPRS = CNAERO [COS (ALPHA)] - CAPRS [SIN (ALPHA)] 

(Eq.H-48) 

15. The rolling moment coefficient in the stability axis with propeller and jet 
thrust removed is 


CRPRS = CRAERO [COS (ALPHA)] + CYAERO [SIN (ALPHA)] 

(Eq. H-49) 

16. The pitching moment coefficient in the stability axis with propeller and jet 
exhaust removed is 





CMPRS = CMAERO 


(Eq. H-50) 


17. The yawing moment coefficient in the stability axis with propeller and jet 
exhaust thrust removed is 

CYPRS = CYAERO [COS (ALPHA)] - CRAERO [SIN (ALPHA)] 

(Eq. H-51) 

18. The drag coefficient in the wind axis with propeller and jet thrust removed 
is 
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CDPEW = CDPRS [COS (BETA)] - CSPRS [SIN (BETA)] 


19. The side force coefficient in the wind axis with propeller and jet thrust 
removed is 

CDPRW = CSPRS [COS (BETA)] + CDPRS [SIN (BETA)] 

(Eq. H-53) 

20. The lift coefficient in the wind axis with propeller and jet exhaust thrust 


removed is 


CLPRW = CLPRS 


(Eq. H-54) 


21. The rolling moment coefficient in the wind axis with propeller and jet 
exhaust thrust removed is 


CRPRW = CRPRS [COS (BETA)] + CMPRS [SIN (BETA)] 


(Eq. H-55) 


22. The pitching moment coefficient in the wind axis with propeller and jet 
exhaust thrust removed is 


CMPRW = CMPRS [COS (BETA)] - CRPRS [SIN (BETA)] 


(Eq. H-56) 


23. The yawing moment coefficient in the wind axis with propeller and jet 
exhaust thrust removed is 


CYPRW = CYPRS 


(Eq. H-57) 
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Inlet Distortion 

Nomenclature 

Required Constants 

Calculation of inlet weight flow/mass flow 

Engine face data 

Bleed flow computations 

Calculation of bypass weight flow/mass flow.. 

Total airflow 

Bad tube substitution scheme 
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SYMBOL 


c 

[ 


b 




L 


ADBPS(J) 

Inlet bypass flow controller steady state distortion, where J = 1 for 
left engine and J = 2 for right engine. 

ADPS(J) 

Inlet engine face static pressure distortion, where J = 1 for left 
engine and J s 2 for right engine. 

AFBSD(J) 

Engine bypass controller steady state distortion, where J = 1 for 
left engine and J = 2 for right engine. 

AFSD(J) 

Engine face steady state distortion index, where J = 1 for left 
engine and J = 2 for right engine. 

AMFD(J) 

Inlet mass flow plug choked flow effective area, sq. in., where 
J = 1 for left engine and J = 2 for right engine. 

AMFDB(J) 

Engine bypass controller plug choked flow effective area, sq in., 
where J = 1 for left engine and J = 2 for right engine. 

APBR(J) 

Engine bypass controller plug effective area, sq. in., where J = left 
engine and J = 2 for right engine. 

APR(J) 

Mass flow plug effective area, sq. in., where J = 1 for left engine 
and J = 2 for right engine. 

ARXT(JXN) 

Engine area of extent, sq. in., where J = 1 is left engine, J = 2 is 
right engine, and N is the ring number. 

AS(J) 

Engine face Mach number, where J = 1 for left engine and J = 2 
for right engine. 

ASB(J) 

Bypass instrumentation plane Mach number, where J - 1 for left 
engine and J = 2 for right engine. 

CAEFT 

Total engine face axial force coefficient, CAEF1+CAEF2. 

CAEF(J) 

Engine face axial force coefficient, where J = 1 is left engine J = 2 
is right engine. 
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SYMBOL NOMENCLATURE 


C(M) 


CB(M) 


CE(M) 


C(IXM) 

DCI(JXN) 

DPRS(J) 

DRI(JXN) 

EXT(JXN) 

KC(N) 

KR(N) 


Constants to be used in mass flow/weight flow calculation. Where 
M is the constant number, M=1 to 25 for the left engine; M=26 to 50 
for the right engine; M=51 to 56 choked flow for the left engine; 
M=57 to 62 choked flow for the right engine. 

Constants to be used for bleed weight flow calculation. Five for 
each bleed with a maximum of 12 bleeds per engine, where M = 1 
to 60 for the left engine and 61 to 120 for the right engine. 

Constants to be used in mass flow/weight flow calculation for a 
bypass controller. Where M is the constant number, M=1 to 25 for 
the left engine; M=26 to 50 for the right engine; M=51 to 56 choked 
flow for the left engine; M=57 to 62 choked flow for the right 
engine. 

Constants for computing bleed mass flow, where I is bleed 
number and M is a constant number of 1 to 5 for each bleed. 

Inlet circumferential distortion intensity, where J = 1 is left 
engine, J = 2 is right engine, and N is the ring number 
Engine loss in surge pressure ratio, where J = 1 for left engine 
and J = 2 for right engine. 

Engine radial distortion intensity, where J = 1 for left engine, J = 2 
for right engine, and N is ring number. 

Engine theta value of extent, where J = 1 for left engine, J s 2 for 
right engine, and N is ring number. 

Engine circumferential distortion constant, where N is ring 
number. 

Engine radial distortion constant, where N is ring number. 
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SYMBOL 

K(N) 

KERAK(L) 

| KEXIT(L) 

KPPBSCL) 
KPPBTCL) 
KPPT(L) 
KPRMSfL) 
KPSB(L) 
KPTB(L) 
KPTRK(L) 
KTERK(L) 

: / 

r KTTB(JXL) 


NOMENCLATURE 

Engine constant, where N is ring number. 

Constant used to compute engine face static pressure, where L is 
tap number; default is 0.0. 

Constant used to compute flow plug exit static pressures, where L 
is tap number; default is 0.0. 

Constant used to compute bypass flow plug exit static pressure, 
where L is tap number; default is 0.0. 

Constant used to compute bypass flow plug exit total pressure, 
where L is probe number; default is 0.0. 

Constant used to compute flow plug exit total pressure, where L is 
probe number; default is 0.0. 

Constant used to compute flow RMS pressure, where L is probe 
number; default is 0.0. 

Constant used to compute bleed flow exit static pressure, where L 
is tap number, default is 0.0. 

Constant used to compute bleed flow exit total pressure, where L 
is probe number; default is 0.0. 

Constant used to compute bypass instrumentation plane total 
pressure, where L is probe number; default is 0.0. 

Constant used to compute engine face total pressure, where L is 
probe number; default is 0.0. 

Constant used to compute bypass flow plug exit total temperature, 
where L is probe number; default is 0.0 : where J = 1 for left 
engine and J = 2 for right engine. 



j) 
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SYMBOL 



KTTP(L) 

MFLB 

MFLO 

M/MOB(J)(I) 

M/MOBKJ) 

M/MOBP(J) 

M/MOI(J) 

M/MOP(J) 

M/MOT(J) 

MPR(J)(N) 

NBLD(J) 

NBS(JXI) 


Constant used to compute flow plug exit total temperature, where 
L is probe number; default is 0.0. 

Bypass controller flow calculation flag, where 0 = no bypass flow, 
1 = mass flow computation, and 2 = weight flow computation. 
Inlet flow calculation flag, where 0 = weight flow and 1 = mass 
flow. 

Bleed mass flow ratio, where J = 1 for left engine, J = 2 for right 
engine, and I is bleed number. 

Bypass flow plug mass flow ratio based on instrumentation plane 
pressure ratio, where J = 1 for left engine and J = 2 for right 
engine. 

Bypass flow plug nozzle mass flow ratio, where J * 1 for left 
engine and J = 2 for right engine. 

Engine mass flow ratio based on engine face pressure ratio, 
where J = 1 for left engine and J = 2 for right engine. 

Inlet mass flow plug nozzle mass flow ratio, where J = 1 for left 
engine and J = 2 for right engine. 

Inlet total mass flow ratio, where J = 1 for left engine and J = 2 for 
right engine. 

Engine multiples per revolution, where J = 1 for left engine, J = 2 
for right engine, and N is ring number. 

Number of bleeds in engine, where J = 1 for left engine and J = 2 
for right engine. 

Number of static pressures in engine bleed, where J = 1 for left 
engine, J = 2 for right engine, and I is the bleed number. 


5 ? 
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SYMBOL 



NOMENCLATURE 


f 


r 

p 





NBT(JXI) 

NPPBSKJ) 

NPPBT(J) 

NPPS(J) 

NPPT(J) 

NPRING 

NPRMS(J) 

NPSEF(J) 

NPTB(J) 

NPTEF(J) 


Number of total pressures in engine bleed, where J = 1 for left 
engine, J = 2 for right engine, I is the bleed number. 

Number of bypass flow plug exit static pressures, where J = 1 for 
left engine and J = 2 for right engine. 

Number of bypass flow plug exit total pressures, where J = 1 for 
left engine and J = 2 for right engine. 

Number of exit static pressures on left mass flow plug rake, 
where J = 1 for left engine and J = 2 for right engine. 
NPPS1+NPPS2 maximum of 20 

Number of total pressures on mass flow plug rake, where J = 1 for 
left engine and J = 2 for right engine. NPPTl+NPPr2 maximum 
of 100. 

Number of total pressures per ring at engine face. 

Number of engine face RMS pressures, where J = 1 for left engine 
and J = 2 for right engine. NPRMS1+NPRMS2 max of 80. 

Number of static pressures on engine face rake, where J = 1 for 
left engine and J = 2 for right engine. NPSEF1+NPSEF2 
maximum of 80. 

Number of bypass instrumentation plane total pressures, where 
J s 1 for left engine and J = 2 for right engine. NPTB1+NPTB2 
maximum of 80. 

Number of total pressures on engine face rake, where J = 1 for left 
engine and J = 2 for right engine. NPTEF1+NPTEF2 maximum of 
100 . 


- 1 

• 4 

1 

* 
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SYMBOL NOMENCLATURE 


NPTOP 

PAFSD(J) 

PALWR(JXN) 

PAVR(JKN) 

PB(JXL) 

PB(JV 

PTB(JXI) 

PEFRMSA(J) 

PEFRMS<L) 

PERAKE(L) 

PEXBR(J) 


Number of total temperature probes at mass flow plug, where 
J = 1 for left engine and J = 2 for right engine. NPTT1+NPTT2 
maximum of 12. 

Mass flow plug steady state distortion index, where J = 1 for left 
engine and J = 2 for right engine. 

Average engine face low pressure ratio used to compute 
distortion, where J = 1 for left engine, J = 2 for right engine, and N 
is ring number. 

Engine face average total pressure per ring, lbs./sq. in., where 
J = 1 for left engine, J = 2 for right engine, and N is ring number. 
Engine bleed static pressure, lbs./sa in., where J = 1 for left 
engine, J = 2 for right engine, and L is the tap number. 

Ratio of average bleed static pressure to average bleed total 
pressure, where J = 1 for left engine, J = 2 for right engine, and I 
is the bleed number. 

Engine face average RMS pressure, lbs./sq. in., where J = 1 for left 
engine and J = 2 for right engine. 

Engine face RMS pressure, lbs./sq. in., where L is probe number 
with a maximum of 40. 

Engine face rake static pressure, lbs./sq. in., where L is the tap 
number. 

Average bypass flow plug exit static pressure divided by average 
flow plug exit total pressure, where J = 1 for left engine and J = 2 
for right engine. 
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SYMBOL 


NOMENCLATURE 




PEXIT(L) 

PEXTR(J) 

PFSDB(J) 

PI(J)/PTO 

PPBS(J) 

PPBS(L) 

PPBT(J) 

PPBT(L) 

PPS(J)/PTO 

PPT(L) 


Mass flow plug exit static pressure, lbsVsq. in., where L is tap 
number. 

Average flow plug exit static pressure divided by average flow 
plug exit total pressure, where J = 1 for left engine and J = 2 for 
right engine. 

Bypass flow plug steady state distortion index, where J = 1 for left 
engine and J = 2 for right engine. 

Average engine face static pressure divided by freestream total 
pressure, where J = 1 for left engine and J = 2 for right engine. 
Ratio of average bypass flow plug exit static pressure to free 
stream total pressure, where J = 1 for left engine and J = 2 for 
right engine. 

Bypass flow plug exit static pressure, lbsVsq. in., where L is tap 
number. 

Ratio of average bypass flow exit total pressure to free stream total 
pressure, where J = 1 for left engine and J = 2 for right engine. 
Bypass flow plug exit total pressure, lbsVsq. in. where L is probe 
number. 

Average mass flow plug static pressure divided by free stre am 
total pressure, where J = 1 for left engine and J = 2 for right 
engine. 

Mass flow plug exit total pressure, lbs./sq. in., where L is the 
probe number. 
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SYMBOL NOMENCLATURE 


PPTCJVPTO 

PRMS(L) 

PRMS(LV 

PTI(J) 

PSAB(JXI) 

PSB(L) 

PSEKJyPTO 

PTAB(JXI) 

PTB(L) 

PTBD(J) 

PTERAK(L) 


Average mass flow plug exit total pressure divided by freestream 
total pressure, where J = 1 for left engine and J = 2 for right 
engine. 

Engine face RMS pressure, lbs./sq. in., where L is tap number 
Engine face RMS pressure divided by engine face average total 
pressure, where J = 1 for left engine, J = 2 for right engine, and L 
is tap number 

Ratio of average engine bleed static pressure to free stream total 
pressure, where J = 1 for left engine, J = 2 for right engine, and I 
is bleed number. 

Engine bleed static pressure, lbs./sq. in., where L is the tap 
number. 

Average flow plug static pressure divided by free stream total 
pressure, where J = 1 for left engine and J = 2 for right engine. 
Ratio of average engine bleed total pressure to free stream total 
pressure, where J = 1 for left engine, J = 2 for right engine, and I 
is bleed number. 

Engine bleed total pressure, lbs./sq. in., where L is the probe 
number. 

Ratio of average bypass instrumentation plane total pressure to 
free stream total pressure, where J * 1 for left engine and J = 2 for 
right engine. 

Engine face total pressure, lbs./sq. in., where L is the probe 
number. 


SYMBOL NOMENCLATURE 


PTKJVPTO 

PTRAKB(L) 

SAB(JXI) 

SANF(J) 

SCAP(J) 

SANFB(J) 

SCAPB(J) 

SEF(J) 

SEFB(J) 

TFRMS(J) 

TTBP(JXL) 

TTBP(J) 


Average engine face total pressure divided by freestream total 
pressure, where J = 1 for left engine and J = 2 for right engine. 
Bypass instrumentation plane total pressure, lbs7sq. in., where L 
is probe number. 

Engine bleed exit area, sq. in., where J = 1 for left engine, J = 2 for 
right engine, and I is bleed number. 

Inlet engine face annular flow area, sq. in., where J = 1 for left 
engine and J * 2 for right engine. 

Inlet capture area, sq. in., where J s 1 for left engine and J * 2 for 
right engine. 

Bypass annular flow area, sq. in., where J * 1 for left engine and 
J = 2 for right engine. 

Bypass capture area, sq. in., where J = 1 for left engine and J = 2 
for right engine. 

Engine face area, sq. in., where J = 1 for left engine and J = 2 for 
right engine. 

Bypass instrumentation plane area, sq. in., where J = 1 for left 
engine and J = 2 for right engine. 

Engine face RMS turbulence, where J = 1 for left engine and J = 2 
for right engine. 

Bypass flow plug total temperature, *F, where J = 1 for left engine, 
J s 2 for right engine, and L is probe number. 

Average bypass flow plug total temperature, *F, where J * 1 for 
left engine and J = 2 for right engine. 
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SYMBOL 


NOMENCLATURE 


TTB(J) 

TTE(J) 

TTP(J) 

TTP(L) 

WCBR'JXI) 

WCB(JXI) 

WCCR(J) 

WCCBR(J) 

WCP(J) 

WCR(J) 

WCT(J) 

WC/WCC(J) 


Bypass instrumentation plane computed total temperature, *R, 
where J = 1 for left engine and J = 2 for right engine. 

Engine face computed total temperature, *R, where J = 1 for left 
engine and J = 2 for right engine. 

Mass flow plug total temperature, *F, where J = 1 for left engine 
and J = 2 for right engine. 

Mass flow plug total temperature, *F, where L is the probe 
number with a maximum of 8. 

Inlet bleed corrected mass flow, lbs/sec, where J = 1 for left 
engine, J = 2 for right engine, and I is bleed number. 

Inlet bleed corrected mass flow, lbs/sec, where J s 1 for left 
engine, J = 2 for right engine, and I is bleed number. 

Engine face choked weight flow, lbs/sec, where J = 1 for left engine 
and J = 2 for right engine. 

Bypass instrumentation plane choked weight flow, lbs/sec, where 
J = 1 for left engine and J = 2 for right engine. 

Engine weight flow based on engine face pressure ratio, lbs/sec, 
where J = 1 for left engine and J = 2 for right engine, where J = 1 
for left engine, J = 2 for right engine, and I is bleed number. 
Engine face corrected weight flow, lbs/sec, where J = 1 for left 
engine and J * 2 for right engine. 

Engine face total corrected weight flow, lbs/sec, where J = 1 for left 
engine and J s 2 for right engine. 

Engine corrected weight flow ratio based on flow plug calibration, 
where J = 1 for left engine and J = 2 for right engine. 
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SYMBOL 

NOMENCLATURE 

WPB(JXI) 

Engine bleed physical weight flow, lbs/sec, where J - 1 for left 
engine, J = 2 for right engine, and I is bleed number. 

WPT(J) 

Engine face total physical airflow, lbs/sec, where J = 1 for left 
engine and J = 2 for right engine. 

WP(J) 

Engine face physical air flow, lbs/sec, where J = 1 for left engine 
and J = 2 for right engine. 

XMBE(JXI) 

Engine bleed exit Mach number, where J = 1 for left engine, J = 2 
for right engine, and I is bleed number. 

XMEF(J) 

Engine face Mach number, where J = 1 for left engine and J = 2 
for right engine. 

XMFB(J) 

Bypass instrumentation plane Mach number, where J s 1 for left 
engine and J = 2 for right engine. 

XMFP(J) 

Mass flow plug nozzle mass flow function, where J = 1 for left 
engine and J - 2 for right engine. 

XMMBP(J) 

Bypass flow plug nozzle Mach number, where J = 1 for left engine 
and J = 2 for right engine. 

XMP(J) 

Mass flow plug nozzle Mach number, where J = 1 for left engine 
and J = 2 for right engine. 

XPBR(.J) 

Bypass flow plug axial position, in., where J = 1 for left engine and 
J = 2 for right engine. 

XPR(J) 

Mass flow plug axial position, in., where J * 1 for left engine and 
J = 2 for right engine. 
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AmiMU 

Module 1 
Inlet Distortion 


A. Required 

The constants for inlet distortion calculations are given in the nomenclatures. 
Constants of the same name are also described by modules B and E. All 
constants are initialized to a value of 0.0. 



1. I RAKE - Rake code 

where IRAKE * 5, inlet distortion using rotating rake 

IRAKE = 6, inlet distortion using nonrotating rake 
If IRAKE s 0, 1-4 skip module I 


l 


B - Calculation of Inlet Weight Flow/Maas Flnw 

The same general equation can be used for the calculation of either an area 
from which mass flow is computed or for the direct calculation of weight flow. 
The difference will be in the way the calculated terms are used in succeeding 
calculations. The calculation path will be determined by the flag MFLO. 

APR(J) or WCR(J) = Cl + C2*XPR(J) + C3*PEXTR(J) + C4*XPR(J) 2 + 

C5*XPR(J)*PEXTR(J) + C6*PEXTR(J) 2 + C7*XPR(J) 3 + 
C8*XPR(J) 2 *PEXTR(J) + C9*XPR(J)+ PEXTR(J) 2 + 
C10*PEXTR(J) 3 + C11*XPR(J) 4 + C12*XPR(J) 3 *PEXTR(J) 
+ C13*XPR(J) 2 *PEXTR(J) 2 + C14*XPR(J)*PEXTR(J) 3 + 
C15*PEXTR(J) 4 + C16*XPR(J) 5 + C17*XPR(J) 4 *PEXTR(J) 
+ C18*XPR(J) 3 *PEXTR(J) 2 + C19*XPR(J) 2 *PEXTR(J) 3 + 
C20*XPR(J) 3 *PEXTR(J) 3 + C21*XPR(J)*PEXTR(J) 4 + 
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C22*PEXTR(J) 5 + C23*XPR(J) 5 *PEXTR(J) + 
C24*XPR(J)*PEXTR(J) 5 + C25*XPR(J) 5 * PEXTR(J) 5 



(Gq, 1-1) \\ 

i 


t 

Calculation for choked flow: 


AMFD(J) OR WCCR(J) = C51 + C52*XPR(J) + C53*XPR(J) 2 + C54*XPR(J) 3 + i 

C55*XPR(J) 4 + C56 *XPR(J) 5 


(Eq. 1-2) 


where PEXTR(J) = PPS(J)/PTO/PPT(J)/PTO 

(Eq. 1-3) 

tv , L =NPPS(J) L-NPPS<J) 

PPS( J) / PTO = I PEXIT(L) * KEXIT(L) / XKEXIT(L)*PTO 

and L=i l=i 


(Eq. 1-4) 


If MFLO is equal to 1 the area terms APR(<> ) and AMFD(J) are calculated and 
the following equations are used. 


If PEXTR(J) is ^ 0.6 then use the choked flow value to calculate a mass flow 
ratio as follows: 


L=NPTT(J) L=NPT T(J) 

TTP(J)= X TTP(L) * KTTP(L) / XKTTP(L) 

L-l L=1 

M/MOP(J) = (AMFD(J)*PPT(J)/PTO)/(XMO*SCAP(J)) 

where XMO=0.9189*MACH*(1+0.2*MACH 2 )' 3 
and MACH is obtained from Appendix A 


(Eq. 1-5) 

(Eq. 1-6) 
(Eq. 1-7) 
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If PEXTR( J)>0.6 then this equation is used: 


M/MOP(J) = (APR(J)*XMFPNL*PPT(J)/PTO)/(XMFFO* 

((TTP(J) + 459.67W(TTO + 459.67) *SCAP(J)) 

(Eq. 1-8) 


where XMPL = V5*(PEXTK(J)* 2/7 -l) 

and XMFPNL = 0.9189*XMPL*(1+0.2*XMPL 2 )' 3 

(Eq. 1-9) 


Average flow plug exit total pressure ratio 

L=NPPT(J) L=NPPT(J) 

PPT( J) / PTO = z PPT(L) * KPPT(L) / IKPPT(L) 

L=1 L=1 

(Eq. MO) 

Mass flow plug steady state distortion index 

PAFSD(J) = ((PPT(L)/PTO) max - (PPT(L)/PTO) min )/PPT(J)/PTO 

(Eq. I- 11) 

Engine face Mach number. 

AS(J) = (XMO*(SANF(J)/M/MOP(J))*PTI(J)/PTO)/SCAP(J) 

(Eq. 1-12) 

where XMO is defined in Eq. 1-7. 

If AS(J) is £ 0.53177 then XMEF(J) is 1.0. Otherwise compute engine face 
Mach number by iteration from the following 


(Eq. M3) 


AS(J) = 0.9189 *XMEF(J)*(1+0.2*XMEF(J) 2 )* 3 


If MFLO is equal to 0 (the default) the WC values are calculated and the 
following equations are used. 

WCAVCC(J) = WCR(J)AVCCR(J) 

(Eq. 1-14) 


WP(J) = (WCR(J)*(PTI(J)/PTO)*PTO)/(VTrO+459.67 *RC) 

(Eq. 1-15) 

where RC= 14.696/V518.68 


wliich is a reference pressure divided by a reference temperature and is a 
constant used oecause of the flow plug calibration. 


WCP(J) = MA*SEF*RC 


(Eq. 1-16) 


where MA = 0.9189*M2*(1+0.2*M2 2 )' 3 

and M2 = Vs^O'TKJi/PTO/PKJJ/PTO)' 277 -1) 


M/MOP(J) = (WCR(J)*PTI(J)/PTO)/(XMO*SCAP(J)*RC) 


where XMO is defined in Eq. 1-7. 


(Eq. 1-17) 


M/MOI(J) = (WCP(J)*PTI(J)/PTO)/(XMO*SCAP(J)*RC) 

(Eq. 1-18) 
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Engine face Mach number 


EFMAC(J) = WCR(J)/(RC*SANF(J)) 


(Eq. 1-19) 


If EFMAC(J) £ 0.53177 set XMEF(J) equal to 1.0. If EFMAC(J) is < 0.53177 
then compute XMEF(J) by iteration from the following expression 

EFMAC(J) = 0.9189 *XMEF(J)*(1+0.2*XMEF 2 )* 3 

(Eq 1-20) 


Computed engine face total temperature 

TTE(J) = (TT0+459.67V(1+0.2*XMEF(J) 2 ) 


(Eq. 1-21) 


C. Engine face data 

Average engine face total pressure ratio 


L=NPTEF(J) L=NPTEF(J) 

PTI(J)/PTO= I PTERAK(L) * KTERK(L) / IKTERK(L)* PTO 

L=1 L=1 


(Eq. 1-22) 


Average engine face static pressure ratio 


L=NPSEF(J) 

PI(J) / PTO = £ PERAKE(L) * KERAK(L) / 

L=1 


L=NPSEF(J) 

£ KERAK(L) * PTO 


L=1 


(eq. 1-23) 


RMS pressure ratio 
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[ 





PRMS(LVPTI(J) * PRMS(L)/PTI(jyPTO*PTO 


Engine face RMS turbulence 


(Eq. 1-24) 


L=NPRMS(J) L=NPRMS(J) 

TFRMS(J) = X PRMS(L) * KPRMS(L) / ZKPRMS(L)*PTO 

L=1 L=1 


(Eq. 1-25) 


Engine face static pressure distortion 

ADPS(J) = ((PERAKE(L)/PTO) max - (PERAKE(L)/PTO) min )/PI(J)/PTO 

(Eq. 1-26) 


Engine face total pressure distortion 

AFSD(J) = ((PTERAK(L)/PTO) max - (PTERAK(L)/PTO) min )/PTI(J)/PTO 

(Eq. 1-27) 

Aerodynamic interface plane (AIP) theta extents computed by the 
ARP 1420 methodology. 5 Using these areas and theta values the 
following parameters are computed: 

Circumferential distortion by ring 

DCI(JXN) = (PAVR(JXN) - PALWR(JXN)VPAVR(J)(N) 

(Eq. 1-28) 

where: 







PALWR(JXN) = PAVR(J)(N) - (ARXT(JXN)/EXT(J)(N)) 


(Eq. 1-29) 


5 Gas Turbine Engine Inlet Flow Distortion Guidelines, ARP Paper 1420, Society of Automotive 
Engineers. Inc. 
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Multiples per revolution by ring 




L-NUMEXT 

MPR(JXN)= ' X AREA(I) / AREAMAX 
1=1 

(Eq. 1-30) 

Radial distortion by ring 

DRI(JXN) = ((PTI(J)/PTO)*PTO - PAVR(JXN))/(PTI(JVPTO)*PTO 

(Eq. 1-31) 

If the inlet data is for a specific engine, then the engine constants may be 
input to compute the loss in surge pressure ratio 

N=NRING(J) 

DPRS(J) = X(KC(N) * DCI(JXN) + KR(N) * DRI( JXN) + K(N)] 

N=1 

(Eq. 1-32) 

Engine face stream forces 


CAEF(J) = SANF(J)/(QO*SAREA(I))*((PTI(J)/PTO*PTO 

*(l+GAMJ*XMEF(J) 2 ))/((l+0.2*XMEF(J) 2 ) 2/7 )-PO 

(Eq. 1-33) 


where SAREA(I) is obtained from Appendix D. 



(Eq. 1-34) 


) 


D. Bleed Flow Computations 

All bleed flow computations are to be made for a maximum of 12 separate 
bleeds in each engine. 

Average bleed total pressure 

L=NBT(JXI) L=NBT(JXI) 

PTAB(JXI) = I PTB(L) * KPTB(L) / IKPTB(L) * PTO 

L=1 L=1 

(Eq. 1-35) 



Average bleed static pressure divided by freestream total pressure 


L=NBS(JXI) L=NBS(JXI) 

PSAB(JXI) = I PSB(L) * KPSB(L) / IKPSB(L)*PT0 

L=1 L=1 


(Eq. 1-36) 


Ratio of average bleed static to average bleed total pressure 


PB(JXI)/PTB(JXI) = PSAB(JXI)/PTAB(J)(I) 


(Eq. 1-37) 


Bleed corrected mass flow from calibration 


WCB(JXI) = CB1(M) + CB2(M)*PB(J)(I)/PTB(J)(I) + CB3(M)* 
(PB(JXiyPTB(J)d)) 2 + CB4(M)*(PB(J)(I)/PTB(J)(I)) 3 + 
CB5(M)*(PB(JXI)/PTB)JXl)) 4 

(Eq. 1-38) 

Bleed exit weight flow corrected to engine face 


WCPB(JXI) = (WCB(JXI)*RC*V TTO+459.67 )/PTI(J)/PTO*PTO 

(Eq. 1-39) 
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Bleed mass flow ratio 




M/MOBCJXI) = (RC*WCB(JXI^PTAB(JXI)*PTOVMACH*SCAP(J) 

(Eq. 1-40) 


1 



where MACH is obtained from Appendix A 
Bleed exit Mach number 

XMBE(JXI) = V 5*((l/PB(JXiyPTB(J)(I))' 2/7 -D 




(Eq. 1-41) 


If exit area is provided instead of a calibration then weight flow is computed in 
the following manner. 

WCB(JXI) = (MB*SAB(JXD*PTAB(JXI)*PTOWTTO+459.67 

(Eq. 1-42) 


where MB = 0.9189 *XMBE(JXI)*(1+0.2*XMBE(JXI) 2 )' 3 


(Eq. 1-43) 


Weight flow corrected to engine face 

t 

WCPB(JXI) = WCB(JXI)*RC*V'ITO+459l67/PTI(J)/PTO*PTO 

(Eq. 1-44) 


Bleed mass flow ratio is then computed 


M/MOB(JXI) = (WCB(J)(lWTTO+459.67)/(SCAP(J)*MO*PTO) 

(Eq. 1-45) 
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E. Calculation of Bypass Weight Flow/Mass Flow 

If there is a separate controlled bypass or ejector flow from the inlet, many of 
the same values must be calculated for the bypass controller as for the inlet 
flow controller. The same general equation can be used for the calculation of 
either an area from which mass flow is computed or for the direct calculation 
of weight flow. The difference will be in the way the calculated terms are used 
in succeeding calculations. The calculation path will be determined by the 
flag MFLB. 

APBR(J) = CE1 + CE2*XPBR(J) + CE3*PEXBR(J) + CE4*XPBR(J) 2 + 
or CE5*XPBR(J)*PEXBR(J) + CE6*PEXBR(J) 2 + CE7*XPBR1 3 

WCBR(J) + CE8*XPBR(J) 2 *PEXBR(J) + CE9*XPBR(J)* PEXBR(J) 2 

+ CE10*PEXBR(J) 3 + CE11*XPBR(J) 4 + CE12*XPBR(J) 3 
*PEXBR(J)+CE13*XPBR(J) 2 *PEXBR(J) 2 +CE14*XPBR(J) 
*PEXBR(J) 3 +CE15*PEXBR(J) 4 + CE16*XPBR(J) 5 
+ CE17*XPBR(J) 4 *PEXBR(J) + CE18*XPBR(J) 3 
*PEXBR(J) 2 + CE19*XPBR(J) 2 *PEXBR(J) 3 
+ CE20*XPBR(J) 3 *PEXBR(J) 3 *CE21*XPBR(J)*PEXBR(J) 4 
+ CE22*PEXBR(J) 5 + CE23*XPBR(J) 5 *PEXBR(J) 

+ CE24*XPBR(J)*PEXBR(J) 5 + CE25*XPBR(J) 5 * PEXBR(J) 5 

(Eq. 1-46) 


Calculation for choked flow: 

AMFDB(J) = CE51 + CE52*XPBR(J) + CE53*XPBR(J) 2 

or + C54*XPBR(J) 3 + CE55*XPBR(J) 4 + CE56 *XPBR(J) 5 

WCCBR(J) 

(Eq. 1-47) 
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where PEXBRGJ) = PPBS(jyPPBT(J) 


W 


(Eq. 1-48) 


L=NPPBS(J) L=NPPBS(J) 

PSBD( J) = I PPBS(L) * KPPBS(L) / IIO>PBS(L)*PTO 

L=1 L=1 


(Eq. 1-49) 


and 


PPBT(J) = 


L=NPPBT(J) 

I 

L=1 


L=NPPBT(J) 

PPBT(L) * KPPBT(L) / 1 KPPBT(L) * PTO 


L=1 


(Eq. 1-50) 

If MFLB is equal to 2 the area terms APBR(J) and AMFDB(J) are calculated 
and the following equations are used. 


If PEXBR(J) is £ 0.6 then use the choked flow value to calculate a mass flow 
ratio as follows: 


L=NPTTB(J) L=NPTT B(J) 

TTBP(J) = I TTBP(J)(L)*KTTB(JXL)/ IKITB(JXL) 

L=1 L=1 

(Eq. 1-51) 

M/MOBP(J) =(AMFDB(J)*PPBT(J)/PTO)/(XMO*SCAP(J)) 

(Eq. 1-52) 

where XMO is defined in Eq. 1-7. 

If PEXBR(J)>0.6 then this equation is used: 

M/MOBP(J) = (AP(J)*XMFPBNL*PPBT(J)/PTO)/(XMO* 

((VTTPB(J)+459.67 'VTTO+459.67)*SCAP(J)) 

(Eq. 1-53) 
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where XMPBL = V5*(PEXBR(J)' 2/7 -l) 

(Eq. 1-54) 

and XMFPBNL = 0.9189*XMPBL*(1+0.2*XMPBL 2 )' 3 

(Eq. 1-55) 

Average bypass flow plug total pressure ratio 

L=NPPBT(J) L=NPPBT(J) 

PTBD(J) = £ PTRAKB(L) * KPTRK (L) / XKPTRK(L)*PTO 

L=1 L=1 

(Eq. 1-56) 

Bypass mass flow plug steady state distortion index 

PFSDB(J) = (PPBT(LVPTO) max - (PPBT(L)/PTO) min /PPBT(J)/PTO 

(Eq. 1-57) 

Bypass instrumentation plane Mach number. 

ASB(J) = (XMO*(SANF(J)/M/MOP(J))*PTBD(J)/P)/SCAPB(J) 

(Eq. 1-58) 


where XMO is defined in Eq. 1-7. 


(Eq. 1-59) 


If ASB(J) is £ 0.53177 then XMBF(J) is 1.0. Otherwise compute 
instrumentation plane Mach number by iteration from the following 


ASB(J) = 0.9189 *XMFB(J)*(1+0.2*XMFB(J) 2 )' 3 

(Eq. 1-60) 

If MFLO is equal to 0 (the default) the WC values are calculated and the 
following equations are used. 
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WC/WCCB(J) * WCBR(J)AVCCBR(J) 


(Eq. 1-61) 

WBP(J) = (WCBR( J)*(PTBD( J)/PTO)*PTO)/(' / TTO+459.67 *RC) 

(Eq. 1-62) 


where RC= 14.696^518.68 

which is a reference pressure divided by a reference temperature and is a 
constant used because of the flow plug calibration. 


WCPB(J) = MA*SEFB*RC 


where MA *= 0.9189*M3*(1+0.2*M3 2 )' 3 

and M3 = V 5*((PTBD(J VPPBS(J)) 2/7 -1) 


(Eq. 1-63) 


M/MOBP(J) = (WCBR(J)*PTBD(J)/PTO)/(XMO*SCAPB(J)*RC) 

(Eq. 1-64) 


where XMO is defined in Eq. 1-7. 


M/MOBKJ) = (WCPB(J)*PTBD(J)/PTO)/(MO*SCAPB(J)*RC) 


Bypass instrumentation plane Mach number 

XMFB(J) = WCBR(J)/(RC*SANFB(J)) 


(Eq. 1-65) 


(Eq. 1-66) 

If EFMACB(J) 2 0.53177 set XMFB(J) equal to 1.0. If EFMACB(J) is < 0.53177 
then compute XMFB(J) by iteration from the following expression 


EFMACB(J) = 0.9189 *XMFB(J)*(1+0.2*XMFB(J) 2 )' 3 

(Eq. 1-67) 
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TTB(J) = (TT0+459.67y(l+0.2*XMFB(Jr) 


(Eq. 1-68) 


Total corrected mass flow ratio 


I=NBLD(J) 

M / MOT(J) = M / MOP(J) + £M / MOi3( J)(I) + M / MOBP(J) 


(Eq. 1-69) 


Total corrected weight flow 


I=NBLD(J) 

WCT(J) = WCR(J) + £ WCPB(JXI) + WCBR(J) 


(Eq. 1-70) 


Total physical weight flow at engine face 


I=NBLD(J) 

WPT(J) = WP(J) + £ WCB(J)(I) + WBP(J) 


(Eq. 1-71) 


The calculations that use the engine face total pressures require that there be 
values for all pressures on the rake. Since it is an unusual occurrence to get 
through a test without losing some pressures, a probe substitution scheme 
has been added to the data reduction. This system works from a ring rake 
probe numbering scheme as shown in figure 1-1. All engine face total 
pressures must be numbered and arranged in this order. Although the 
figure shows a 40 probe arrangement with five rings and eight rakes, it will 
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also work for a different arrangement of more or less rings with more or less 
probes as long as all the rings have the same number of pressures. In the 
following explanation 'n' will denote ring number and 'm' will denote rake 
number. The KTERK values are used to activate the probe substitution. When 
a KTERK value is zero, the program determines where the bad probe is 
located and substitutes for it by the following criteria. All probe substitutions 
are made before any values using the total pressures are computed. 


1. One bad probe on an interior ring (rings 2, 3, or 4). This value will be taken as 
an average of the four surrounding values, two from the same ring, two from 
the same rake. 


(n,m) 


^• P (n.m+1) P (n.m-1) * P (n+l,m) P (n-l,m)^ 


1 


(Eq. 1-72) 


2. One bad probe on an I. D. or O. b. ring (rings 1 or 5). This value will be 
computed as an average of the three surrounding values. 

p F P (l.m+l) + P (l.m-1) * P (2.m) j 

r (l,m> ~ 3 

p _t P (5,m+l) + P (5,m-1) + P (4,m)J 
r (5,m) " 3 


3. Two adjacent probes bad on an interior ring. Values for this case are 

computed using the average of the three surrounding values, two from the 
same rake, and one from the same ring. 
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_*■ P (n+l.m) + P (n-l.m) + P (n,n>-1)^ 

*(n,m) ” 3 

p £ P (n+l.m+l) * P (n-l.m+l) 4 P (n,m+2p 

* (n,m+l) “ 3 

(Eq. 1-73) 


4. Two adjacent probes bad on an I. D. or O. D. ring. These values are computed 
from the average of the two surrounding values, one from the same rake, and 


one from the same ring. 

p [P (l,m-l) + P (2,m) J 

r (l,m) " 2 

J P (l.m+2) 4 P (2.m+1) ] 

Ml,m+1) — 2 

p ^ P (5.m-1) 4 P (4,m) } 

M5,m) " 2 

p Jjj5jn+2) 4 P (4.m-H ) } 

N5,m+1) ~ 2 

(Eq. 1-74) 

5. Two adjacent interior probes bad on the same rake. These values are 

computed as the average of two probes from the same ring, and one from the 


same rake. 

p P (n,m+1) 4 P (n.m-1) 4 P (n-l.m) \ 

*(n,m) “ 3 

p _} P (n+l.m+l) 4 P (n+l.Tn-l) 4 P (n+2,m) j 

*(n+l,m) “ 3 


(Eq. 1-75) 
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6. Two adjacent probes bad on the same rake, with one being on an I D. or O. D. 
ring. 

^ [P (5.m+l) + P (5,m-1) ] 

p [ P (4.m+1) * P (4 t m-l)‘ f P (3,m) } 

*(4.rn) “ 3 

p _t P d.m,l) + P (l,m-1) j 

r (l.m) “ 2 

p J P (2,nH-l) + P (2.m-1)+ P (3.m) j 
^(2,m) " 3 

(Eq. 1-76) 
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